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PAPERS 


TRANSPORTATION SUSPENDED SEDIMENT 
WATER 


Measurements the sediment and velocity distributions laboratory 
flume were made for various values rate flow, slope channel, and size 
and amount suspended load. The experimental sediment-distribution 
measurements were compared with the theoretical distribution derived the 
assumption that the coefficient for the turbulent transfer suspended sediment 
the same the coefficient for the turbulent transfer momentum. 

The following are the main findings: 

(1) The measured sediment distributions have the same form the 
theoretical distribution but not agree quantitatively with them. This 
shows how invalid the basic assumption that the transfer coefficient for 
sediment was equal that for momentum. general, the measured relative 
concentrations were larger than those given the theory, the difference 
being greatest for the fine sediment and diminishing the suspended load 
became coarser. 

(2) Suspended sediment flow tends reduce appreciably the 
transfer momentum and hence the resistance flow, allowing the sediment- 
laden water flow more rapidly than comparable clear water flow. This 
effect varies with the total sediment load and the grain size the sediment. 


(3) Suspended load flow may factor causing the secondary 
circulation observed rivers. 


INTRODUCTION 


Many important practical problems today involve the movement 
sediments water, and their complete solution requires more knowledge 


comments are invited for immediate publication; insure publication the last dis- 
cussion should submitted November 1944 


Project Supervisor, Cooperative Laboratory, Dept. Agriculture, SCS, Asst. Prof. Hy- 
draulics, California Inst. Technology, Pasadena, Calif. 
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this subject than now exists. Flood control and channel maintenance 
streams, both large and small, are just much problems sediment movement 
water flow. Established practices, fact, are successful handling 
the water flow floods, but the problem passing the sediment load still 
unsolved. recent years the work the United States Department 
Agriculture developing methods reducing accelerated erosion agri- 
cultural lands has encountered many new sediment-transportation problems, 
which have shown the urgent need for additional fundamental knowledge 
this subject. connection with this need that this study was undertaken. 

Sediment may transported flowing water two essentially different 
ways—by rolling sliding along the bed the stream channel suspension 
the body the fluid. The material transported called, respectively, 
the load the stream,” simply load,” and the “suspended 
load.” Just above the bed continual interchange material occurring 
between the bed and the overlying fluid and this locality obviously difficult 
distinguish between the bed load and suspended load. The two types 
transportation are means independent. They are separated only for 
convenience studying them and them.? This paper deals only 
with suspended load. 


Previous 


theories suspended-load transportation are based 
developments the mechanics turbulent flow fluids. Consider two- 
dimensional turbulent flow the z-direction with velocity u(y), which 
normal Then, according Osborne the shear parallel 


which the mass density the fluid, and and are the turbulent 
velocity fluctuations the and directions, respectively, and the bar denotes 
mean value. 

Boussinesq‘ introduced expression for the shear the form, 


which the quantity (often called the eddy viscosity) analogous 
the coefficient viscosity the expression for viscous shear; and the 
mean velocity the z-direction. The quantity has the dimensions 
kinematic viscosity but unlike its counterpart function flow con- 
ditions and, therefore, varies from point point. The quantity 
cient expressing the exchange momentum between neighboring filaments. 


Distinction Between Bed-Load and Suspended Load Natural Streams,’’ Einstein, 
Anderson, and Johnson, Transactions, Am. Geophysical Union, Pt. II, pp. 628-633. 


the Dynamical Theory Incompressible Viscous Fluids and the Determination the Cri- 
Osborne Reynolds, Catalogue Scientific Papers (1800-1900), compiled the Royal Soc. 
London, Cambridge Univ. Press, Vol. pp. 535-577; also Philosophical Transactions, Royal Soc. 
Vol. 186, 1895, 123. 


sur Theorie des Eaux Boussinesq, Mémoires, Presente par divers 
Savants, France, Vol. 23, pp. 1-680. 


| 
ea 
| 
; 
> 
| 
i 
q 
3 


June, 1944 SUSPENDED SEDIMENT 795 


This quantity called the turbulent momentum transfer coefficient. Obviously 
has the dimensions, length velocity. 
Theodor von Am. Soc. E., introduced the equation, 


which correlation coefficient, the average the absolute values 
the fluctuations normal the main flow, and the mixing length. From 
Eq. seen that the turbulent momentum transfer coefficient, has 
the value, 


and the dimensions, velocity length. 

Just Eq. expresses the rate transfer momentum due turbulent 
mixing, possible write similar equation the transfer suspended 
particles, any other property that carried the fluid. Thus, the rate 
transfer mass suspended particles per unit area 


which the concentration suspended material mass per unit volume. 
Professor von pointed out that the values v’, and 
and need not necessarily the same. 

For steady conditions the upward transport material due turbulence 
balanced the settling due the force gravity. This expressed 


which now the concentration the fraction the sediment with 
settling velocity Transposing, the differential equation for suspended 
sediment becomes, 


which the sediment transfer coefficient that similar but not 
necessarily equal it. Solution Eq. gives, 


constant over the depth, that is, constant, Eq. may 
integrated give 


Aspects the Turbulence Theodor von Proceedings, 4th International 
Cong. Applied Mechanics, 1934, pp. 54-91. 
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which the base the natural logarithms; and, simplify typography, 


Eq. was apparently first given Wilhelm who applied 
the study suspension dust particles the atmosphere. Eq. was derived 
more direct the work Professor von 

can integrated. Since the expression for not known the assumption 
made that and Eq. integrated this basis. From and 


which was first presented Professor von 1934. 
uniform open-channel flow whose ratio width depth large, 
these relations will apply. The expression for the shear stress given 


Introducing Eq. 12a into Eq. 11, 


loge 


expression for the velocity gradient given the von 
universal velocity defect law, 


U — Umax 


p 
Massenaustauch freier Luft und verwandte Wilhelm Schmidt, Probleme 
der kosmischen Physik, 1925, Pt. 


the Theory Turbulent Flow and Its Relation Sediment-Transportation,” 
O’Brien, Transactions, Am. Geophysical Union, Section Hydrology, 1933, pp. 


and Skin Theodor von Journal the Aeronautical Sciences, 
Vol. No. 1934 pp. 1-20. 


Introducing Eq. into Eq. yields the relation 
(12a) 
which the shear the bottom, or, 
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(0.40 for pipes) and the distance from the wall the pipe radius ro. 
For wide open channels, discussed here, Eq. replaced with ym, 
the depth the flow. The derivative the velocity then becomes, 


which, when introduced into Eq. 13, gives 


Integrating Eq. 


which, simplify typography, 


From possible plot curves relative sediment distribution 
was first done about 1936 Arthur Ippen, Assoc. Am. Soc. E., 
the suggestion Professor von 

Eq. 16, well the equations preceding it, gives only relative concen- 
tration and therefore does not make possible the calculation the total 
transport the stream. make this possible value for needed, and 
the problem obtaining the transporting capacity flow may simply 
stated one determining this value.. discussing the magnitude 
Professor von states that probably depends the size the 
sediment and the magnitude the shearing stress acting the ground. 
Hunter has suggested that the limiting value Eq. made 
equal the height the roughness elements further suggests that 
the corresponding concentration for this level calculated from Eq. 
that the concentration set equal that the bed material. 
Lane, Am. E., and Kalinske, Assoc. Am. Soc. E., 


356, Vereines Deutscher Ingenieure, Vol. 1932, 36. 


Rauhen Nikuradse, Forschungsheft 361, Vereines Deutscher 
Ingenieure, Vol. 1933, pp. 


Turbulent Flow Open Channels,” Garbis Keulegan, Journal Research, National 
Bureau Standards, Dept. Commerce, Vol. 21, 1938, pp. 


Conceptions the Mechanics Fluid Turbulence,” Hunter Rouse, Transactions, 
Am. Soe. E., Vol. 102 (1937), Fig. 30, 536. 


Analysis Sediment Transportation the Light Fluid Turbulence,” Hunter Rouse, 
Mimeograph Publication 1939. 
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approaching the problem statistically, that only those particles having 
settling velocities less than the vertical fluctuations the bed velocity 
will picked up. assuming normal distribution and setting the 
root mean-square fluctuation proportional the friction velocity 


they proceed calculate the probable rate suspension material. 


result analysis the energy balance stream carrying sus- 
pended load, Robert Am. Soc. E., advanced the hypothesis 
that the capacity stream limited only the settling velocity the 
sediment greater than the rate fall the stream which the 
average velocity flow and the slope the channel. Furthermore, 
deduced that when the settling velocity less than the material contributes 
energy the flow, thus making possible pick additional load until 
extremely high concentrations are reached. 

Verification Theories—Although comprehensive suspended-load 
experiments have been made which the important variables have been 
investigated systematically, some noteworthy experimental work has been done. 

Mr. appears have first applied the turbulence theory the 
problem suspended load his studies the suspension dust particles 
air. obtained Eq. which assumes constant value for the exchange 
coefficient. 

1929 Jun. Am. Soc. E., reported series experiments 
the suspension sediment cylindrical column water agitated uni- 
formly series propellers much the same those used Joule his 
experiments the mechanical equivalent heat. various levels his 
apparatus measured the concentrations three well-graded sands (mean 
size 0.9 mm, 0.4 mm, and 0.2 mm) for several propeller speeds and found that 
the distribution sediment followed exponential relationship the form 
Eq. Mr. Hurst was also able derive his relationship from analysis 
based the kinetic theory gases. 

Professor performed some experiments similar those Mr. 
Hurst using agitator coarse mesh bars that was oscillated the vertical 
direction with simple harmonic motion. The results Professor Rouse’s 
work indicate that Eq. valid, and verify the theory. Professor Rouse 
observed systematic deviation the value obtained for the coarse sand 
from that obtained from fine sand. The deviation was attributed differ- 
ence the correlation factor Eq. for mixing the fluid and the 
sediment. 

1932 John outlined method calculating the exchange 
coefficient for natural streams and discussed its effect suspending sediment 


Relation Suspended Bed Material Rivers,’’ Kalinske and Lane, Transac- 
tions, Am. Geophysical Union, Pt. 1939, pp. 637-641. 

1938, pp. 

Vol. 124, 1929, pp. 196-201. 

the Mechanics Sediment Suspension,” Hunter Rouse, Proceedings, 5th 
International Cong. for Applied Mechanics, 1938, pp. 

Theory the Morphologic Significance Turbulence the Flow Water Streams, 
John Leighly, University California Publications Geography, Vol. No. 1932, pp. 1-22. 
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the body the flow. paper Mr. calculated the exchange 
coefficient from some sediment measurements natural streams using 
Eq. also calculated the sediment distribution from the exchange 
coefficient step integration. 

1935 Assoc. Am. Soc. E., attempted check 
Eq. applying sediment measurements made Harry Blaney 
and the late Samuel Fortier, Members, Am. Soc. E., the Colorado River 
and Imperial general the observations agreed with Eq. but 
the results were not entirely conclusive since complete data were not available 
and was necessary estimate such quantities the slope the canals and 
the settling velocity the sediment before calculations could made. 

made some simultaneous sediment and velocity measure- 
ments flume wide, deep, and 150 long which led him 
conclude that the vertical distribution sediment nearly exponential. 
This agrees with Eq. and would indicate that the mixing coefficient had 
approximately constant value over the depth. Such condition could 
approximated channel downstream from disturbance such might 
occur the inlet flume. later this same author reports the 
concludes that the boundary region the concentration suspended sediment 
inversely proportional and that the open stream varies exponentially, 
much expressed Eq. 


APPARATUS AND PROCEDURE 


Main Flume.—The experiments were performed flume 33} in. wide 


long, arranged the slope and discharge could varied. The slope 
was adjusted means four pairs screw jacks, and the discharge was 
varied from per sec about 5.5 per sec changing the speed 
the pump the system through variable-ratio transmission the electric 
motor drive. The sides the channel were made 12-in. structural channels 
lined with rubber. The bottom was steel plate artificially 
roughened with sand. making the circulation system closed circuit with 
the minimum velocity the flume, was possible recirculate the suspended 
sediment with the water without danger depositing material the piping 
system. 

The flume first constructed was described considerable detail* 
Robert Knapp and Robert Oaks and shown diagrammatically 
Fig. Experiments with this circuit yielded fair results but indicated that 


and the Transportation Rock Debris Streams,” John Leighly, The Geo- 
graphical Review, Vol. 24, No. 1934, pp. 453-464. 


Silt Open Channels,” Christiansen, Transactions, Am. Geophysical 
Union, Pt. 1935, pp. 


Technical Bulletin No. 67, U.S.D.A., 1928. 


Science, Vol. 17, No. 1934, pp. 769-783. 


Suspension Solids Turbulent Stream,” Richardson, Proceedings, Royal Soc. 
London, Series Vol. 162, 1937, pp. 


‘Closed-Circuit’ Flume for Suspended-Load Studies,’’ Robert Oaks and Robert Knapp, 
Transactions, Am. Geophysical Union, Pt. 1936, pp. 525-528. 
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the disturbances due the pump and inlet were exerting appreciable 
effect flow the flume. 

The flow was improved revising the circuit shown Figs. 1(b) and 
Moving the pump and venturi meter the discharge end the flume provided 
the flow before entered the transition the 
flume section. eliminate the possibil- 
ity separation the included angle between the 
sides the transition was made 5°. 

Fig. shows general view the flume from 
the inlet end. The beam spanning the channel 
about midlength carries the instruments used 
making measurements. The diffuser lid the 
flume inlet shown the left. 

Adjustment flow the flume 
was readily adjusted given depth pump 
ing the correct volume water into the system 
and varying the discharge until uniform 
was reached. Uniform flow for given dis- 
charge was obtained easily (once the flow was set with the water manometer 
the venturi meter) adjusting the amount water the closed system. 
difficulty was ever experienced with the closed system because pendula- 
tion any other instability. The flow was adjusted easily and remained 
steady for indefinite lengths time. 

Depth depths were measured with point gage 
which can read 0.001 means vernier scale. Depth measure- 
ments were made first adjusting the elevation the bracket holding the 
instrument that the scale read zero when the point was resting the 
bottom and then lifting the point that just came contact with the 
surface. The average depth for the channel was obtained from measurements 

Velocity were measured with pitot static tube 
(diameter, in.) the standard Prandtl design. The differential pressure 
the tube was read 0.001 water manometer and converted velocity 
assuming coefficient unity for the tube. (That this assumption 

Sediment distribution sediment was determined from 
samples siphoned from the flume through brass-tube pipette (outside diam- 
eter, in.) shaped much like pitot tube. The tip the sampler was 
flattened its inside dimensions were 0.490 in. 0.139 in. The average 
velocity which the sediment-laden fluid entered the flat-tip section (sampling 
velocity) was made equal the stream velocity the sampling point 
adjusting the head the siphon. 

The time remove sample was calculated for each point from the 
measured stream velocity and the area the flat-tip sampler (0.064 in.), 


Peters, Handbuch der Experimentalphysik, Leipzig, Vol. Pt. 1931, 
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the head the siphon required for the rate flow was determined trial, 
and three 1-liter samples were taken. The sediment was removed from the 
sample filtering and was dried and weighed. The average weight the 
three samples was taken the concentration the point. Sample weights 
more than were determined cg, and the smaller quantities were 
determined mg. 

Water removed from the flume sampling was replaced order keep 
the depth uniform. was found that removing adding little gal 
water made measurable difference the tailwater depth. 


The effect the sampling velocity the amount sediment removed with 
the sample was studied experimentally. Several 1-liter samples were removed 
from point sampling velocities equal 1.00, 0.86, and 0.75 times the 
local stream velocity the sampling point. The sampling point was 0.5 
from the center the flume, from the inlet (that is, station 44), and 
0.1 from the bottom. The local velocity was 3.19 per sec, the flume was 
set slope 0.0025, and the water was flowing equilibrium depth 
0.495 ft. The suspended load had mean settling rate still water 0.058 
per sec. 
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The results this experiment, presented Table show that, the 
sampling velocity decreases, the quantity sediment collected per given 
volume fluid increases. When the sampling velocity lower than the local 
stream velocity, the fluid stream lines curve away from the sampler tip. 


TABLE DETERMINATION THE EFFECT THE SAMPLING 
VELOCITY THE AMOUNT SEDIMENT COLLECTED 


1.03 1.08 1.12 
1.00 1.05 1.09 the ratio sampling velocity the 


velocity the stream the sampling point; 


0.98 1.13 the average sediment concentration, grams 
1.04 per liter, for the series samples; the 

.06 0.99 0.97 stream velocity; and C/C the ratio the sedi- 
ment concentration for individual sample 

1.10 


Because their greater inertia the sediment particles follow paths less 
curvature than the fluid, with the result that sediment transferred the 
fluid entering the sampler and the apparent concentration increased. When 
the sampling velocity greater than the stream velocity the curvature 
toward the sampler, and tendency decrease the amount sediment 
collected would therefore expected. 

The values C/C line Table give some idea the fluctuation 
concentration point the flow. Although the deviation the weights 
individual samples from the mean was high 10%, interesting 
note that the deviation the average weight three consecutive samples 
from the mean does not exceed 3%. Furthermore, the length filament 
fluid having volume liter and area equal that the sampler tip 
well excess the length the flume, which ft. 

Sediment Used Suspended used suspended load were 
prepared grading silica sand wind-tunnel classifier 
George Otto and Professor Rouse. Fig. 4(a) shows logarithmic probability 
chart the size distributions the three sands used suspended load. 
The parameters that describe size and grading these sands are listed 
Fig. 4(b). 

The quantity the geometric mean sieve diameter and read directly 
from the chart. For normal distribution (that is, when the curve straight 
line), the geometric mean and median diameters are equal. The quantity 
(Fig. known the geometric standard deviation the sieve diameters. 
the ratio the geometric mean size the size corresponding the 
ordinate value 15.9%, which turn also equal the ratio the size corre- 
sponding the ordinate value 84.1% Dm. The method plotting the 


for Sand and George Otto and Hunter Rouse, Civil Engineering, 
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size distribution and calculating these parameters was shown Mr. 
The quantity weighted mean sedimentation diameter, the sedi- 
mentation diameter particle being defined the diameter sphere 
having the same density and the same settling velocity the particle. The 
mean settling velocity the 
particles was determined from 
observed settling times indi- 
vidual particles represent- 
ative sample the sediment. 
The mean sedimentation diam- 
eter then the diameter 
the quartz sphere having the 
settling velocity water. 
Calculation the sedimentation 
diameter was based the re- 
sistance coefficients for spheres 
The mean sedimentation di- 
ameter was also determined 
from sieve analyses taking 
weighted mean the sedimen- 
tation diameters the material 
held each sieve, which were 
determined from 
ling-velocity measurements 
samples taken from the sieve 
0.01 

0.061 0.074 0.208 0,246 The average sediment con- 
centration for the entire flow 
was determined from samples 
siphoned from the collection 
tank through pipe and 

Geometric Geometric Mean rubber hose. The entrance 

cally upward the downward 

0.147 0.160 flow the tank could enter with- 

0.091 0.100 out changing direction. During 

sampling the pipe was moved 

around obtain fluid from all 

parts the cross section. 
Bottom mentioned previously, the sides the flume were 
lined with rubber, and the bottom was roughened with The two bottom 


Modified Logarithmic Probability Graph for the Interpretation Mechanical Analyses 
George Otto, Journal Sedimentary Petrology, Vol. No. 1939, pp. 62-76. 
for the Settling Velocity Hunter Rouse, Exhibit the Annual 


Report the Committee Sedimentation, National Research Council, Div. Geology and Geography, 
Washington, C., 1936-1937, pp. 57-64. 
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roughnesses used the experiments were designated sands and The 
mean diameters for sands and were, respectively, 0.47 and 0.88 mm, 
and the corresponding standard deviations, were 1.12 and 1.10. These 
sands were graded the wind-tunnel classifier described previously. The 
sand was stuck the bed the flume with water-resistant bitumastic paint. 
The metal surface was first painted and the sand immediately applied 
sufficient quantity cover the surface one layer deep. Any excess material 
was swept off after the paint was dry. 


General Outline experiments this study may 
divided into two series. the first series, the flume circuit shown 
Fig. 1(a), the slope the flume 0.0025) was kept constant and only. 
one bottom roughness (sand and one suspended load (sand were used. 


Only the rate which water was discharged and the amount sediment 


the system were varied. mentioned under the heading, “Apparatus and 
Procedure,” this circuit gave disturbed condition the entrance and was 
not considered entirely satisfactory. The disturbances were due the pump, 
which was near the inlet the flume; and separation the flow the 
transition from the pipe the flume. The latter separation the flow was 
remedied, least part, baffles, which, course, introduced considerable 
resistance and, therefore, turbulence. good index the performance the 
system given the head loss the head developed the pump, which 

The second series experiments was made the revised circuit shown 
Fig. separation occurred the transition section, and the velocity 
distribution was uniform and symmetrical. Only about head was 
required circulate flow per sec the revised system. this 
series, complete measurements were made with one bottom roughness (sand 5), 
three sizes suspended load, and two depths each two slopes. 

Velocity for the distribution sediment over 
vertical section stream was derived for flow with large width-depth 
ratio—that is, two-dimensional conditions—and for the case logarithmic 
velocity distribution. Before the relation can applied the present 
experiments necessary show that the assumed conditions obtain. Fig. 
5(a), with Table 2(a), shows the velocity distribution across the flume for the 
first series experiments (runs velocity distribution curves 
for series (runs 22) are shown Fig. with Table 2(b). Comparison 
once the improvement the flow effected series changing the system 
from that shown Fig. the one Fig. 1(6). velocity distribu- 
tions also show that the influence the walls the channel does not extend 
the center the flume and that the assumption two-dimensional condi- 
tions near the center the flume justified. The only way the resistance 
the walls can transferred the center the channel turbulent shear 
momentum transfer, which can exist only there horizontal velocity 
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Distance from Feet 


Distance from Center, in Feet 


2 3 4 52 3 4 53 4 52 3 4 5 
Velocity, wu, in Feet per Second 
(a) SERIES 
Station 407 
2 3 42 3 4 51 2 3 41 2 3 


Velocity, u, in Feet per Second 


SERIES 
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gradient. The fact that this gradient zero very small permits the con- 
clusion that near the center the flume the influence the walls negligible. 

Fig. 5(a), distributions and Fig. distributions give 
the velocity various distances from the inlet, and show the gradual stabili- 
zation the velocity. Here again may noted the improved flow conditions 
series (Fig. 5(b)). Despite the higher discharge, the equilibrium velocity 
established short distance and obtains over the remaining length the 
flume. The distortion the velocity depth 0.5 station (Fig. 
distribution due backwater from the discharge end, which station 
57. Fig. also shows the effect the sediment the flow, discussed 
subsequently. 


Distribu- Depth Flow Velocity Concentration 

A 1 46 0.481 4.80 3.60 0 

B 1 46 0.484 4.75 3.54 0.75 
Cc 13 40 0.473 4.87 3.72 3.2 

F 1 46 0.485 4.85 3.61 0 

G 1 15 0.480 4.87 3.66 0.75 
H 1 26 0.480 4.87 3.66 0.75 
0.480 4.87 3.66 0.75 

B 14 40 
D 14 50 
538 5.24 3.51 1.3 

H 14 50 

K 18 45 0.462 2.82 2.20 0 

L 17 45 0.269 1.31 1.76 0.05 
M 20 45 0.462 4.07 3.18 1.2 
N 21 45 0.234 1.42 2.19 1.2 
1) 18 45 0.462 3.11 2.43 1.2 
P 22 45 0.295 2.13 2.60 1.3 


Figs. and show dimensionless semilogarithmic graphs velocity profiles 


for series and respectively. The relative distance from the bottom 


plotted the logarithmic ordinate scale against the dimensionless velocity 


the abscissa scale. The quantity the so-called friction 


velocity, commonly written scale shifted for each curve, 


the abscissa value being indicated for each curve. The curves Fig. 
are arranged four groups according the average sediment content the 
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Those Fig. are grouped according size the sediment being 
transported. 

all runs the velocity distribution fits the logarithmic law very well, 
although near the bottom there systematic deviation from the law, which 
not understood. departure expected near the bottom since the viscous 
forces are longer negligible and the conditions for which the logarithmic law 
was derived longer obtain. The viscosity tends suppress the turbulence 
and reduce the exchange coefficient Eq. 1b, and, order develop 
the shear, the velocity gradient must increase. This would cause the velocity 
decrease more rapidly than indicated the logarithmic law, and the 
result opposite the results shown Figs. and the present experi- 
ments, the origin the depth coordinate was assumed the top the 
sand grains forming the bottom roughness. inspection the velocity 
profile curves will show that the fit could improved assuming the origin 
below this level although the justification for doing not clear. The 
deviation from the law was first attributed errors measurement due 
interference between the pitot static tube (diameter, in.) and the bottom. 
Tubes small in. diameter gave substantially the same readings 
the larger instruments and showed that the measurements were reliable. 
apparent deviation from the law also noted the free surface. How- 
ever, because surface waves and the disturbance the pitot tube itself, 
felt that these measurements are not too reliable and that conclusions 
based them must drawn with reservation. 

The fact that the velocity profiles fit the logarithmic law well means 


that the average velocity will always occur the relative 


average velocities determined this method agree within with 
those determined graphical integration. The velocities Table were 
determined this method. 

The velocity curves Figs. and are the form, 


which universal constant. The quantity 2.3/k the slope the 
curve the semilogarithmic graph, and should constant does not vary. 
close inspection the curves will show that the slopes vary, the size 
and concentration sediment are varied. The values for each curve are 
shown Figs. and andin Table general, decreases the concen- 
tration increased the size sediment suspension decreased. 

Fig. shows the velocity profiles for series and respectively, plotted 
according the von velocity defect law (Eq. 14). Fig. 8(a), 
the data Fig. have been combined according the sediment content 
the flow and show clearly the variation Fig. the data Fig. 


Distribution Open Channels,” Vito Vanoni, Civil Engineering, June, 
356-357. 
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have been grouped according the values The curves show that 
decreased the presence the sediment, but the relation between size 
sediment and the universal constant not clear. Fig. does show, however, 
that the logarithmic velocity law valid and that its use justified de- 
veloping Eq. for suspended load distribution. 


Values of Log H 


0.01 

0.006 
0.004 
0.003 


Average Suspended Load, in Grams per Liter ream 


Relative Sediment Concentration, 


Distribution Suspended plotted logarithmically straight 
line with slope, 


which may read directly from the graph. Figs. and show the sediment 
distribution data for series and respectively, plotted logarithmically; 
that is, log plotted against log The origin the abscissa scale 
(log H), which occurs shifted for each curve. (The results 
experiment series were first reported the writer 


Experiments the Transportation Suspended Load,” Vito Vanoni, Transactions, 
Am. Geophysical Union, Pt. III, 1941, pp. 608-621 
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The value may also calculated from Eq. 17, which the settling 
velocity the sediment, the von universal constant, and 


the friction velocity. (In calculating the value the settling velocity 
the center mass group particles representative the suspended 
load was used.) distinguish between the value obtained Eq. and 
that obtained fitting the experimental data Figs. and 10, the latter 


value denoted Obviously, the theory correct, the two 


will agree. 

The intersection each curve with the line log 1.0 indicated 
short horizontal dash. The vertical distance from the origin 
the dash, measured cycles the logarithmic scale the value Values 
both and are listed Table The two values agree well for runs 
17, which were all made with the coarsest suspended load (sand 1). For 
all the runs with finer sand the value from about 20% 40% greater 
than 

From inspection Figs. and once apparent that the results 
series are much closer agreement with the distribution law (Eq. 16) 
than those series The disagreement occurs the upper levels the 
flow—that is, for small values H—where more sediment found than 
would expected from the law. Results the two series may compared 
only for those runs made under comparable conditions with sediment the 
same size. Runs series were made with the same sediment the 
runs series and may therefore compared with runs that series. These 
runs (14 17) also show more sediment the upper levels than given 
Eq. 16, but the difference slight. The rather large departures from the 
distribution law series are doubt due disturbances the flow the 
entrance, which, stated previously, were appreciable compared with those 
occurring the revised circulation system (Fig. used series 

Figs. and show that the departures the measured relative concen- 
trations from the distribution law are small, and the curves also show con- 
clusively that least the form the distribution function correct. This 
conclusion confirmed Fig. 11, which shows plotted against 
log for series and Here again the disagreement the upper levels 
between measurement and the distribution function noticeable, also 
the better agreement the results series 

Fig. shows the relative sediment distribution function the relative 
depth measured upward from the reference level 0.05 ym. The measure- 
ments are indicated circles and the theoretical distribution shown 
solid line. The agreement for series (Fig. 12(b)), previously seen, good 
over the entire depth, whereas, series (Fig. 12(a)) shows the departure 
the higher levels noted before. 

well note that, calculating the sediment distribution, the experi- 
mental values, the exponents were used. Table shows agreement 
between and for only one sediment not possible say that the 
distribution law (Eq. 16) gives the correct values the relative concentration, 
What can said that the distribution function the correct form. 
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Turbulent Transfer turbulent transfer coefficient for 
momentum given Eq.10. Introducing 12b and 15a into Eq. gives 


Eq. 20a assumes logarithmic velocity distribution and linear distribution 
shear and gives the coefficient for momentum transfer. deriving Eq.16, 
the assumption was made that The measured sediment distribution 


now makes possible calculate the coefficient for the sediment applying 
Eq. which gives 


The transfer coefficients may reduced dimensionless form dividing 
the depth times the friction velocity follows: 


and 


Fig. shows both the coefficient for momentum (Eq. 21a) and for sediment 
(Eq. for four flow conditions. Figs. 13(a) and are for series with 
average suspended loads 1.34 and 3.36 per liter (Table 3), respectively, 
and general are typical this series. Figs. 13(c) and 13(d) are typical 
results obtained series the former being for the coarsest sediment (sand 1). 

The agreement between the two mixing coefficients for series (Figs. 13(a) 
and rather good. For the coarse sediment series the sediment 
coefficient smaller than the momentum coefficient, whereas the opposite 
true for the fine sediment. Eq. shows immediately that small values 
give small values the concentration, also apparent from Eq. 
that increasing and therefore will tend make the concentration more 
uniform, that is, will decrease These effects may noted Figs. 13(c) 
and 13(d). Fig. 13(c), less than and the measured exponent 
greater than whereas Fig. 13(d) the opposite less than 


RESULTS 


Effect Suspended Load Velocity Distribution—The von 
universal constant (Eqs. and 18) depends only the mechanism 
turbulent mixing, and change means that the mixing mechanism has 
some way been altered. From and 15a, 
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the momentum transfer directly proportional decrease 


decreases and the shear (depth and slope) kept constant the 
velocity gradient will have increase compensate for the reduced 


Run Number 8 Run Number 14 
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FOR MOMENTUM TRANSFER AND SEDIMENT TRANSFER 


This will result increase the velocity the flow. Thus seen that 
also factor establishing the velocity distribution. 

The values obtained for various sediments and flow conditions are 
shown Table The effect upon increasing the sediment load 
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The results series show that the coarse sediment (0.160 mm), although 
reducing somewhat, has less effect than the finer sediment. general, 
also appears vary systematically with the friction velocity 
However, there appears simple relationship between the several 
factors. From these results one would expect function the 
settling velocity the sediment, the average sediment load, and the friction 
velocity 

Attention called the lack agreement between the values obtained 
series and for similar flows. For instance, the average for clear water 
for series was 0.372 whereas for series was 0.394. straightforward 
explanation these results can offered. However, the difference most 
likely due the difference the inlet conditions for the two series experi- 
ments. The k-value obtained for series which was run with very good 
flow conditions, checks very well with that obtained for pipes, and reason- 
able assume that the results this series are more reliable than those 
series 

The Effect Suspended Load the Resistance Flow.—As stated, 
reduction results increased velocity, which will also give reduction 
the resistance coefficient the channel. The Chézy formula for flow 
conduits is, 


which the mean velocity the cross section, the Chézy coefficient, 
the hydraulic mean depth (area cross section divided wetted per- 
imeter), and the slope pressure gradient. For open channels, 
given the Manning formula, 


which roughness coefficient. The Chézy coefficient can also 
expressed terms the resistance coefficient, 


Eqs. and can solved for and with the following result: 
1.486 


(25a) 
and 


Several values and for flows with and without sediment are shown 
Table 

The effect the suspended load reducing the resistance seen clearly. 
The reduction both and appears about the same for the two sedi- 
ments but greater for the low flows than for the high flows. The first three 
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lines Table indicate that the resistance decreases the sediment load 
increases. 

The increase velocity decrease channel resistance, suspended load 
added the flow, can explained the effect the sediment reducing 


TABLE AND ROUGHNESS COEFFICIENTS FLUME WITH 
AND SUSPENDED FLOW 


anning clear flow coefficient, clear flow diameter, concentra- 
a 


silt flow silt flow (mm) 


Feet per second. Grams per liter. 


the turbulence. keep the sediment from settling, work must done it. 
The energy for this work can provided only the turbulence which 
damped and reduced intensity when gives energy. This means that 
the momentum transfer coefficient also decreased thus allowing the velocity 
and velocity gradient increase according Eq. 1b. (The late Henry 
Eakin, formerly head the Sedimentation Division Research, 
Soil Conservation Service, had anticipated the effect the sediment the 
flow and advanced his idea the subject the writer 1936 during dis- 
cussion the suspended-load problem.) 

Buckley* observed, from measurements the Nile River Egypt, 
that suspended load decreased the roughness coefficient the channel. The 
reduction resistance that obtained for sediment loads from 1.0 
per liter was general greater than that obtained the experiments described 
this paper. Mr. Buckley attributed the reduction resistance the 
damping effect upon the turbulence heavily loaded layer water near 
the bed the stream. expected the greatest reduction when stream was 
depositing, for then the layer near the bottom would most dense. 

Distribution Suspended Load.—All the measurements suspended-load 
distribution followed very closely the form the sediment distribution law, 
Eq. 16, but agreement between the calculated and measured value the 
exponent occurred only for the coarsest suspended load (D, 0.160 mm). 
For the finer sediments the value the exponent that fits the data, the 
measured value, was about 20% less than the value calculated from the theory. 
(Measurements suspended load natural stream made recently Alvin 


Assoc. Am. Soc. E., show that was about half 
Influence Silt the Velocity Flowing Water Open Channels,” Buckley, 
Minutes Proceedings, Inst. E., Vol. 216, pp. 183-211. 


Am. Geophysical Union, Pt. II, 1942, pp. 678-683. 


Run 
No. (ft) 
0.484 3.59 0.0104 1.04 0.0179 1.07 0.160 1.6 
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large asz.) noted previously, the smaller the value the more uniformly 
the sediment distributed over the depth. 

The mixing coefficient given Eq. and has the same form for both 
momentum transfer and sediment transfer but not necessarily the same 
For fine sediment the mixing coefficient for momentum was smaller than the 
corresponding coefficient for sediment transfer. can readily seen that 
increasing will tend produce more uniform sediment distribution, which 
results being smaller than 

The fact that the coefficients sediment transfer and momentum 
transfer are not the same may explained two differences the transfer 
mechanisms for sediment and momentum. The first difference occurs because, 
for momentum transfer, correlation required between the horizontal and 
vertical turbulent velocity fluctuations and uncorrelated random fluctuations 
not contribute the momentum transfer. For sediment transfer this 
kind correlation not necessary and the random fluctuations also transfer 
sediment and tend make larger than The second difference exists 
because the between the fluid and the sediment. The theory assumes 
that the instantaneous vertical velocity the sediment relative the fluid 
always equal its settling velocity still fluid. This true only for steady 
conditions and during acceleration the relative velocity must greater than 
and the sediment carried shorter distance than assumed. This will 
tend make less than 

The experiments show that flows carrying the coarsest sediment the two 
transfer coefficients are about equal and that, finer sediment introduced, 
becomes increasingly larger relative Apparently for the fine sediment 
the effect the random turbulence increasing greater than the effect 
the between the particles and fluid decreasing For the coarse 
sediment the two effects are about balanced and still coarser sediment were 
used resultant effect opposite sign expected and should less than 
The effect the variation the transfer coefficients relative one another 
reflected the relative values and which may noted Table 
For flows carrying the fine sediment where greater than less than 
and when the two transfer coefficients are equal, and are also equal. 

The tendency for the sediment content near the surface the flow 
exceed that indicated the distribution formula, Eq. 16, can explained 
the action the random turbulence. Although the mixing coefficient for 
momentum goes zero the surface, the random velocity fluctuations carry 
sediment right the very surface, thus giving the mixing coefficient for 
sediment finite value the surface. This can noted Fig. 13. 

Attention called the wide disagreement between the values and 
the first two sediment profiles run (Table 3). Considerable difficulty 
was experienced this run because the sediment was not distributed 
regular manner over the cross section. Since the flow was very low and the 
sand coarse, considerable material was moving along the bottom. Any irreg- 
ularity the distribution was therefore easily noticed. The first measure- 
ment (Table and curve Fig. 10) was made over point where considerable 
sediment was motion narrow band the bottom. The abnormally 
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large value could accounted for the occurrence transfer 
sediment the horizontal direction away from the measuring point, brought 
about the existence concentration gradient across the flow. the 
same way possible account for the low value obtained the second 
measurement run 17. This measurement was taken over point where 
very little sediment was motion the bed but few inches one side 
the section considerable sand was being carried. The transverse transfer 
sediment into the measuring section would decrease the vertical concentration 
gradient and therefore the value surprising that even under such 
conditions the distribution should follow the law expressed 

Secondary Flow Due Suspended effect the suspended load 
was first noticed measuring the velocity distribution over cross sections 
the flume. Symmetrical velocity distributions clear flows became distorted 
the addition very small amounts sand. This effect may noted 
Fig. 5(a), distributions and which show the velocity distributions for 
clear water and for average sediment concentrations 0.75 and 3.0 per liter, 
respectively. give idea how small the concentrations actually were, 
might well remark that the first concentration (Cm 0.75) was obtained 
adding only 12.5 sand the entire system. The effect this small 
concentration the flow may also seen comparing the center column 
Fig. 5(a) with the right column. The effect small for this small concen- 
tration, but noticeable. Incidentally, attention called the uniform 
spacing the velocity lines Fig. 5(a), distribution This shows once 
that the velocity distribution logarithmic since the distances from the 
bottom the flume are spaced geometrically. 

The disturbing effect the sediment also shown Fig. 
tions Fig. for clear water, may compared with distributions 
Fig. 5(b), which show the velocities under like conditions except 
that there sediment the flow. Distributions and Fig. 5(b), show 
case which the sediment apparently has stabilizing effect. this one 
instance, the conditions were apparently just right establish good flow 
this particular station. The fact that the distorted velocity distribution with 
clear water was made more uniform the addition small amount 
sediment serves illustrate again that the sediment has considerable influence 
the flow. 

The sediment appeared transported clouds ribbons parallel 
the flow. low flows when appreciable material was moving along the 
bottom, these ribbons appeared bands sand the bottom. When proper 
flow conditions had been attained, these ribbons were symmetrically spaced 
about the center line and were stable regardless the flow rate. Attempts 
break them artificially were only partly successful. Fig. photo- 
graph these bands sediment taken after run 15, looking upstream. The 
photograph was taken after the flow had been stopped and the water slowly 
drained off. Material suspension settled the bottom, thus showing how 
was distributed across the flow. Three main streaks occur the center 
portion the flume and two less pronounced streaks occur the walls. 
This pattern, which was symmetrical, occurred over the entire length the 
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flume. Measurement suspended load distance 0.02 from the 
bottom showed that the sediment was very nearly uniformly distributed over 
the center half the flume, that then suddenly dropped about 25% 
the center concentration the clear areas near the walls (Fig. 14) were ap- 
proached, and finally increased about 50% the center concentration near 
the walls. The concentration for the center in. was about less than 
that over the concentrated streaks the side. (The flume 33} in. wide 
in. deep.) 

The cause these longitudinal bands sediment not known, although 
believed that they are the result secondary circulation and disturbances 
due the presence suspended load. Even the slightest departure from 
uniform sediment distribution will cause horizontal density gradients that can 
set secondary flows. This distribution shows concentration gradient 
across the flow, and one could then expect circulation with upward flow 
the points low concentration and downward flow where the load was 
high. Since there appear four points low concentration, four pairs 
circulations might expected. direct observations were made attempt 
check this hypothesis, and doubtful that this could done without 
considerable work. However, this possible explanation the action 
the sediment causing secondary circulations the flow. Further investi- 
gation will show whether not this idea valid. 

circulation could caused sediment uniform and straight 
laboratory channel, there very good possibility that will also occur 
natural streams where irregularities and disturbances always exist. Therefore, 
some measurements suspended load rivers were examined. Figs. 15(a) 


and show measurements velocity and suspended load cross 


Note the high concentration suspended load about the one-quarter 
points the section. The circulation expected from this distribution shown 
Fig. 15(c). The writer not able state whether this distribution 
typical rivers. However, these measurements the River are evidence 
circulation may caused the sediment. The circulation that 
would induced the same that advanced explanation for the 
depression the maximum velocity natural streams. Although this 
circulation explains the velocity distribution, satisfactory explanation 
was ever found for the circulation itself. The action sediment river 
flow may well prove important factor determining the velocity 
distribution. 

Sediment-Transporting idea that sediment suspension 
reduces the turbulence furnishes possible mechanism for establishing the 
capacity Since the turbulence supports the sediment reduction 
the turbulence adding sediment means that the stream then less 
capable supporting sediment than before load was added. The capacity 
load probably determined continually decreasing the turbulence more 
sediment taken into suspension. The fact that the addition sediment 


Congres International Navigation, Venice, September, 1931, pp. 
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flow increases the velocity, therefore, does not mean that the sediment 
transporting capacity also increased. the contrary the capacity 
actually decreased. The increase velocity also means that the depth and, 
hence, the boundary shear decreased and that particles the bed are prob- 
ably subjected less entrainment force than when the velocity was lower. 

seen not only that turbulence involved the sediment transportation 
problem, but that the sediment changes the turbulence itself, the amount 
change apparently depending upon the amount and size the suspended load. 
clear flow the velocity can expressed terms certain constants which 
characterize the turbulence mechanism. sediment-laden flows these con- 
stants are now functions the sediment load and sediment characteristics. 


CONCLUSIONS 


The distribution the relative concentration suspended load has the 
form the derived distribution law (Eq. 16) but does not agree quantitatively 
with it; that is, the value the exponent given the theory does not agree 
with the exponent that fits the experimental results. the suspended 
sediment fine, less than This shows that the relative concentration 
more uniformly distributed and has greater average value than that 
indicated the theory. other words, flow can support and, therefore, 
transport fine sediment more effectively than the theory indicates. the 
size the sediment increased, the agreement between the two exponents 
improves certain point, beyond which the value may expected 
exceed 

The disagreement between the calculated and measured distributions 
cited conclusion attributed (a) the action the random turbulent 
fluctuations suspending sediment, and (b) the between the fluid and 


.the sediment the sediment accelerated. The theory based the 


assumption that the coefficient sediment transfer equal the coefficient 
momentum transfer, and therefore does not account for the action the 
purely random turbulence and the slip. Such disagreement was anticipated, 
aud the present results confirm it. 

For fine material the coefficient sediment transfer tends exceed the 
coefficient momentum transfer; for coarser sediment the opposite tendency 
found. Conclusion may considered corollary conclusion the 
physical reason for the conclusions being set conclusion 

Suspended load decreases the coefficient momentum transfer, the 
reduction varying with the amount the load and being largest for fine 
sediment. 

Suspended load reduces the value the von universal constant 
turbulent exchange which characterizes the effectiveness the turbulence 
transferring momentum. Reduction means that the mixing less 
effective and would indicate that the sediment tends suppress damp out 
the turbulence. This conclusion follows directly from conclusion 

Suspended load reduces the resistance flow, thus causing sediment- 
laden water flow faster than clear water. This follows directly from con 


e, 
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clusions and the laboratory channel average suspended load 
1.2 per liter reduced the friction factor, much 20%, and the corre- 
sponding Manning roughness coefficient was reduced much 10%. 

Conclusions can all related essentially three effects that 
occur the flow the presence the sediment: (1) The sediment appears 
damp out the turbulence such way that the momentum transfer 
reduced, (2) random turbulence, which not factor the transfer mo- 
mentum, contributes the transfer sediment, and (3) the “slip” between 
the fluid and the sediment tends make the sediment transfer coefficient 
less than the momentum transfer coefficient. 

The idea that suspended sediment damps the turbulence suggests that 
the transporting capacity stream determined through the reduction 
the turbulence the sediment. 

The increase velocity, due the addition sediment clear 
stream, probably does not give rise increased sediment-transporting 
capacity since accompanied reduction turbulence and the 
entrainment force the stream bed. 

10. Suspended load tends cause flow become unevenly distributed. 
Unsymmetrical sediment distribution within the flow may cause secondary 
circulation. 

11. density gradients due suspended load may important 
cause the secondary circulation, which has the effect causing the maximum 
velocity occur below the water surface. 

12. The velocity distribution near the center the flume used these 
experiments was the logarithmic form and followed the von universal 
velocity defect law. 
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PAPERS 


CHARACTERISTIC REDUNDANTS USED FOR 
ANALYZING STATICALLY INDETERMINATE 
STRUCTURES 


Orthogonal functions are commonly used the solution vibration prob- 
lems. lesser degree such functions have been applied some the 
more complicated static analyses, as, for example, plate analysis and the 
determination beam deflections trigonometric series. this paper, 
forces acting relatively simple structures are reduced characteristic forces 
which are orthogonal, and methods analysis based these characteristic 
forces are discussed. 

The heart the multiply indeterminate solution may interpreted 
means the determination so-called characteristic forces. After these 
forces have been found, redundant stresses corresponding any condition 
loading distortion may evaluated without the further solution simul- 
taneous equations, regardless the degree which the structure may 
statically indeterminate. 

not claimed that the methods suggested are either more quickly applied 
more easily understood than other methods which are commonly used. 
hoped, however, that they will interest because they differ radically 
from the usual methods. Being different, possible that other applica- 
tions the principles involved may lead methods solution which con- 
stitute improvement the methods outlined herein. 


The symbols used this paper are defined where they first appear, either 
text diagram, and are assembled for convenience reference the 
Appendix. 


comments are invited for immediate publication; insure publication the last 
discussion should submitted November 1944. 


Prof., Structural Eng., Mass. Inst. Tech., Cambridge, Mass. 
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CHARACTERISTIC FoRCES 


Suppose elastic structure acted upon number concentrated 
forces Ps, ---, Pn, which are applied points ---, respectively, 


Py, 


Let values ---, determined such manner that 


(1) 


having the same value each Eqs. 
set forces which satisfy Eqs. will called 


DETERMINATION CHARACTERISTIC FORCES 


value P,, then, superposition: 


Since the right-hand side Eqs. zero, the numerator determinant will 
zero. Hence, order for the forces ---, have values other than 
zero, necessary that the denominator determinant formed the coeffi- 
cients Ps, the left side these equations also equal 
zero. This condition represented 


€ 
i 
n P,, 
| 
A, 7 
TANS 
be. . . . . . . . . 


~ ' 
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The expansion the determinant Eq. leads nth degree algebraic 
equation terms forces are being considered, there will values 
which will satisfy this equation. These values will denoted 
Kn, taken order descending magnitude. 

Each these values may now substituted separately into Eqs. 
The simultaneous solution equations Eqs. will then permit 
one these forces such Pam. nth equation Eqs. may used 
check both and the solution the equations. Since there are 
values km, sets characteristic forces exist, which are determined the 
foregoing procedure. The complete set characteristic forces corresponding 

There will result then, the following sets characteristic forces: 


For 

for 

and, for 


Poa Pas; Pin Pos: Pus @ Aes Pen 


Amn having been determined the simultaneous solution equations 
Eqs. 


Section methods were presented whereby each set characteristic 
forces could expressed terms any one force the set, such Pam. 
The relations Eqs. not define Pam, however, which may given any 
value without altering the characteristic aspect the loads. matter 
mathematical convenience the development follow, Pam chosen arbi- 
trarily that: 


satisfy Eq. substitute values from Eqs. and solve for Pam: 


1 + (Asm)? + + (Anm)? 


When values are determined this manner, set characteristic 
forces said “normalized.” the following treatment, normalized 
characteristic forces will denoted Thus represents the normalized 
characteristic force corresponding kn. Similarly, represents the 
complete set normalized characteristic forces for points corre- 
sponding km. Using this new notation, Eqs. and respectively, 


| 
é 
’ 
’ 
) | 
r 
. 
| 
: 
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should rewritten 


——— ..(8 
For 
for 
for Kn, 
and (see Eq. 6): 


ORTHOGONALITY NORMALIZED CHARACTERISTIC FORCES 


Suppose that structure acted upon one set normalized characteristic 
forces subjected virtual distortion corresponding the deformation 
that occurs when second set normalized characteristic forces applied 
the structure. The external virtual work done the forces Q,, may 
expressed 


Now suppose that the same structure, when acted upon the forces 
subjected virtual distortion corresponding the deformation that occurs 
when the forces Q,, are applied the structure. The external virtual work 
done forces may expressed 


Betti’s law? states that the virtual work done group external forces 
during distortion due group external forces equal the 
work done the forces during distortion due the forces 
Hence, and can equated from Eqs. and 12, respectively, lead- 
ing 


Since the parenthetical expressions the two sides Eq. are identical, 
and, since Km, follows that the equation can satisfied only when 


Qam Qan + Qom Qin + + 0 (14) 


Thus, for any two sets normalized characteristic forces, the sum the cross 
products forces applied the same point vanishes. Thus the two sets 


Co., New York, Y., 1937, 39. 


i 
: 
| 
| 
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normalized characteristic forces are orthogonal; and orthogonality does not 
depend normalizing. 


First Work THEOREM 


Because the orthogonality normalized characteristic forces, follows 


Thus the external virtual work done any one set normalized charac- 
teristic forces Qm, during distortion caused any other set normalized 
characteristic forces equal zero. This conclusion will called the 
first virtual work theorem for normalized characteristic forces. Its validity 
does not depend normalizing. 


Assuming that follows from Eq. that 


since the forces are normalized. Thus the virtual work done given 
set normalized characteristic forces during distortion caused the 
same set normalized characteristic forces Qn, equal km. This conclu- 
sion will called second virtual work theorem for normalized charac- 
teristic Its validity does depend normalizing. 


CHARACTERISTIC 


Referring Fig. let ---, any set actual forces applied 
points ---, respectively, along the lines action the characteristic 
forces Ps, respectively. These actual forces may expressed 
terms normalized characteristic forces superposition follows: 


which B,, may called the participation factors for the first, 
second, ---, and nth sets normalized characteristic forces, respectively, with 
respect forces Ro, ---, are known, the participa- 
tion factors may evaluated follows: 
With the forces ---, acting the structure, suppose the struc- 
ture undergo virtual distortion due the further application the mth 
set normalized characteristic forces Qm. The virtual movements points 
Anm The virtual work done forces Ra, Rs, given 


| 
q 
7 
q 
? 
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Substituting the relations Eqs. into Eq. simplifying: 


+ B, (Qa2 Gus + Qom + + Qne Gun) 


+ Ba (Qan Qam + Qon Qom + Gin cee 
From the relations Eqs. and follows that 


Equating Eqs. and for and solving for 


With values determined 20, the superposition Eqs. com- 
pletely defined. 


RECIPROCAL RELATIONS 
Letting B,, defined Eq. substituted into the mth equation 
17, and simplifying, 


Rin Ra (Qa Qmi + Qaz Q m2 + + Qmn) 
+ Ry (Qo: Qmi Qme + + Qon Qua) 


+ R, (Qui Qmi + Qno Qm2 + sat + & ale (21) 


signing value R,, but making all other values equal zero: 


assigned value, but all other values taken zero: 


Qni Qmi + Qne + + Qan Qmn = 0 


sense, Eqs. and are reciprocal expressions Eqs. and 14. 


DEFLECTIONS 


Referring Section let required compute the total deflection 
these actual forces may replaced the characteristic forces 


10. CHARACTERISTIC REDUNDANTS 


Let required find the redundant stresses ---, inde- 
terminate structure such that shown Fig. 2(a), due any condition 


Papers 
(19a) 
: 
q 
ar, 
+ 
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external loading imposed strains such those caused temperature 
changes pier settlements. These redundant stresses may considered 
external loads acting statically determinate primary structure, shown 
which represent the relative movements points respectively, due 
all causes, are each instance equal zero. 


The redundants Xs, ---, Xn, however, may expressed terms 
normalized characteristic forces, accordance with follows: 


+ Br Qon 


Thus one may replace the redundant forces shown Fig. with the char- 
acteristic forces shown Fig. 2(c), wherein, the normalized characteristic 
forces ---, have previously been determined, the participation fac- 
tors B,, are the unknowns that must determined order that 
Xa, may computed from Eqs. 25. 


cannot computed Eq. 20. One may proceed follows, however: 

evaluate apply the mth set normalized characteristic forces 


| 
(c) 
ne 
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forces that act the primary structure during its distortion due all causes 
combined. Note that consists pair equal and opposite forces 
applied points etc. 

The distortion the primary structure may considered due three 
classes causes: 


(a) Distortion due external forces, temperature changes, pier settle- 
ments, etc., including all causes except the redundant stresses themselves; 

Distortion due the redundant B,, (which consists forces Qam, 
Bn Qom, Bn Quam); and 

(c) Distortion due all other redundants—that is, Qn, 
but not including Qn. 


Let Sno the relative deflections points ---, respec- 
tively, due the first class causes, these deflections being taken positive 
when they act the directions ---, Xn, respectively, shown 
Fig. Then, assuming the virtual forces act the directions shown 
Fig. 2(c), the external virtual work done the virtual forces during 
this part the distortion given 


Due redundant the external virtual work done the virtual forces 
seen the second virtual work theorem (Section given 


Due each the other redundants, the external virtual work done the 
virtual forces seen the first virtual work theorem (Section 
given 


The total external virtual work done the virtual forces given 
This total must equal zero since points ---, undergo 
relative movements when all causes are considered. Hence 


For structures direct stress, the stress any bar the primary 
structure due unit value Q,m, and the change length any 
bar the primary structure due external loads temperature changes, 


= DS [AL (Qam + Qom Fom + + Fum) ] (28a) 
but Fam Fom Fam equals the stress any bar the 


primary structure due all the virtual forces Q,, acting simultaneously, and 
may denoted Fm, whence 


| 
d 
q 
7 
= 
: 
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Considering only the action external forces, given 


which the stress any bar the primary structure due external 
forces only. Therefore, 


and Eq. becomes 


may noted that, although the determination characteristic loads 
involves the solution Eqs. the evaluation redundant stresses corre- 
sponding any condition loading distortion involves further solution 
simultaneous equations, regardless the degree which the structure 
statically indeterminate. 


12. SOLUTION 


The use characteristic redundants analyzing statically indeterminate 
structure illustrated example evaluating the stresses and 


TENCE. 


the structure shown Fig. For convenience, assume that for each 


member, 25. 


Determination Normalized Characteristic steps are required 
normalize the characteristic loads: (a) Evaluation (b) evaluation 
Km; (c) determination characteristic loads; and (d) normalizing. 


(0) TX S ; 
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(a) Evaluation (By Method Virtual Fam, Fom, and 
the stresses any bar due the condition loading shown Figs. 
3(c), and 3(d), respectively. The computations for are arranged shown 


TABLE 


Fun L 


Fon Foe 


(b) Evaluation (Using Eq. becomes— 


which reduces 


(c) Determination Characteristic Loads (Using Eqs. 3).— 
For 


Bar 


o 


32) 


3a) 


June, 1944 839 


For 


For K3, 


Hence the characteristic loads are defined by: 


Normalizing Characteristic Loads.—Substituting Eq. 


Eq. the normalized characteristic loads are given by: 


Evaluation first step determine the value 
Let and the stress any bar due the condition loading shown 
Figs. 3(e) and 3(f), respectively. Note that the F,,-stresses may obtained 
from the Fam, and stresses given Table superposition. The 
arrangement the computation shown Table which supplies the re- 
maining values substituted Eq. 31, thus: 


™ 


n 
(35) 
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0.4890 and 0.1601 The redundants are then determined 
substitution Eqs. 25, thus: 


Hence, 


TABLE 


will noted that the additional computation required investigate 
second condition loading not very great, since modifying the 
will change only part the computation for and consequently only 
part the final step. 


13. ALTERNATE ANALYZING STATICALLY INDETERMINATE 


the number independent components displacement for structure 
less than the degree which the structure statically indeterminate, 
may desirable modify the details the analysis. 

Characteristic forces which correspond the independent components 
distortion may determined and normalized. For the actual external Icad- 
ing condition investigated, participation factors for the normalized char- 


forces may evaluated. Deflections the points application 


these forces may then computed, following which, the evaluation stresses 
the structure straightforward. 

Suppose, for example, required that the end moments the frame 
Fig. 4(a) determined, due the load acting shown. Corresponding 


| 
Bar 12 Fo 10 Fi 
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the three independent components distortion and there will 
considered set characteristic forces the external moments and 
the external force respectively, acting shown Fig. 


statics, the following three equations may written: 


Expressing end moments the slope deflection equations, these reduce 


—d=0 


and X», whence 


Eqs. correspond Eqs. and Eq. becomes 
—4k 


Successively substituting these values into Eqs. 39, the characteristic 
forces aré determined follows: 


Using Eq. normalize the characteristic forces, 


= Xai} Xu = Xa = 19.77 Xai 


| 
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istic forces are, therefore, 


The participation factors applied each these sets normalized 
characteristic forces order that they may correspond the loading shown 
this condition loading are now computed Eq. 24: 


Therefore, and the end moments the frame are 


follows: 


14. INTERDEPENDENCY CHARACTERISTIC 


For the general case, Eqs. consist simultaneous equations with 
unknowns. The solution then requires the evaluation the nth order deter- 
minant Eq. and, for each the roots thus obtained, the simultaneous 

Under certain circumstances, however, the points application and the 
directions characteristic forces may such that the equations Eqs. 
may break down into two more groups, each which may solved 
itself. For example, one might obtain for Eqs. for five forces Pe: 


(45a) 


This condition would arise, for example, neither the forces caused 
displacements points along the lines action any the forces 


7 
‘ 
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P., Pa, respectively. Eqs. there are two separate groups. Eqs. 
45a and 45b may solved together, involving only second-order deter- 
minant. From these two equations two sets characteristic forces are found, 
for each which and have val- 

Eqs. 45c, 45d, and 45e may then 
solved together, involving third-order 
determinant. From these three equa- 


tions three sets characteristic forces 

are found, for each which P., Pa, and 

have values, whereas and are 

zero. 
When break down this 


manner, considerable saving compu- 
tation possible, illustrated with 
reference the structure Fig. 
taking the redundants shown 
Fig. Values may deter- 
mined the usual methods structural analysis, leading 


Eqs. become: 


First, solving Eq. directly, 116.6 and Next, solving 
47a and 47c simultaneously, 


+24 


’ 
N 
o 
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Thus the normalized characteristic forces are: 


and, for 6.00, 


the usual methods structural analysis, the following deflections due 


Eq. 27: 


B, = 


and, Eqs. 25: 


The end moments the members the frame may now computed 
statics, and will found equal those obtained Section 13. 


the analysis fixed end arch, well known that the redundants 
may advantageously taken acting the elastic center the arch. The 
one wherein the redundants are 
characteristic forces, chosen such 
manner that Eqs. break down into 
three independent equations. illus- 
trate the foregoing, consider the struc- 
ture Fig. 4(a) taking the redun- 
dants acting the elastic center 
shown Fig. and applying the pro- 
cedures developed this paper. For 
this structure, the elastic center lies 
the vertical axis symmetry for the 
and 4.00 below the girder 

the structural analysis, leading to: 


El=5 


Center Line 


12! 


Seca = 0 ; => 0 Gos = 6 00 


7 
q 


Tm 


| 
| 
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become: 


Solving these equations separately: 
Ka = + 96.0; Pas = + Ps ee (54) 


Application Eq. leads the following sets normalized characteristic 


usual methods structural analysis, the following deflections due load 
27: 


and, Eqs. 25, 
(56) 


The end moments the members the frame may now computed 
statics, and will found equal those obtained Section 13. 


16. ReMARKS CONCERNING SOLUTION Eas. 


Section was assumed that all values were real and separate. 
all solutions made the writer, values have been real. Should imagi- 
nary values encountered, believed that the physical conditions 
stated Eqs. are impossible fulfilment, and that result the analysis 
cannot properly applied. considered doubtful that this condition will 
arise. 

possible for two more values equal numerically, however, 
although this condition infrequent. Suppose that the values 
given value appears times. will then found that only equa- 
tions Eqs. are independent when substituted, that only values 
can expressed terms the remaining values This leaves 
degrees freedom for each value x’, that m(m further 
conditions must imposed. These conditions must lead separate sets 
characteristic forces, which must mutually orthogonal. Within the 
limitation these requirements, arbitrary conditions may imposed, since 


q 
)- 
7 
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Eqs. will satisfied. The foregoing may illustrated the following 
example: Consider the equations 


The determinant formed the coefficients forces reduces to: 16.00) 
substituted into Eqs. 57, only two equations are formed, which 
reduce 


Since 58a not define the relations between and any assumption 
may made regarding this relation, and Eqs. will satisfied. Suppose 
that for one set characteristic forces, one such assumption made arbi- 
assumption defines the following set characteristic forces: 


For the second set characteristic forces, second assumption such 
may made, which another numerical factor which cannot, 
however, chosen arbitrarily. This defines the second set characteristic 
forces as: 


The two sets characteristic forces defined Eqs. and must 
satisfy the orthogonal requirements Eq. 14, whence 


and Thus for for example, one may choose arbi- 
trarily; that is, one may select thus defining the first set char- 


acteristic forces as: 


Then for xs, must computed from the relation 1/C 
thus leading the second set characteristic forces: 


Since these two sets characteristic forces are both associated with the same 
numerical value the distinction which shall given the subscript 
and which the subscript arbitrary. different assumption the value 
would have led different sets characteristic forces, but believed 
that the over-all solution would not modified thereby. 


4 
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17. ALTERNATE FOR ORTHOGONALITY 


the treatment this subject, sets characteristic forces have been 
used, which were orthogonal with respect each other. not necessary 
that the orthogonality exist between sets forces; alternate basis for pro- 
cedure would operate sets forces chosen that the gth set forces 
orthogonal with respect the deflections caused the rth set forces. 
This condition may illustrated the equation 


61, the force applied point the set such forces, denoted 
Amr denotes the total deflection point due the rth set such 
forces. 

such sets forces, each composed forces, are determined, 
there are all forces which values must assigned. For the first, 
forces, respectively, the only limitation assigning these values being that 
each set forces shall independent with respect the other sets. The 
remaining forces, which consist forces for the first, second, 
sets forces, respectively, must computed the basis Eq. 61; 
this will necessitate the respective solution 0,1, sets simulta- 


neous equations. the judicious choice the arbitrary values 
assigned (such assigning values and 0), the amount com- 
putation arriving the dependent values may greatly reduced. 


truss statically indeterminate the nth degree, let the redundant 
follows: 


= Ci Sa + C2 + + Cn San 


which C2, are participation factors for the Ist, 2d, ---, nth sets 
forces, respectively. Let the set forces Sag, Sng denoted 
let the total stress any bar due acting the primary structure 
denoted F,. Let the total stress any bar due external loads 
acting the primary structure. 

Suppose the primary structure, acted upon the external loads and 
all the redundants ---, Sa, subjected virtual distortion 
due increase the value ACm. Then the redundant 
will external virtual work which may expressed 


The other redundants will external virtual work accordance with 
Eq. 61. The external loads will external virtual work, which may ex- 
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pressed 


Thus the total external virtual work, which equals zero, given 


whence 
Am) 


Hence sets forces satisfying Eq. can obtained, Eqs. and 
lead the evaluation Xs, ---, The manner which these sets 
forces may obtained, well the application this entire procedure, 
may illustrated consideration the structure Fig. 3(a), for which 
member. the usual methods structural analysis (see 


Section 12) the following values are obtained: 


Referring Fig. the first set forces 


100 that Eq. obviously satisfied. For the third set 
satisfy Eq. 61. Considering the first and third sets forces: 
0.160. For the forces: 


summarize, the three sets forces are: 
For 
for 
Sa2 = 0 ; Spo = 0 > Se = (68) 
and, for Ss, 
Sas =o 0.160; Sos = i Ses 0.160 | 


4 7 
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The total deflections due these forces are, respectively: 


For Si, 

for So, 

Ago = 0; Ase = 16 ; Aw = (69) 
and, for Ss, 


Values are obtained from Table wherein values may 


TABLE 3.—EVALUATION ALTERNATE 


LeLs —0.600 +0.096 
—0.833 +1.000 —0.160 —0.833 +0.133 
>= 


obtained statics. Values F,, are also obtainable statics; they may 

obtained superposition from stresses due unit values and 

which will already available, since they will have been used (see Section 12) 


evaluating For substitution into Eq. 64: 1.667 (25) 


Hence, Eq. 64: 


and 


q 
4 
| 
| 
| 
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Finally, Eqs. 62: 


further condition could have been imposed with respect each set 


forces: 


This would amount normalizing each set forces with respect the deflec- 
tions caused that set forces. The detail the solution would then 
change slightly; replaces indicate that Eq. has been observed, 
and interpreted the stress any bar due acting the primary 
structure, Eq. 63a becomes and Eq. 63b remains unchanged. 
The denominator then drops out Eq. 64. Eqs. 62, should 

This type solution possible for any degree indeterminacy struc- 
ture, although admittedly the labor involved determining sets forces 
which satisfy Eq. would increase rather rapidly with the number un- 
knowns. Once such sets forces have been established, however, the redun- 
dant stresses may determined without further recourse simultaneous 
equations. This the same advantage that gained use the elastic 
center, which, however, limited its direct benefits structures with three 
fewer redundants. 


18. 


The purpose this paper has been present some the principles in- 
volved making use orthogonal forces the analysis statically indeter- 
minate structures. examples demonstrate these principles, but 
the structures analyzed were chosen because their suitability for illustrating 
principles rather than because their solution the methods discussed was 
advantageous. 

Structural engineers have been using orthogonal forces for some time 
various analyses, but, some instances, perhaps, without realizing that they 
were doing so. hoped that the presentation this paper will stimulate 
further consideration the possible advantages that may obtained 
their use. 


APPENDIX 


The following letter symbols, use the paper, conform essen- 
tially American Standard Letter Symbols for Mechanics, Structural Engi- 
neering and Testing and American Standard Letter Symbols for 


ASA—Z10a—1932. 
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Mechanics,‘ prepared committees the American Standards Association, 
with Society representation, and approved the Association 1932 and 


1942, 


> 


respectively: 


also (see Eq. 30), cross-section area structural member; 

participation factor Eq. 17; ---, are participation factors 
for the first, second, and nth sets normalized characteristic forces, 
Ra, Ro, Ra, respectively; 

numerical factor Eq. also, Ci, C2, ---, are participation fac- 
tors (Eq. 62) for the sets forces S,, S2, ---, respectively; 

numerical factor Eq. 58c; 

lateral distortion the top column (Fig. 

modulus elasticity; 

total internal stress member; thus, the stress any bar 
the primary structure due unit value the stress 
any bar the primary structure due external forces only; 

moment inertia; 

length structural member; 

number times that the value appears equation; also used 
general integer; 


external load (Eq. 2); 


normalized characteristic force; thus, Qam (Eq. represents the com- 
plete set normalized characteristic forces for point corresponding 
Km; 

particular set forces (Eq. 61); 

resultant external forces (Eq. 18); 

particular set forces (Eq. 61); 

set forces; Smg the force applied point the gth set such 
forces; 


set forces replace Eq. 72; load applied the primary 


structure which produces stress 

external virtual work; the external virtual work done the 
forces Q,, when structure acted upon one set forces sub- 
jected virtual distortion corresponding the deformation that 
occurs when second set forces applied the structure; 


redundant total stress; 
total deflection; the total deflection point the direction 


due forces Ps, Pa; 

deflection due unit load; the deflection any point the 

angular distortion, the subscript denoting point rotation; 


proportionality factor; given value that occurs times for 
rotation factor Eq. 44. 
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PAPERS 


TRANSLATORY WAVES NATURAL 
CHANNELS 


The phenomena translatory wave travel natural channels have been 
studied for many years without entirely satisfactory results. Published studies 
have been largely theoretical, supplemented laboratory tests some in- 
stances, but with few observed data waves natural channels. Stream 
flow regulation storage reservoirs the Tennessee River Basin during the 
past few years has developed numerous translatory waves various magni- 
tudes. Under such conditions, waves can selected for study which travel 
considerable distances without being appreciably affected tributary inflows. 

This paper presents data natural channel translatory waves the 
Clinch and the lower Tennessee rivers, and Wheeler Reservoir Alabama 
for considerable range flow each type channel. Theoretical velocities 
are computed for these waves formulas and are compared with observed 
velocities. Results indicate that wave velocities the Clinch River can 
computed with reasonable accuracy, and that velocities for the rising face 
wave the lower Tennessee River can approximated. Computed veloci- 
ties for waves Wheeler Reservoir are poor agreement with observed 
velocities. 


TRANSLATORY WAVES THE CLINCH RIVER 


Regulated discharges Norris Dam Tennessee have developed num- 
ber well-defined translatory waves which were not affected appreciably 
tributary inflow the reach the Clinch River between Norris Dam and 
Wheat stream gage 65.3 miles downstream (see Fig. 1). The Clinch River 
Norris Dam drains mountainous area and has average annual flow 
4,600 per sec. Below the dam the river drains hilly country and has 


comments are invited for immediate publication; insure publication the last dis- 
cussion should submitted November 1944. 
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average slope about 1.6 mile the reach from the dam the 

Wheat gage. 
Continuous stage recorders located intervals throughout this reach pro- 

vide records from which gage-height graphs Fig. are plotted. Norris (2) 


TABLE CHANNEL AREAS, 


height 


Gage heights Norris, see Fig. 


area (miles per area (miles per area 
(sq ft) hr) (sq ft) hr) (ou per 
(2) (3) (4) (5) (6) (7) (8) (9) 
420 100 1,260 0.9 120 
790 1,150 2.1 1,120 1,750 1.8 2,030 
2,210 1,550 2.8 1,850 2,010 2.1 3,680 2,400 3.7 
4,050 1,950 3.5 2,710 2,270 2.4 5,980 2,770 4.5 
6,310 2,360 4.0 3,690 2,540 2.6 8,500 3,150 4.5 
9,000 2,790 4.4 4,790 2,830 2.8 11,050 3,550 
11,860 3,220 6,010 3,130 3.0 13,600 3,950 
14,760 3,650 7,350 3,440 3.2 16,150 4,360 
17,660 4,100 8,810 3,750 3.4 18,700 4,780 
20,560 4,550 10,390 4,050 3.6 21,250 5,220 3.9 
23,500 5,000 12,090 4,360 3.8 23,800 
26,400 5,450 13,910 4,670 4.0 
29,500 5,900 4.5 15,850 4,990 4.2 
32,800 6,350 17,910 5,310 
35,900 6,800 20,090 5,640 4.6 
22,390 5,980 4.7 
27,350 6,680 4.8 


and Wheat gages have been rated the Geological Survey, and ratings 
for Clinton and Edgemoor gages are developed from these ratings. Discharge 
ratings and corresponding channel cross-section areas for these stations are 


Breedenton 
Gage 


€ 
Scale Miles 


givenin The gage-height relation between Wheat gage and Gallahers 
Ferry gage also given. discharge measurements for rating Wheat gage 
are made Gallahers Ferry, data the channel cross section Wheat 
gage are not available. Pertinent data these gages are given Table 


@ 
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Velocity Translatory Wave.—To determine the velocity translatory 
wave from observed data, necessary identify identical points the 
wave successive stations. this study (1) the beginning-of-rise and 
(2) the midpoint-of-rise are used for this purpose rising waves, and (3) the 


DGEM0g 


5,220 


5,670 


Norris 


Wheat gage Wheat gage 
(cu per Relative (miles per Wave Water 
height ft) (miles per (miles per 
(11) (12) (13) (14) (15) (16) 
: 2.8 0 2.05 1,030 aa 2.0 0.4 2 
590 3.52 1,460 1.5 2.7 0.8 
1,780 4.93 1,890 2.2 3.2 1.2 
3,300 6.27 2,310 2.6 3.6 1.5 
864.3 4,900 7.58 2,730 2.9 3.8 1.8 6 
ps 64,1 6,610 8.75 3,110 3.1 3.9 2.1 7 
p 64.0 8,360 9.88 3,500 3.2 4.0 2.3 8 
10,160 11.01 3,880 3.4 4.1 2.5 
3.9 12,000 12.07 4,250 3.5 4.15 2.7 10 
63.8 13,940 13.16 4,630 a 4.2 2.8 1l 
15,940 14.21 5,010 3.6 
17,940 15.25 5,380 
20,040 16.28 5,770 3.7 14 
q 22,140 17.31 6,150 15 
4 24,340 18.30 6,530 3.8 16 


beginning-of-fall and (4) the midpoint-of-fall falling waves. Velocities indi- 
cated the points fully developed rise and fully developed fall were inves- 
tigated but were found unsatisfactory, partly because these points could 
not identified accurately and partly because these points are modified 
considerably 


ences. TABLE River 
The following formulas? are 
Drainage 
used for computing theoretical Station River mile area 
velocities from channel charac- miles) 
teristics: Norris Dam........ 79.8 2,912 
stics Norris 78.7 2'913 
Clinton gage....... 59.0 778.54 3,056 
0.682 Edgemoor 47.6 768.08 3,089 
0.682 (1b) Elevation gage zero above Mean Gulf Level. 


the average velocity the wave, taking into account the channel discharge 


Hydraulics Flood Movements Rivers,” Harold Thomas, Engineering Bulletin, 
Carnegie Inst. Tech., 1934, pp. 


ris 
urea A 
(sq ft) 
(9) 
1,350 
1,690 
2,040 
2,400 
2,770 
3,150 
3,550 
3,950 
4,360 
4,780 
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characteristics throughout the entire rise the wave; and are the dis- 
charges corresponding gage heights and respectively; incre- 
the channel; and and respectively, are the cross-section areas below 
stages and 

applicable only for determining the velocity wave traveling 
stream given stage. This becomes evident when the derivation 
the equation examined. applicable for determining the average 
velocity the wave, takes into account the channel discharge charac- 
teristics throughout the entire rise the wave. Results from both equations 
will the same when the discharge characteristic the channel, that is, 


constant throughout the range stage covered the wave. will 
noted that function the slope the rating curve given stage, 
whereas function the slope secant the rating curve between two 
stages, the average slope the rating curve. The constant 0.682 Eqs. 
inserted convert miles per hour. 

Using Eqs. wave velocities are computed for the four identifiable points 
each wave each the four stage recording stations. These computed 
velocities are assumed exist half the distance adjacent stations deter- 


TABLE (Eq. la) For 


Observed 


Station (ft) (ft) (Eq. la) (miles) (hr) 
(1) (2) (3) (4) (5) (6) (7) (8) 
Norris 4.0 417 210 0.1 3.43 9.85 2.87 
Clinton 5.5 282 110 0.1 2.66 15.55 5.85 
Edgemoor 3.2 369 210 0.1 3.88 22.25 5.72 
Wheat 5.8 139 0.1 2.92 16.55 5.67 


1.—Stations which data gage heights are available. 

Col. 2.—The river stage the beginning-of-rise, taken from the gage-height graph for each station 
Fig. data the channel Wheat are not available, the stage Gallahers corresponding 
the observed stage determined from the gage-height relation Table for computing the velocity 
this station. 

Col. 3.—Channel widths these stages are determined from the rate change channel cross-section 

Col. 4.—The rates change discharges, dQ, per 0.1 change stage, these stages are determined 
from the discharge rating curves Table 

Col. 5.—0.1 change stage. 

Col. 7.—The reach length miles the distance over which the computed velocity each station 

Col. 8.—The number hours required for the wave traverse this distance traveling 
The average velocity ci, from Norris Wheat, then equal 64.2 miles divided the total number 

hours required for the wave traverse these four reaches. 


mining the theoretical time required for the wave traverse each the three 
reaches between stage recording stations. The average velocity between 
Norris gage and Wheat gage then computed from the total time required 
for the wave traverse the three intermediate reaches. 


0} 


fro 

Gal 

Col 

Col 

Col 

Col 

Col 

Col 

Col 

Col 

ti 

0 

| 


June, 1944 TRANSLATORY WAVES 859 


The method applying Eqs. illustrated Tables and using data 
from Wave October 1937. 

curves for values can developed for each station shown Fig. 
curve for the average velocity for the Norris-to-Wheat reach can also 


(1) (2) (3) (4) (6) (7) (8) (9) (10) (11) (12) 


Norris 4.0 4,050 3,560 1,950 1,610 4.30 9.85 2.28 
Clinton 5.5 4,210 4,830 2,680 2,150 3.38 15.55 4.61 
Edgemoor 3.2 14,620 4,080 4,110 2,470 1,640 4.38 22.25 5.08 
Gallahers 7.3 14,540 4,590 4,736 2,642 2,094 


Col. 1.—Stations which data gage heights are available. 

Col. 2.—Gage height end-of-rise from Fig. 2(a) (observed). 

Col. 3.—Gage height beginning-of-rise (observed). 

Col. 4.—Discharge corresponding gage height determined from Table 

Col. 6.—Channel cross-section area A2, corresponding gage height determined from Table 

Col. corresponding gage height 

Col. from Cols. and 

Col. from Cols. and 

Col. 10.—Wave velocity computed substituting Cols. and 1b. 

Cols. and and (Cols. 10, 11, and 12) are determined the same manner corresponding 


developed determining stages successive stations corresponding 
selected stage Norris gage mean relations and computing 
the average velocity Such curve indicates average values plotted 
against Norris stages for the Norris-to-Wheat reach under average conditions 
local inflow. gives the same value computed Table This 
occurs only when average conditions local inflow exist and the stages all 
stations agree with the mean gage-height relation used for constructing this 
curve. 

For the purpose comparison, curve showing the water velocity the 
Norris-to-Wheat reach also shown Fig. (water velocity the velocity 
indicated current meter). This curve developed from the discharge 
rating and channel cross-section areas each station for stages determined 
gage-height relations, assuming that each station velocity exists halfway 
adjacent stations. The water velocity about half the wave velocity 
ordinary stages and less than half low stages. 

charges and successive stations will noted. This can caused 
either error the station ratings, errors observed stages, irregularities 
flow. believed that this inconsistency will not affect results materially 
long the area difference, corresponds the discharge difference, 

Theoretical velocities and for the beginning and midpoint rising 

waves, and for the beginning and.midpoint falling waves shown Fig. 
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together with the velocity corresponding the observed time travel 
from Norris Wheat, are shown Table 

Comparison Observed and Computed Wave shows veloci- 
ties Co, and plotted against the gage heights observed Norris gage, with 
mean curves drawn for each. Mean curves computed velocities for the 


Gage 


Whea 


Gage Height Norris 


Feet 


Norris Gage 


Height, 


Gage 


Velocity, Miles per Hour 


Norris Dam WHEAT GAGE 


midpoint-of-rise (Fig. 4(b)) and the midpoint-of-fall (Fig. 4(d)) agree very 
well with the mean curves for observed velocities. will noted that the 
range gage heights covered these data small, and, the gage heights 
being midpoints, not extend either extremely low high stages. Results 
for the beginning-of-rise (Fig. 4(a)) and the beginning-of-fall (Fig. not 
agree nearly well, computed velocities being smaller than the observed 
both cases. the latter plottings, velocities are computed from the total 
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rise fall the wave and are plotted against respectively. Mean 
curves for observed velocities for both the beginning-of-rise and the beginning- 
of-fall are higher given stage than observed velocities for the midpoints. 

Wave February wave February 1938 (Wave Table 5), 
developed when the quantity released Norris Dam was increased during 
comparatively low and steady flow the Clinch and upper Tennessee rivers. 


TABLE AND WAVE VELOCITIES BETWEEN Norris 
AND WHEAT River PER Hour 


FALLING WAVES 
1937 
1 Oct. 7.8] 3.8 | 3.2 5.9 5.6)4.0] 4.6 | 3.8 4.9 |3.1) 3.6/3.5 
Oct. 9 |6.8] 7.9] 4.3 | 4.0 7.4 |3.9|4.0/4.0] 4.3 | 4.0 §.2 |3.513.713.5 
1938 
10 Dec. 14 | 5.4] 6.0] 3.7 | 3.8 5.7 6.0)5.5] 4.4 | 3.8 §.8 13.6] 3.8] 3.7 
15 |5.5] 6.0) 3.6 | 3.8 5.7 |3.213.8]3.7] 6.0)5.4] 4.1 | 3.8 5.7 | 3.7] 3.8/3.8 
12 Feb. 8 [1.2] 6.7| 2.4 | 1.2 4.0 |2.8/2.8/3.0| 6.8}1.2| 4.9 | 3.9 4.0 |3.0] 3.2/3.0 


the wave velocity corresponding the observed time travel. 


The release was gradually increased from per sec 24,300 per 
sec hr, was held about steady for hr, and then was gradually reduced 
about the original release hr. Gage-height graphs the resulting 
wave Clinch River stations and Chickamauga Dam site, near Chatta- 
nooga, Tenn., are shown Fig. 2(c), together with the 
water samples from the Chattanooga City Pumping This intake 
about six miles below Chickamauga Dam site and the water travels this 
distance about hr. Prior the arrival the wave from Norris Dam, 
soap-hardness tests Chattanooga showed fairly constant value 60. 
However about after the arrival the wave the hardness increased 
value for 2-day period. has been found that the hardness Norris 
Reservoir water fairly constant ppm, whereas the hardness Tennessee 
River water Chattanooga averages about ppm, Norris water being con- 
siderably harder than water from most the tributaries above Chattanooga. 
There marked similarity between the pattern the points for soap hard- 
ness and the shape the Chickamauga Dam gage-height graph during the 
arrival the wave from Norris Dam, the peak hardness occurring about 
two days later than the peak flow. unfortunate that data are not avail- 
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able the hardness the water the several tributaries coincident with the 
passage the wave. 

the hardness peak indicates the arrival actual water released from 
Norris Dam, requires four days, corresponding velocity about miles 
per hr, for water travel 183.3 miles from Norris Dam Chattanooga, 
compared with (3.5 miles per hr) for the wave. The flow the river 
Chattanooga before the arrival the wave was 23,000 per sec. This 
would indicate that the actual water the wave when arrived Chicka- 
mauga Dam site was the river front the original wave developed 
Norris Dam. 


TRANSLATORY WAVES WHEELER RESERVOIR 


Translatory waves Wheeler Reservoir (Fig. have quite different char- 
acteristics from those tributary streams, such the Clinch River. The 


TENNESSEE 
WHEELER 
DAM 
Mallard Creek 
M300 
Whitesburg 
Decatur 
GUNTERSVILLE 
Scale Miles DAM 


Fic. RESERVOIR 


channel has greater depth and width, together with flatter water-surface 
slope. Even with abrupt changes considerable magnitude the rate 
inflow, translatory waves are quickly smoothed out. 

Continuous stage recorders located Alabama Guntersville lower 
lock gage (mile 349.0), Whitesburg (mile 333.7), Decatur (mile 305.1), and 
Mallard Creek (mile 293.4), together with hourly observations Wheeler 
Dam (mile 274.9), provide data from which elevation graphs Figs. and 
were prepared. The accompanying hydrographs are prepared from computed 
discharges Guntersville and Wheeler dams. The drainage area above 
Guntersville Dam 24,450 miles with mean annual discharge 41,400 
per sec, whereas the drainage area above Wheeler Dam 29,590 miles 
with mean annual discharge 50,000 per sec. 

Theoretical discharge Guntersville Dam, which 
developed these waves, were made quickly enough considered sudden 
for the purpose analysis. Am. Soc. E., has given 
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that applicable for computing velocities such waves: 


approximation this formula, based the assumption that the differ- 
ence between and negligible, and also that the channel uniform 
cross section, given Professor 


which the velocity wave created suddenly, the effect the differ- 
ence between and being taken into consideration; initial water 
velocity; the acceleration gravity; and are the depths the center 
pressure and respectively; and are channel cross-section 
areas occupied the wave beginning and end, respectively; the veloc- 
ity wave created suddenly, the difference between and being con- 
sidered negligible; and the mean depth water. The constant 0.682 
inserted convert miles per hour. 

The times travel observed translatory waves reservoirs are deter- 
mined from the beginning-of-rise rising waves and the beginning-of-fall 
falling waves. Observed velocities, co, and velocities and computed 
Eqs. are shown Table for each four rising and four falling waves, 
for each the four reaches between intermediate gages, and also for the over- 
all reach between Guntersville and Wheeler dams. 

Velocities Translatory Waves Wheeler computed 
Eqs. (see Table agree fairly well with the observed velocities some 
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cases and differ widely others. 


For rising waves the Guntersville Dam 


Whitesburg reach, computed velocities Eq. and Eq. both agree 
closely with the observed velocities except for the wave January 


served velocities for this wave 
are lower all reaches than 
observed velocities for the other 
three waves. this reach 
computed velocities and 
are practically the same for all 
waves. 

For rising waves 
Whitesburg Decatur reach, 
computed velocities are about 
the same for the previous 
reach, whereas computed veloc- 
ities are consistently lower. 
The average velocities, cs, 
practically the same 
average For rising waves 
the Decatur.to Mallard Creek 
reach, both computed velocities 
give about the same results, 
which are higher than observed 
velocities all cases. The av- 
erage only little more 
than half the observed velocity 
for the previous reach. For 
rising waves the Mallard 
Creek Wheeler Dam reach, 
both computed velocities again 
give about the same results, 
but are about double the aver- 
age the observed velocities. 
the over-all reach from 
Guntersville Dam Wheeler 
Dam, rising wave computed 
velocities 
lower than computed velocities 
but are about one third larger 
than the observed velocity. Ob- 
served and computed velocities 
for falling waves are similar 
those for rising waves. 

Considerable difficulty en- 
countered determining the 
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effective average channel cross section for the application Eqs. the Whites- 
burg-Decatur and the Decatur-Mallard Creek reaches. Channel cross sections 
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taken about one-mile intervals are available throughout the length Wheeler 
Reservoir. Upstream from Whitesburg the reservoir largely within the old 
river channel and has fairly regular cross sections, while downstream from Mal- 
lard Creek the reservoir covers the entire river bottom. Between gages 
there gradual transition from river channel river bottom sections. This 
may account part for the wide discrepancy between observed and computed 
velocities these two intermediate reaches, but would not account for the 
even wider discrepancy the results for Mallard Creek Wheeler Dam reach. 
probable that when the reservoir widens, the lateral spread the wave 
reduces its downstream velocity. Another probable source error the 
inability determine discharges intermediate stations. Discharges 
Guntersville and Wheeler dams are determined from gate settings and should 
reasonably accurate. Discharges intermediate stations are estimated 
from these discharges taking into account local inflow from tributaries. 


GUNTERSVILLE Dam 


WHITESBURG 
wave, 
1939 Observed Computed Observed Computed Observed 


Jan. 1.50 10.0 19.7 19.1 2.50 11.6 7.6 9.9 2.00 6.0 
May 0.75 20.0 20.8 20.4 1.50 19.3 17.4 20.9 1.00 12.0 
June 0.75 20.0 19.0 19.1 1.25 23.2 16.8 20.1 1.25 9.6 
Aug. 0.75 20.0 19.4 18.6 1.75 16.6 18.3 19.5 0.92 13.0 
Average 17.5 19.7 19.3 17.7 7.5 0.1 10.2 

(b) Wav 
Jan. 0.83 0.7 20.2 1.92 15.1 18.3 0.5 1.25 
May 1.00 15.0 20.1 20.7 1.50 19.3 20.8 1.00 12.0 
June 0.83 18.0 21.3 20.0 9.2 19.1 20.3 1.00 12.0 
Aug. 0.83 18.0 19.1 19.7 2.67 10.9 17.9 20.2 1.00 12.0 
Average 0.3 20.2 13.6 18.1 0.4 114 


general, observed velocities the reaches upstream from Decatur where 
the pool largely within the original river channel are materially higher than 
the reaches below Decatur where the pool submerges the original river banks. 
Two bridge embankments Decatur restrict the flow the original channel 
these points. When wave reaches the upstream embankment portion 
rebounds form upstream wave that reduces the magnitude the 
wave downstream points. 

Observed velocities determined from beginning-of-rise and 
points not differ materially given reach, both being about miles 
per upstream from Decatur and about miles per downstream. 

Surges Wheeler will noted Fig. that, when the 
release Guntersville Dam was reduced suddenly and then held steady for 


7 
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time, series surges occurred fairly regular intervals about hr. 


This interval about the time required for wave travel from Guntersville 
Dam Wheeler Dam and return. Usually the peaks Guntersville Dam 
coincide with valleys Wheeler Dam. These surges have much smaller 
magnitude Decatur than stations upstream downstream. Water- 
surface profiles the reservoir plotted for elevations observed successive 
hours show that the reservoir water surface appears pivot Decatur. 

The water surface falls upstream stations when rising downstream 

stations. 


TRANSLATORY WAVES THE TENNESSEE RIVER BELOW 
Pickwick LANDING Dam 


Translatory waves the Tennessee River below Pickwick Landing Dam 
(see Fig. are typical waves developed regulated releases natural 


Hour, WHEELER RESERVOIR 


MALLARD CREEK GUNTERSVILLE Dam 


CREEK WHEELER Dam WHEELER Dam 
Date 
wave, 
Computed Observed Computed Observed Computed 1939 


13.4 14.0 2.00 9.0 24. 23.0 8.00 9.3 18.4 19.1 Jan. 
17.2 14.6 1.50 12.0 23.8 23.6 4.75 15.6 19.2 19.9 May 
14.7 14.3 1.00 18.0 21.6 23.3 4.25 17.4 17.7 19.3 June 
13.8 13.4 1.83 9.8 23.5 22.7 14.1 18.5 18.5 Aug. 
14.8 1 deo 12.2 23 23.2 iad 14.1 18.5 19.2 Average 

1.50 2.0 23.2 5.50 13.5 19.4 19.6 Jan. 
18.0 14.9 1.50 12.0 20.2 23.5 5.00 14.8 20.4 20.0 May 

14.3 1.00 18.0 16.4 23.1 6.00 12.3 19.5 June 

16.6 14.0 1.25 14.4 25.3 23.0 5.75 12.9 19.3 19.3 Aug. 

16.2 4.1 22.3 23.2 13.4 19.7 19.6 Average 


river channel having considerable discharge capacity,and comparatively 
flat slope. The mean annual discharge this section the river about 
60,000 per sec, and the fall 207 miles, 0.3 per mile. 
Waves for which gage-height graphs are shown Fig. are selected 
cover considerable range flow conditions and have discontinuities that 
can identified successive downstream stations. Pertinent data lower 
Tennessee River gaging stations are given Table Gage-height graphs 
are shown for Pickwick lower lock (Tenn.), Savannah (Tenn.), Perryville 
(Tenn.), Johnsonville (Tenn.), Buchanan (Tenn.), and Gilbertsville (Ky.) 
gages whenever the records for these stations are available for the wave selected. 
Stages Gilbertsville are affected Ohio River stages all times, and high 
stages the Ohio River affect stages far upstream Johnsonville. 
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Wave Velocities velocities for the reach from Savannah Buchanan 
are computed the same formulas are used for the Clinch River waves. 
These results, together with observed velocities determined from the beginning- 
of-rise, midpoint-of-rise, and beginning-of-fall points, are listed Table 


TABLE 7.—Lower TENNESSEE RIVER STATIONS 


Station 


Drainage area 
(see Fig. 


River mile miles) 


134.9 325.27 34,500 
Johnsonville 100.5 319.61 38,520 
21.7 299.56 40,200 


Elevation gage zero above Mean Gulf Level. 


Discharge ratings, area, and channel width curves for each station are given 
Table The Pickwick lower lock Savannah reach was excluded from 
the velocity comparison because the channel and discharge character neither 
these stations are representative the intervening reach, and the Buchanan 
Gilbertsville reach was excluded because points very few waves can 
identified satisfactorily Gilbertsville due the effect the Ohio River. 
The method applying formulas for wave velocities similar that pre- 
viously described for waves the Clinch River. 
Computed wave velocities and for the beginning-of-rise point are 
generally lower than the observed velocity the average about 
1.5 times and 1.4 times However, the variation too large con- 
sidered consistent relation. For the midpoint-of-rise, computed velocities 
general agree much better with the observed. There also better agree- 
ment between and than computed velocities for the beginning-of-rise 
points. For beginning-of-fall points, the agreement between computed and 
observed velocities about the same for the midpoint-of-rise. 
study velocities midpoint-of-fall shows agreement between 
observed and computed velocities, the computed velocities being much higher 
than the observed. When falling wave travels down the lower Tennessee 
River, its shape gradually modified until finally becomes recession that 
follows the characteristic curve water draining out the channel. 
Factors Affecting Results—During periods rapidly changing stage, such 
exist during the passage wave, the slope the water surface and the 
discharge differ considerably from the steady flow conditions. Due the 
flat slopes that normally exist the lower Tennessee River, these variations 
have relatively greater effect the discharge than the Clinch River 
which has steeper slope. This factor accounts part for results being less 
satisfactory for the lower Tennessee River than for the Clinch River, and 
would most pronounced the beginning-of-rise where the poorest agree- 
ment was found. 


65 


| 


June, 1944 TRANSLATORY WAVES 873 


Rising waves the lower Tennessee River have tendency flatten and 
lengthen successive downstream stations contrast the tendency 
Clinch River rising waves maintain their height and shorten. Falling 
waves both rivers have tendency recede characteristic pattern 
which depends the volume the river channel. 

The wide variation and lack consistency observed velocities are due 
partly inability identify selected points the wave with any degree 


TABLE 8.—TRANSLATORY WAVES THE TENNESSEE RIVER FROM 
SAVANNAH (MILE 189.9) BucHANAN (MILE 66.2) 


1938 
May 5.2 10.6 6.0 3.9 4.2 9.5 12.8 
May 4.8 12.5 6.1 4.0 4.5 10.3 13.8 
May 7.7 11.7 5.4 4.2 4.6 13.6 15.5 
July 16.4 6.8 4.3 4.8 11.3 15.0 
July 3.1 22.2 4.5 3.7 4.9 10.0 21.5 
4.1 6.7 6.4 3.7 3.9 8.0 9.0 
Jan. 9.1 12.9 4.6 5.1 13.4 14.2 
Jan. 6.2 16.4 7.2 4.5 5.1 14.5 22.1 

For 

1938 
May 5.2 10. 4.4 4.4 4.2 9.5 12.8 
July 3.1 22.2 5.0 4.9 10.0 21.5 
Dec. 4.1 4.2 3.8 3.9 8.0 9.0 

1939 
Jan. 6.2 16.4 4.2 5.0 5.1 14.5 22.1 

1988 
May 10.6 4.8 5.1 4.6 4.3 12.8 10.3 
May 12.5 7.8 3.9 4.8 4.9 13.8 13.6 
July 16.3 9.4 4.5 5.0 5.1 15.0 13.8 
6.7 4.1 3.9 4.0 3.9 9.0 8.1 
Jan. 12.9 6.2 5.1 4.8 5.0 14.2 14.0 
Jan. 17.4 10.4 8.7 5.4 5.2 22.1 17.0 


certainty, particularly the two downstream stations. The midpoint-of-rise 
can determined with the most certainty when the wave consists one 
continuous rise; but unfortunately this occurs very few waves. 


CHARACTERISTICS TRANSLATORY WAVES 


When rising wave developed fairly uniform increase discharge, 
the case most the Clinch River waves, the rising face the wave 
tends retain the initial rate increase discharge successive down- 
stream stations. When the wave develops empty channel, Fig. 
2(c), the beginning-of-rise continues abrupt successive stations; but 
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when the wave develops considerable initial flow, Fig. 2(a), fillet 
develops the beginning-of-rise which enlarges the wave travels down- 
stream. 

similar fillet occurs the wave approaches fully developed rise. the 
wave progresses downstream these two fillets continue enlarge until the 
upper limb the beginning-of-rise fillet meets the downward limb the end- 
of-rise fillet. Although these fillets are enlarging, the duration rise the 


TABLE CHANNEL AND 


SAVANNAH 
area miles per area (miles per area 
(sq ft) hr) (sq ft) hr) (sq 

(1) (2) (3) (4) (6) (7) (8) 

2 6,000 3,620 2.2 14,000 10,200 3.4 ak 
12,800 5,800 2.5 21,000 11,500 3.8 11,900 18,00 
21,500 8,000 2.8 28,750 12,700 4.1 18,200 19,200 
8 31,500 10,300 3.1 37,500 14,000 4.5 25,600 20,400 
10 42,400 12,500 3.3 46,500 15,300 4.8 34,200 21,800 
54,400 14,910 3.6 55,500 16,700 4.9 43,600 
66,400 17,200 3.8 65,500 18,100 5.0 54,000 
78,400 19,570 3.9 75,500 19,500 5.0 64,700 25,700 
90,600 21,900 4.0 86,000 21,000 5.0 76,000 
102,800 24,290 4.1 97,000 22,400 5.0 87,600 28,400 
115,200 26,800 4.2 108,000 23,800 4.9 100,000 29,800 
127,600 29,060 4.3 119,000 25,300 4.9 113,000 31,200 


initial rate becomes shorter and finally disappears. the Clinch River (see 
Fig. this seems occur some point above Wheat, miles below Norris, 
and the lower Tennessee (see Fig. point above Perryville, miles 
below Pickwick. 

the case falling wave, characteristic recession tends develop 
regardless the initial rate reduction discharge. This recession re- 
tained successive stations the Clinch River until fillets the beginning- 
of-fall and the end-of-fall enlarge cover the entire falling face the wave. 
This again appears occur the time that waves reach Wheat Perryville. 

When the period maximum discharge the crest wave continues 
for only short period, Fig. 2(c), the spreading fillets the end-of-rise 
and the beginning-of-fall gradually shorten this period the wave pro- 
gresses downstream. Eventually the two fillets will overlap and the discharge 
will begin fall off before the full rise has been developed. When this occurs 
flattening the wave begins and the magnitude rise becomes progressively 
smaller. valley between two rises affected similar manner, resulting 
gradual reduction the depth the valley shown Figs. 
and 2(f). 


Velocities translatory waves, both rising and falling, can approxi- 
mated with either Eq. Eq. natural channels such the Clinch and 
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lower Tennessee rivers when rated stations representative sections are 
located suitable intervals throughout the reach. 

Midpoints rise fall are best suited for determining the velocity 
observed wave. 

Eqs. and give similar results for velocities translatory waves 
Wheeler Reservoir, but these results are consistently higher than observed 
velocities. 


JOHNSONVILLE BucHANAN 
Gage 
(miles per (cu area (miles per area miles per 
per (cu per 
hr) sec) (sq ft) hr) sec) (sq ft) hr) 
(10) (11) (12) (13) (14) (15) (16) (1) 
3.4 21,250 13,200 3.2 7,000 16,300 ‘a 4 
4.0 31,750 15,200 3.6 14,500 18,200 2.8 
4.5 44,000 17,300 3.9 23,500 20,000 3.2 
7,000 19,400 4.2 33,000 21,900 3.7 
5.2 70,000 21,600 4.4 44,500 23,800 4.0 
5.4 84,000 23,800 4.5 56,500 25,800 4.3 
5.6 98,000 25,900 4.7 69,500 27,800 4.6 
5.8 112,000 28,100 4.9 83,000 29,900 4.7 
6.0 126,000 30,200 5.0 96,000 31,800 
6.2 142,000 32,500 5.1 112,000 33,900 5.0 
6.3 158,000 34,800 5.1 127,000 36,000 4.9 


Within the range these data the volume flow appears have little 
effect the velocity translatory waves reservoir. 
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AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 1852 


PAPERS 


WEDGE-SHAPED STRUCTURAL MEMBERS 
UNDER DIRECT STRESS AND BENDING 


Presented this paper are expressions for computing compressive and 
shearing unit stresses the concrete, tensile unit stresses the reinforcement, 
and bond unit stresses between the and concrete for wedge- 
shaped reinforced concrete members subjected direct stress and bending. 
Curves expedite the use the expressions are also presented. Only 
approximate solutions are offered since the materials not follow strictly the 
assumptions made the derivations. 


INTRODUCTION 


The subject stresses wedge-shaped reinforced concrete members sub- 
jected direct stress and bending has been treated relatively few papers. 
presented derivation pertaining the shearing unit stress for the 
case which the compressive surface normal the section which the 
stress considered acting, and the reinforcement inclined that section. 
treatment serves basis for the more extensive derivations 
expressions for shearing unit stress and bond unit stress found this paper. 
presented expressions with which compute the shearing unit 
stress any point along section the member, especially the compressive 
side the neutral axis. However, the shearing unit stress used design 
measure diagonal tension, and for this reason only the stress the tensile 
side the neutral axis interest. the compressive side the neutral 
axis the actual shearing unit stress (maximum) would found equal 
one half the maximum compressive unit stress and would acting planes 
inclined 45° the compressive surface. 


comments are invited for immediate publication; insure publication the last dis- 
cussion should submitted November 1944 


Structures Consultant, Goodyear Aircraft Corp., Akron, Ohio. 
Eisenbetonbau,” Konrad Wittwer, 1929, Vol. Pt. 290. 


Schubspannungen bei aussermittiger Beanspruchung gerader mit Quer- 
Beton und Eisen, Vol. 37, 1938, 137. 
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symbols used this paper are defined where they first 
appear, either text diagram. They conform essentially American 
Standard Letter Symbols for Mechanics, Structural Engineering and Testing 
Materials‘ and American Standard Letter Symbols for prepared 
committees the American Standards Association, with Society representa- 
tion, and approved the Association 1932 and 1942, respectively. 

Assumptions.—The usual assumptions that are made deriving the ordi- 
nary expressions reinforced concrete design are: 


(1) section plane before bending remains plane after bending; 

(2) The reinforcement resists all the tension; 

(3) The unit stress the concrete proportional the distance from the 
neutral axis and the reinforcement times the proportionality for the 
concrete, which the ratio between modulus elasticity steel and 
modulus elasticity concrete; and 

(4) Perfect adhesion exists between the concrete and the reinforcement. 


BENDING STRESSES 


Fig. shows part wedge-shaped reinforced concrete member acted upon 
external force, This force can resolved into components 
(normal) and (shear) shown. The total compression the section 


indicated and the component, normal the section the total 
tension the reinforcement The compressive unit stress normal 
and the edge the section shown fo. 

Fig. assume that plane A—A’ and bound part wedge-shaped 
concrete member. The dashed line the side A’B’ represents 
the reinforcement. Assume that plane held position while plane 
rotated through very small angle The deformation, change 


length, the reinforcement for length tu. Using the law sines one 
may write 


ASA—Z10a—1932. 
ASA—Z10.3—1942. 
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Eqs. and the distance from the compression surface the member 


90+a-8 
6 


--B' 


BETWEEN ANGLE CHANGE BENDING AND DEFORMATION STEEL 


the center the tensile reinforcement measured the plane section 
and the ratio the distance from the neutral axis the outer compressive 
fiber reinforced concrete member the distance, 


tan 


BETWEEN DEFORMATION CONCRETE AND DEFORMATION STEEL 


Let Fig. represent the deformations length Ay. The angle may 


expressed terms the compressive unit stress the concrete, fo, normal 
section 


(3) 
and, multiplying both sides Eq. k), Eq. becomes 


Now 
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and 


Therefore, the unit strain the reinforcement can expressed 


From consideration equilibrium (namely, the vertical 
forces equal zero), the following expression may written: 


(9) 


which the width the member; and the ratio the area the 
tensile reinforcement the effective area the concrete, bd. The area 
the reinforcement taken plane normal the axis the reinforcement. 

Upon taking moments, about the point application the normal com- 
ponent the resultant tensile stress the reinforcement, the forces acting 
the section found that 


Solving for Eq. and substituting Eq. 


Using Eq. and Eq. and eliminating one can write 


Eq. can rearranged and solution for computed follows: 


and the unit stress the reinforcement 
f. 8 8 E. k f k ( ) : 
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1.40 


2. 


cos 6 
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Eq. has been used plot the curves for various cos shown 


Fig. The values are the abscissas and the k-values the ordinates. 


The following expressions are presented for use computing the maximum 
normal stresses: 


and 


Eq. 16a will noted that the maximum compressive unit stress the 
concrete given This expression® for the unit stress 
plane normal the compression face the member. 

the case tee-shaped section which the flange compression, 


the flange. One may then use the following approximate expressions when 
computing the unit stresses due direct stress and bending: 


and 


which width the compressive flange; and area effective tensile 
reinforcement, for this case width 


SHEAR AND STRESSES 


Imagine Fig. part reinforced concrete member with variable 
cross-section area. Assume that two sections are passed parallel the hori- 
zontal plane distance from one another. Taking moments about the 
intersection the resultant the compressive stresses with the section, 


which: resultant the external forces;r perpendicular distance from 
line action the resultant, the moment center; total tensile stress 
the reinforcement; and perpendicular distance from the line action 
the tensile force the reinforcement the moment center. 


Pressure, Walls and Bins,” William Cain, John Wiley Sons, Inc., New York, 
Ed., 1916, 240. 
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Considering changes take place only and the change tensile 
stress the reinforcement distance may expressed 


omens 


The following expressions for dr, dz’, and may derived from the geo- 
metrical relations (see Fig. 5): 


and 
Substituting Eqs. Eq. 19, the expression for 
The differential tension, the produces unit shearing 
stress, the adjacent concrete, which parallel the reinforcement, and 
equal unit shearing stress right angles the reinforcement. 
the width the member, oné may write 


Substituting Eq. Eq. and simplifying: 


) 

oN 
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From Fig. also noting Eq. 18, may seen that 


and 


or, the moment the external forces with respect the selected moment 


26) 


which the distance, the plane the section, between the resultant 
the compressive stresses and the resultant the tensile stresses. 
desired use the moment about the reinforcement was done 
previously, the expression for becomes 


and substitution Eq. results 


the case tee-shaped cross section, the tee the compressive 
side, the position the resultant the compressive stresses changes slightly 
with respect the tee that tan may taken equal tan without 
appreciable error. 29a can then written 


which the width the stem the tee-shaped member. 

The expression for the unit bond stress, may written 


cos 
since, from Eq, 22, one has bv. For rectangular sections, 


Upon substitution Eqs. Eq. 23, the expression for becomes 

cos 


June, 1944 WEDGE-SHAPED MEMBERS 885 


29a, the expression for becomes 


Eq. may used for tee-sections substituting for the last quantity. 


tan 


approximate expression for may derived the following 
manner: Since varies slightly with especially over distance Ay, one may 
write (noting Fig. 6), 


Dividing the expression for tan becomes 


ILLUSTRATIVE EXAMPLE 


Assume that desired determine the unit stresses section 
the part the structure shown Fig. force, per foot member 


N=30000 Lb 


d=48 In. 
e=72 In. 


normal the plane the paper, resolved into vertical (N) and horizontal 
(V) components. Other data are The computation procedure 


> 
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From Eq. 16a, 596 1.072 639 per in.; and, from Eq. 


The shearing unit stress (as measure the diagonal tension) 


can computed using Eq. 29a; thus, 


(0.903 persqin. The bond unit stress can 


Formulas have been developed for compressive bond and shearing stresses 
the concrete and for tensile stress the reinforcement for wedge-shaped 
reinforced concrete members. These expressions are applicable members 
subject bending and direct stress. The formulas are correct far the 
assumptions used the ordinary flexural theory are applicable reinforced 
concrete. 
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Founded November 1852 


REPORTS 


PRINCIPLES DEPRECIATION 


REPORT THE SPECIAL COMMITTEE AUTHORIZED 
THE BOARD DIRECTION 
ANALYZE AND DISCUSS 
1943 REPORT THE NATIONAL ASSOCIATION 
RAILROAD AND UTILITIES COMMISSIONERS’ 
COMMITTEE DEPRECIATION 


FoREWORD 


For many years the Society has supported busy schedule studies 
improve and crystallize the best practice the application depreciation 
theories related civil engineering projects. 1911 Special Committee 
“To Formulate Principles and Methods for the Valuation Railroad Property 
and Other Public was appointed. presented Progress Report 
1914, and 1916 brought out extensive Final Report. This was widely 
discussed and thus was developed into final form appears 310 pages 
the 1917 Transactions. This treatment was widely acclaimed setting 
fundamental engineering principles applied such matters maintenance and 
operating costs, amortization, depreciation, obsolescence, valuation, and similar 
elements. The Society’s interest has since continued, expressed through 
number committee and Division activities. 

Recently there appeared 277-page printed report Committee the 
National Association Railroad and Utilities Commissioners (NARUC) 
which deviated from some the long-standing principles previously enunciated 
Society publications and general use. Having knowledge this report, 
the Board Direction requested that opportunity present the contrary 
views civil engineers afforded. President Pirnie then appointed Special 
Committee with instructions study the report objectively and present its 
findings the Board Direction the Society. order that this committee 
might submit its findings, the NARUC committee granted extension 
time April before taking further official action. 

The tentative report the Society’s Special Committee was submitted 
the Board its St. Louis Meeting, April 24. After considering the findings, 


forward all comments this report Society head ters, West 39th 
New York, Y., for transmission the 
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the Board instructed that this document printed for the benefit all 
members. presented herewith. summary quoted from pages 
XIX appended. 

still further development that the subject valuation, depreciation, 
and rate making made special item technical consideration the 
time the Annual Convention the Society Cleveland, Ohio, July 21, 
1944. view the active diversity views occasioned the NARUC 
report, especially urged that all who can submit written comments 
for the study the Special Committee prior that meeting. Such comments, 
forwarded headquarters, will promptly transmitted the committee. 


REPORT THE SPECIAL COMMITTEE 


Special Committee the Society was appointed February, 1944, for 
the purpose analyzing and commenting upon the report the Committee 
Depreciation the National Association Railroad and Utilities Com- 
missioners (NARUC) issued 1943. The Society’s Special Committee has 
studied that report and respectfully submits its analysis and 

The NARUC report embraces some 277 pages. The Special Committee 
the Society has deemed best limit its report discussion only the 
most fundamental differences which has with the report the NARUC 
committee. The principal divergencies views are follows: 


The Society’s committee takes exception the fundamental concept 
depreciation set forth Conclusion the summary the NARUC com- 
mittee report (Appendix). 

The NARUC report states (Conclusion 26, Appendix): 


“The straight-line method generally for public utility 
accounting and financial purposes and also for the computation both 
depreciation expenses and accrued depreciation for purposes rate making.” 


and again (Conclusion 21, Appendix), 


“The straight-line method very general use and the ascendancy 
the public utility field.” 


(Conclusion 22, Appendix), 


“The method simpler than, and not seriously affected 
errors estimates service lives the interest methods.” 


The NARUC report repeatedly states that actual depreciation conforms 
the depreciation reserve accrued the straight-line method. The following 
excerpts are illustrative (page 169, NARUC report): 


“The depreciation reserve shows the results date the effective 
depreciation 


(on page 170, NARUC report), 
cost the gross base, then the reserve should deducted from 


arriving the rate base, for the reserve alone measures actual de- 
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and (Conclusion 39, Appendix), 


properly computed depreciation reserve the best measure 
accrued existing depreciation, since such reserve reflects that part the 
cost the property service which relates the exhausted expired 
economic service 


The NARUC report states (Conclusion 24, Appendix): 


The total dollar charges customers the long run, however, 
will greater under the sinking-fund compound-interest methods than 
under the straight-line method.” 


numerous places the NARUC report reference made the 
necessity correcting depreciation reserves utilities when such reserves 
differ from those resulting from the application the straight-line method, 
and the recommendation that the deficiency made through with- 
drawals from surplus special section the income account, some other 
source. This illustrated the following quotation (Conclusion the 
Appendix): 


“The depreciation reserve should not readjusted gradually through 
modification the annual depreciation rates when the difference between 
the book reserve and the proper reserve substantial. When the differ- 
erence not substantial satisfactory spread the remaining net cost 
the properties over their remaining 


The foregoing examples and other features the NARUC report are dis- 
cussed more fully subsequently. 


DEFINITION DEPRECIATION 


The Society’s committee believes that the many definitions deprecia- 
tion that Mr. Chief Justice Hughes the Lindheimer case (292 U.S. 151), 
follows: 


“Broadly speaking, depreciation the loss, not restored current 
maintenance, which due all the factors causing the ultimate retirement 
the property. These factors embrace wear and tear, decay, inadequacy 
and obsolescence. Annual depreciation the loss that takes place 
year. 


and that the existing NARUC uniform System Accounts reading: 


applied depreciable utility plant, means loss 
service value not restored current maintenance incurred connection 
with the consumption prospective retirement utility plant the 
course service from causes which are known current operation 
and against which the utility not protected insurance. Among the 
causes given consideration are wear and tear, decay, action the 
elements, inadequacy, obsolescence, changes the art, changes demand 
and requirements public authorities *.” 


—are the best. 
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contrasted with these definitions, which reflect the latest authoritative 
statement the Supreme Court, the 1948 NARUC depreciation com- 
mittee report describes depreciation Conclusion (Appendix) follows: 


“The significance the depreciation phenomenon may described 
follows: 


part, the service life, capacity, utility property resulting 
from the action one more the forces operating bring 
about the retirement such property from service; 

The forces operating include wear and tear, decay, action 
the elements, inadequacy, obsolescence, and public require- 
ments; 

Depreciation results cost service.” 


and again Conclusion (Appendix), 


the using the economic service life asset 
regardless the specific factor which may ultimately cause the asset 
retired. the exhaustion service life, not the particular cause 
retirement, that important.” 


Attention directed the foregoing definitions and the fact that the 
definitions heretofore effect and approved the Supreme Court refer 
service value not restored current which the 
new committee definition referred the “consumption whole 
part the service life capacity utility This change which 
effect predicates the deduction from cost new upon service 
life capacity utility” rather than upon loss service value not restored 
current maintenance. 

The Society’s committee cannot subscribe the implied opinion the 
NARUC committee that remaining service life and remaining service value are 
synonymous, although there usually relationship (not direct, however) 
between the two. 

There are other definitions and terms used the NARUC report which 
the Society’s committee takes exception, but which would seem unneces- 
sary devote time and space this discussion. 

would appear serve useful purpose the discussion the NARUC 
report refer specifically many minor points divergence with the theories 
outlined the NARUC report their application the mathematical 
examples with which the Society’s committee not accord, the unrealistic 
conclusions which the application the principles the NARUC report 
would lead many instances. 


DEPRECIATION 


NARUC report discusses the several aspects depreciation and 
depreciation accounting length. the opinion the Society’s committee 
fails present adequately the fundamental conception that loss value 
property through age and use inescapable and that the extent that loss 
value any date should the basis all accounting and financial policies 
with reference depreciation. 


q ( 
‘ 
< 
| | 
| 
| 
q 
| 
4 
| 
‘ 
‘ 
q 
1 F 


June, 1944 DEPRECIATION 891 


the last analysis the estimate depreciation for any period must rest 
human judgment. The most reliable estimate depreciation one made 
qualified persons who have examined the particular property, studied the 
mortality characteristics the kinds units which composed, and 
considered the prospects future use each the units (or groups like 
units when there are included the property large numbers these). 

Moreover, the property ought re-examined and the depreciation rede- 
termined appropriate intervals order take into account any changes 
the probable future usefulness any the units. 

Accrued depreciation existing any particular property any particular 
time matter fact and not theory. True cannot determined with 
mathematical accuracy formula, but subject determination within 
reasonably close limits men experienced the design, construction, and 
operation utility properties. basically engineering determination 
and the Society’s committee contends that absolutely impossible arrive 
any reasonable answer any other process, and certainly not any 
theoretical calculation. 

The NARUC committee report goes considerable length discredit 
what refers depreciation,” and which explains that 
depreciation which can seen the naked eye. one competent 
determine depreciation contends that depreciation can seen with the naked 
eye. The modern methods determining depreciation include some observa- 
tions the property, but they include much more. This modern method 
might properly called the engineering method estimating actual de- 
preciation. 

The engineering determination accrued depreciation involves many 
operations and analyses addition what has been referred 
it. Operating records are studied, including past and probable future growth 
the system, the types equipment, their functional features such adequacy 
and obsoleteness, and many other elements that would affect probable future 
usefulness. Poles are inspected and butts measured determine their actual 
condition; maintenance records are examined determine the care that has 
been given the upkeep the property; deferred maintenance, any, 
records past retirements are and fact all data having 
bearing the future usefulness the units involved the particular property 
question are studied. 

When the depreciation property has been determined the afore- 
described manner, the actual depreciation thus established factual and affords 
sound basis for appropriate financial and accounting procedures. The 
Society’s committee contends that accrued depreciation can only ascertained 
accurately this method. 


STRAIGHT-LINE DEPRECIATION 


The NARUC report repeatedly expresses preference for the 
method depreciation accounting, ostensibly the grounds simplicity. 
The Society’s committee emphasizes that actual depreciation factual and that 
accounting, correct, must reflect the facts actual depreciation. The 
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depreciation the property—not the accounts. Therefore, inad- 
visable that any method accounting for actual depreciation should adopted 
for other reason than that simple. 

exception taken the use age-life methods help estimating 
the annual charges for depreciation where sufficient retirement data are available 
from any particular utility system enable past useful lives retired property 
determined with fair degree accuracy. Accurate retirement data that 
could considered mass data cover relatively small number the numer- 
ous units involved public utility property. Also, the data that are available, 
generally, cover items that are different and usually those presently 
existing the usual utility property. 

example, most retirement data poles not reflect the service that 
can gotten from pole having modern preservative treatment. Retirement 
data water mains necessity cover pipe installed and used during period 
when modern coatings were not available, and modern practice water treat- 
ment correct for corrosion and incrustation was not yet developed. Electric 
meters that have been retired for other than functional reasons were old types 
developed the early years the electric industry. situation exists 
for gas and water meters, and most the equipment and material all utilities. 
Included the available retirement data are some modern items that have been 
retired due casualty functional reasons, but such data certainly cannot 
considered more than very roughly indicative what may expected the 
future. 


Estimates future useful lives made with the help pertinent retirement 


data may considered indicative where the data are sufficiently compre- 
hensive. The available data may meager nonexistent many the 
items involved, and exactly similar units will have varying life histories 
different properties. large percentage retirements are made for func- 
tional reasons, such inadequacy, obsolescence, action public authorities, 
etc., manifest that factors such the character the community served, 
its rate growth, the cost fuel, labor and material, have definite bearings 
usefulness items involved utilities serving them. operating unit 
that, due size operating costs, would retired one locality, may have 
several times many years useful service some other system. 

Incidentally, the Society’s committee challenges the contention simplicity 
because the data employed computations connection with actual deprecia- 
tion are commonly taken from tables; and, any case, inability make calcu- 
lations simple compound interest should not sufficient grounds for using 
depreciation accounting method that distorts the facts for certain classes 
property. 

That method accounting for depreciation for any particular unit 
property best which most closely harmonizes with the factual accrued depreci- 
ation the unit. general, depreciation emerges from the twilight zone 
maintenance for very short lived units, which are usually financed directly 
through operating expenses, the very small annual requirements units 
property having long lives. With the short lived types (10 years less) 
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not important what annual depreciation method used, but with the long 
lived units the importance greatly increases. 

Straight-line depreciation does not conform the facts for long life property 
whose accrued depreciation experience more nearly parallels the accumulation 
sinking fund. 

Depreciation reserves are accrued meet future obligation and where that 
period long the use the straight-line method unsound economically. 
With long life property, such dams, main line railroad embankments, cast- 
iron water works distribution systems, masonry structures, etc., manifestly 
sound business provide for the obligation maturing that remote time 
annuity rather than uniform allowance straight-line basis. The 
actual depreciation which the long lived unit has suffered will usually much 
less than the accumulation straight-line reserve. 


DEPRECIATION NOT CHEAPEST FOR CUSTOMERS 


The NARUC report purports show that the method 
depreciation accounting will require the lowest contribution from the customers 
any the methods considered. This manifestly incorrect. The cus- 
tomers’ contribution entirely independent the method accounting 
followed. Moreover, all methods depreciation accounting will result 
identical contributions from the customers applied correctly with due regard 
the time value money. Therefore the statement Conclusion 
(Appendix) that the 


greater under the sinking-fund compound-interest methods than under 
the straight-line 


incorrect and misleading and has basis other than the assumption that the 
rate payer would less provident getting interest his money than the 
utility will investing the sinking fund. the depreciation accounting 
correct, the total charge against operation for any unit property its value 
new less net salvage value. 


DEPRECIATION RESERVE 


The depreciation reserve created credits from annual charges against 
operations for depreciation. Retirements are charged the reserve, and the 
balance any one time should least equal the amount accrued de- 
preciation existing the same time. The reserve should normally exceed the 
accrued depreciation reasonable margin, the amount depending upon 
circumstances existing the particular utility. should large enough 
absorb casualties (not covered insurance) which are properly not included 
the determination accrued depreciation. 

The things circumstances that bring about depreciation and retirement 
property units, stated previously this report the Society’s committee, 
vary each property and from year year any one property. entirely 
feasible have the amount accrued depreciation determined regular 
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intervals varying from years, depending upon the property involved, and 
found that the balance the depreciation reserve either too high 


too low, adjustment can made the annual depreciation allowances 


the two into proper relationship. Due the nature the reserve (including, 
should, casualties and other elements that have place measurement 
accrued depreciation) there can justification for deducting accrued 
depreciation. 

submitted that the actual depreciation should control the size the 
reserve, rather than that the reserve should determine the actual accrued 
depreciation. 

The depreciation reserve created take care retirements, and seems 
proper that property retired should retired the cost recorded the books 
and records the property. The reserve created from charges operations, 
depreciation undoubtedly incremental part the cost service. The 
Society’s committee believes that the reserve definite liability and 
sense trust fund, held until such time replacements are necessary, and that 
the case continuing utility should not considered return capital 

Uniformity treatment and administrative convenience dealing with 
depreciation would desirable equitable; but neither uniformity nor 
administrative convenience should control over the facts particular situa- 
tion. The NARUC report under discussion would make abstract theory, rigid 
uniformity, and administrative convenience controlling. The Society’s com- 
mittee believes that the results such action would ruinous the public 
utility business. already stated, actual depreciation existing physical 
property public utility fact. the kind fact that cannot de- 
termined merely from set books tables, any other means than study 
the ground qualified, trained engineers operators with engineering 
training who follow approved methods. 

Abuses the engineering determination accrued depreciation are more 
probable than manipulation the basic figures upon which theoretical age-life 
determinations are made. The possibilities improper use any method 
should not used justification for abandoning produces the correct 
answer. 


RETROACTIVE ADJUSTMENT THE DEPRECIATION RESERVE 


The the NARUC Committee Depreciation recommends the 
setting depreciation reserves the basis straight-line requirements. 
further recommends that where the reserve computed the straight-line 
method substantially exceeds the reserve currently the property that the 
adjustment necessary bring the current reserve the depreciation reserve 
computed the straight-line basis should through charge surplus 
special section the income account. 

The NARUC report recognizes that this adjustment presents many practical 
difficulties. The Society’s committee feels that view the fact that practi- 
cally all utilities have been subject governmental regulation over period 
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two four decades and the reserves currently carried the books the 
utilities have been accumulated, credited and charged compliance with the 
regulation such bodies, the adjustment the current reserve upward 
equal those computed the straight-line basis withdrawal from surplus 
otherwise not only generally financially impracticable but also ethically 
unjust. 

Present reserves should,.and usually do, reflect rates charged and collected 
from consumers. arbitrarily increase the reserves this late date 
withdrawal from surplus require the utility rather than its consumers 
make good the cost property used their service. This procedure 
conflict with the basic principle depreciation well plain justice. 

The reserve requirements should measured the actual needs the 
utility and not arbitrary which the case 
many the best utilities the United States would require reserves sub- 
stantially double those which long background experience under regulation 
has determined adequate. 

the opinion the Society’s committee that Conclusion (Appendix) 
ambiguous, being variously interpreted, and impracticable equitable 
application. consideration many matters other than depreciation. 


GENERAL 


this report the Society’s committee has endeavored cover only few 
the most important features the NARUC committee report. That report 
contains many valuable data, and can considered definite contribution 
the subject depreciation. The Society’s committee has given the report 
careful study and has reached the conclusion that desirable may be, from 
the standpoint regulatory simplicity, impossible determine deprecia- 
tion, either annual accrued, means theoretical actuarial calculations 
because the many variables existing the elements causing depreciation. 
would fine crystal ball could used, and the future accurately de- 
termined, but this has not been possible date and, far can foreseen, 
never will be. 

Depreciation, its causes and effects, are pretty generally understood and 
should accepted phenomenon differing every property, and not subject 
mathematical standardization. 


Respectfully submitted, 

Howson 

Henry Chairman 


Special Depreciation Committee the Board Direction 


April 15, 1944 
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APPENDIX 


SUMMARY THE REPORT THE COMMITTEE DEPRECIATION, 
ASSOCIATION RAILROAD AND 
UTILITIES COMMISSIONERS 


(Pages x1x the NARUC Report) 


The more important conclusions and recommendations contained this 
report are summarized briefly below afford convenient general preview 
the substance the report whole. statements obviously 
not possible capture all the refinements meaning all the particular 
qualifications set forth the main text. Therefore these statements are 
taken condensed expressions the ideas and principles that are developed 
more fully the report. 


DEPRECIATION 


The modern view depreciation emphasizes the concept that the purchase 
capital goods essence purchase future services. 

The significance the depreciation phenomenon may described 
follows: 


(a) Depreciation the expiration consumption, whole part, 
the service life, capacity, utility property resulting from the 
action one more the forces operating bring about the re- 
tirement such property from service; 

(b) The forces operating include wear and tear, decay, action the 
elements, inadequacy, obsolescence, and public requirements; 

(c) Depreciation results cost service. 


Depletion the using wasting assets—natural resources, such 
timber, oil, natural gas, etc. Thus closely related depreciation 
taken unit production basis. 

Amortization applied fixed assets the process writing off the cost 
intangibles, such patents, copyrights, leaseholds, over period 
time. Such assets frequently have contractual legal life over 
which the cost should spread. 

Depreciation accounting the process charging the book cost de- 
preciable property operations over its life. 

The physical, functional, and accounting processes defined the four 
preceding paragraphs are often referred collectively, for convenience, 
depreciation. 


the past several methods have been employed account for the con- 


sumption the service economic life plant assets, including the 
retirement, replacement, appraisal, arbitrary writedown, retirement 
serve, and depreciation methods. Only the depreciation method survives 
general use. 

There are two interrelated aspects depreciation: annual expense and 
accrued depreciation. Annual expense shown the income statement; 
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10. 


12. 


13. 


15. 


16. 


17. 


accrued depreciation recorded the depreciation reserve, balance 
sheet account. 

Depreciation the using the economic service life asset re- 
gardless the specific factor which may ultimately cause the asset 
retired. the exhaustion service life, not the particular cause 
retirement, that important. 

Depreciation should not confused with efficiency. machine may 
highly efficient when near the end its service life and when almost 
fully depreciated fact. 

While the progress depreciation may retarded proper maintenance, 
the view that maintenance prevents unsound. Regardless 
the degree maintenance, economic life cannot prolonged in- 
definitely. 

The depreciation reserve measures that part the cost plant still 
service which has been written off, usually operating expense. 
the depreciation reserve has been properly determined, measures the 
accrued depreciation. The depreciation reserve should preferably 
shown the asset side the balance sheet deduction from de- 
preciable plant. 

The depreciation reserve account contra the plant account. 
not Depreciation accounting normally results the retention 
the business assets which otherwise might distributed divi- 
dends. Ordinarily these assets are not earmarked and, with minor 
exceptions, there reason why they should be. 

The financing replacements not the purpose depreciation accounting. 
Its purpose record cost operations the pro tanto cost 
property consumed therein, thus maintaining the integrity the 
investment whether not replacement occurs. Where depreciation 
computed base higher than the actual cost properties, accretions 
capital result; where the depreciation base lower than the actual 
investment, capital impairment follows. 


DEPRECIABLE AND PROPERTY 


practice distinction made between depreciable and non-depreciable 
assets. Non-depreciable property is, fact, only that property which 
has unlimited life service. Mere longevity the difficulty pre- 
dicting service life should not prevent provision for depreciation ac- 
cording the best information available. 

Overheads—including interest during construction—and indirect costs, 
the extent associated with depreciable plant, represent depreciable costs. 
Such overhead and indirect costs should assigned units property 


that the full costs the property units can spread over their re- 
spective service lives. 


Depreciation should not based upon fair value replacement cost. 
this respect the doctrine followed United Railways Electric Co. 
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20. 


21. 


22. 


23. 


24. 
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West (280 234, 253, 1930) misconstrues the economic nature 
depreciation and would administratively impossible reasonably 
accurate application. 

Depreciation expense should based upon the cost depreciable assets. 
The actual cost utility properties acquired operating units 
systems may differ from their original cost. Under most current systems 
accounts for utilities has been established that depreciation shall 
based original cost. 

Any difference between the cost and original cost less depreciation 
recorded the acquisition adjustment account. Since the proper dispo- 
sition amounts this account must depend the circumstances which 
gave rise them, specific rules for their accounting treatment are 
recommended. 


METHODS 


Several methods have been used compute public utility depreciation, but 
only three are generally recognized accepted today. These are the 
straight-line, compound-interest (“‘modified sinking-fund’’), and sinking- 
fund methods. 

The straight-line method very general use and the ascendancy 
the public utility field. 

The straight-line method simpler than, and not seriously affected 
errors estimates service lives the interest methods. 

all times prior the end the service life plant unit, the straight- 
line method results higher depreciation reserve than the interest 
methods. the case long-lived plant and where high interest rate 
used under the interest methods, the difference between reserve com- 
puted the straight-line method and one computed according the 
interest methods will substantial. Conversely, the reserve computed 
according the interest methods will approach the straight-line reserve 
the interest rate approaches zero. 

The straight-line method results higher charges consumers for depreci- 
ation and return (assuming the same rate return) when property 
young, and lower charges when property old, than does either the 
interest methods. can argued that under the “present-worth” 
theory the real cost consumers the long run will the same, 
regardless the depreciation method, provided extra profit accrued 
the utility because differential between the sinking-fund interest 
rate and the rate return. The total dollar charges customers the 
long run, however, will greater under the sinking-fund compound- 
interest methods than under the straight-line method. 


25. The compound-interest method, which also called the present-worth 


method, and the sinking-fund method, which modification the 
former, yield the same results the aggregate, where the depreciation 
reserve requirement deducted under the compound-interest method 
and where the interest rates under both methods are the same the rate 


return. 
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26. The straight-line method generally recommended for public utility 
accounting and financial purposes for the computation both 
depreciation expenses and accrued depreciation for purposes rate 
making. 

27. Depreciation rates may applied individual units depreciable 
property groups units. The group plan has many advantages 
and should used wherever property can rationally grouped. 

28. The determination depreciation rates according any reasonable pro- 
cedure involves the use judgment the application methods. 
Statistical data drawn from property experience and applied with judg- 
ment based knowledge the property should given much greater 
weight than general judgments based principally observation the 
property. 


CHARGES AND ADJUSTMENTS RESERVES 


29. Current charges for depreciation expense should based upon the best 
possible estimates the amount properly applicable the period 
covered income statement, without modification for excessive 
deficient charges the past. 

30. The depreciative reserve should not readjusted gradually through 
modification the annual depreciation rates when the difference between 
the book reserve and the proper reserve substantial. When the differ- 
ence not substantial satisfactory spread the remaining net cost 
the properties over their remaining lives. 

31. principle any necessary correction depreciation reserves should 
made through surplus special section the income account. How- 
ever, the adjustment inadequate depreciation reserves, while sound 
theory, presents many practical difficulties. Where such deficiences are 
serious, desirable make every effort adjust them, although 
recognized that the application uniform rule without regard what 
equitable and feasible under the circumstances individual cases 
might cause injury security holders out proportion the long-range 
benefits. Therefore, concluded that the objective correcting 
inadequate reserves should approached with appropriate consideration 
the practical effects alternative courses action. 


DEPRECIATION INCOME TAXATION 


32. Because certain special problems and provisions the federal income tax 
laws and regulations, absolute consistency between the treatment 
utility depreciation for tax, corporate, and regulatory purposes 
may not now feasible, but such consistency should sought the 
greatest practicable extent. 

33. Inconsistencies between the treatment public utility depreciation for 
income tax, corporate, and regulatory purposes have sometimes been 
unjustifiably great and often have been inadequately disclosed with 
attendant effects upon the tax liabilities stockholders. Analysis 
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substantial differences the treatment depreciation for these re- 
spective purposes necessary, both for the protection the investor and 
for fully effective regulation. 

34. Since the depreciation reserve builds more rapidly under the straight-line 
method than does under the interest methods, the use the straight- 
line method tends reduce the investors’ risk and therefore lower the 

cost capital. 


DEPRECIATION AND THE INVESTOR 


Failure make adequate provision for depreciation results the impair- 
ment capital, unless appropriate retention surplus made. 
Depreciation policies and practices should, therefore, taken into 
account the consideration questions public utility finance. 

36. Although the immediate interests equity- and debt-holders may some- 
times conflict, and although the short-run and long-run interests 
investors class may sometimes differ, adequate provision for de- 

preciation essential the protection the interest investors, taken 

whole and from long-range point view. 


DEPRECIATION RATE MAKING 


The same factors which cause annual depreciation expense also cause 
accrued depreciation. Accordingly, the same principles should 
applied determining annual depreciation expense and determining 
accrued depreciation. 

38. So-called observed depreciation not accrued existing depreciation. 
Depreciation not visually observable; only the physical effects some 
the forces which operate cause retirement can observed. 

39. properly computed depreciation reserve the best measure accrued 
existing depreciation, since such reserve reflects that part the cost 
the property service which relates the exhausted expired economic 
service life. 

40. fixing public utility rates adequate depreciation expense should 
allowed according the service life basis and, principle, the reserve 
requirement determined the same basis should deducted de- 
termining the rate base. There are cases, however, where equity re- 
quires the deduction the book depreciation reserve, even though 
differs materially from the reserve requirement. The determination 
when this condition exists involves, among other things, consideration 
the history regulatory requirements the jurisdiction and the experi- 
ence and practices the company. other words, individual situations 
should dealt with the commission having jurisdiction accordance 
with the facts each case. 

41. The use the straight-line method determining both depreciation ex- 
pense and accrued depreciation generally recommended for rate- 

making purposes. 
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42. The sinking-fund method may sometimes used rate making when 
impracticable determine the accrued depreciation. When the sinking- 
fund method used rate proceedings, the interest rate employed 
should the same the rate return which applied the unde- 
preciated rate base. 


RESPONSIBILITY MANAGEMENT AND PuBLIC ACCOUNTANTS 


43. Management and public accountants have responsibility indicate 

financial statements the proper depreciation expense and the reserve 
requirement, according their best information after appropriate study, 
those instances where the book figures are inadequate excessive 
material amount. 
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AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 1852 


DISCUSSIONS 


FLOOD FORMULAS BASED DRAINAGE 
BASIN CHARACTERISTICS 


Discussion 


CLARENCE JARVIS 


Am. Soc. the many hours required 
for the writer grope through the maze interrelationships disclosed this 
paper may taken index the authors’ labor exploring, analyzing, 
and then formulating the wealth material the paper, the profession has 
received the harvest months not years earnest effort. Recognizing the 
practical value the authors’ masterly approach and the comparatively unique 
not new relationships presented, the writer somewhat reluctantly undertakes 
examine background foundational features, trace the derivations, and 
analyze the adopted factors throughout the practicable range, -thus intro- 
ducing both alternate and divergent views regarding major well minor 
issues—all intended constructive suggestions throughout. 

Translation the quantities listed Cols. and Table into 
terms percentage ratings according the Modified Myers formula‘ (Eq. 
readily discloses the following relationships: 


Quantities listed Table proposed supplement Table form 
geometric series with the following characteristics— 


(a) Percentage ratings for major peaks (Col. 10) are usually from 1.5 
2.5 times those for minor peaks (Col. 9), averaging nearly 2.0. 

(b) Corresponding ratios between percentages listed for rare peaks 
(Col. 11) and for major peaks (Col. 10) are usually from 1.5 
2.0, averaging nearly 1.7. 

(c) Between maximum peaks (Col. 12) and rare peaks (Col. 11), the 
ratios are usually from 3.2 4.2, averaging nearly 3.6. 


paper Kinnison, Am. Soc. E., and Colby, Esq., was published 
March, 1944, Proceedings. 


Senior Civ. Engr., Intelligence Div., Office, Chf. Engrs., Washington, 
Received the Secretary May 1944. 
Flow Jarvis, Transactions, Am. Soc. E., Vol. (1926), pp. 
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JARVIS FLOOD FORMULAS Discussions 


(Maximum Observed Flood Ratings Are Underlined) 


single setting slide rule per item). 


Maximum known peak discharges correspond values Col. Table 
Peak values computed from flood formulas (see Cols. 12, Table 1); percentage 


which peak cubic feet per second; drainage area square miles (solved from 


107 4.2 10.6 15.7 18.0 63.1 
405 1.8 1.9 7.0 24.3 
124 1.7 2.3 3.3 5.5 21.5 
23.3 1.9 3.5 6.2 26.1 
183 1.6 2.3 3.5 5.5 21.1 
35.2 1.8 3.8 6.7 27.0 
211 2.1 4.0 6.9 26.6 
42.2 1.6 3.3 5.4 22.2 
31.3 4.5 5.1 8.9 34.8 
139 2.4 6.4 10.8 37.6 
417 3.1 8.1 14.2 47.6 
224 8.0 9.4 17.0 28.3 95.5 
308 16.6 26.3 29.8 43.3 
71.2 4.9 6.2 118 19.6 
420 8.8 9.7 19.5 31.7 104.0 
37.6 5.1 10.0 11.7 20.1 73.3 
36.4 4.2 9.9 10.0 17.7 65.3 
83.8 2.4 5.8 11.7 37.8 
186 2.7 6.6 11.6 40.9 
370 3.7 9.9 15.1 16.6 55.0 
19.0 1.8 3.9 6.0 6.9 28.9 
18.7 1.6 3.8 6.2 6.9 25.0 
49.9 3.0 7.3 7.6 13.1 48.2 
12.2 2.0 4.4 5.2 8.7 34.6 
362 21.1 29.6 31.6 51.0 171.0 
178 17.5 26.7 42.2 46.0 159.0 
88.4 10.5 18.6 32.2 114.0 
89.4 7.1 14.6 23.8 84.2 
53.3 5.4 12.1 20.4 73.5 
100 2.6 6.5 11.3 14.0 40.8 
201 3.1 8.6 13.7 16.0 46.3 
702 5.6 13.7 17.0 22.8 72.8 
188 2.9 5.5 7.6 12.6 43.4 
149 2.8 6.9 7.7 12.6 
162 10.0 18.3 29.8 30.5 104.4 
52.6 9.1 16.0 25.9 27.4 93.5 
93.7 10.2 19.2 22.6 32.2 115.8 
98.4 2.8 6.8 11.3 41.5 
75.3 6.6 14.3 21.3 25.2 113.8 
4.1 1.9 3.2 3.5 6.1 27.0 
74.5 2.4 6.0 6.4 10.6 38.7 
57.1 5.1 8.5 11.9 20.7 74.3 
280 5.0 6.9 12.4 19.8 65.1 
31.4 9.1 9.7 16.8 30.3 113.0 
39.1 11.5 16.0 19.0 34.0 125.4 
See Fig. 
Years used the study flood peaks (see Col. Table 1). 
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(d) notable degree consistency pervades most the series, mani- 
fested tendency adhere either average values 
above-average below-average values throughout. 

(e) Random selection items for checking computations, mainly 
log-log slide rule, resulted satisfactory checks, within 
less, differences being partly attributable the mechanics the 
procedures and cumulative errors. 

(f) Considering all the factors concerned with flood frequencies and 
magnitudes, seems advisable re-examine all items Table 
with ratings according the Modified Myers formula for ‘‘maxi- 
mum” peaks (Col. 12) greater than 90%. Very few, any, 
these stations seem deserve such high ratings according the 
writer’s interpretation all available data pertinent the 
locality involved. 


Factors involved Eqs. have been platted Fig. semiloga- 
rithmic paper throughout the practicable range for most rivers. depict 
all the trends single drawing, was necessary convert the three fac- 
tors dealing with pondage percentage values, multiplying 100. 
Furthermore, was required use another factor for the authors’ expression, 
take best advantage the chart space. The other two factors were sug- 
gested the writer more nearly expressive the probable influence 
pondage. was also necessary divide the quantities the factor involving 
keep the right-hand curve within the limits the drawing. 
The outstanding characteristics curves relating length runoff travel, 
drainage area, and altitude above the outlet are increased curvature and di- 
vergence above 100. 

The general form Eqs. with fractional power drainage area 
the numerator and with lesser power runoff travel distance the denomi- 
nator seems recall the derivation the formula derived 
Pettis, Am. Soc. E., directly from the Modified Myers Eq. 
Through the introduction “compactness factor” Colonel Pettis obtained 
the expression, 


which represents the rainfall depth inches for the critical adopted 
period and the average width the drainage basin miles, 


review and summary Colonel Pettis’ paper readily available.) 

Within the ordinary range observed drainage areas, where length 
some multiple the width, possible convert the authors’ fractional 
proaching that used Colonel Pettis, with some advantage the way 
simplification, but with appropriate changes coefficients. 


New Theory River Flood Flow,” Pettis, 1927 (privately published). 


the United States; Magnitude and Frequency,” Water-Supply No. 771, 
Geological Survey, 1936, pp. 36-37. 
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With 681 station-years record included for forty-eight items Table 
giving average 14.2 years per station and maximum years per 
station, and with only four the observed maxima exceeding the quantities 
listed Col. (the amount such excess being 25% less), 
difficult reconcile the idea adding another 300% more attain the 
“maximum determined the authors. Addition 150% might 
well suffice. 

Critically examining the various factors Eqs. graphically outlined 
Fig. (with such alterations seemed essential keep them within 
limits the drawing), the following comments are submitted— 


(a) The values vary between 3.09 and 15.5 for miles and 
miles, respectively, 77.5 for 500 miles. Thus, the L-factor 
appears the denominator each formula, its influence would 
such decrease flood peaks 1/25 for such range 
variation. For the alternate Eq. 6b, with instead, the 
total range between and 500 would instead 25—a more 
readily defensible factor according extensive records and 
observations. 

(b) The s-factor for minor floods, with added constant, amounts 
126.27 for 100 altitude, 137 for 200 ft, 237 for 500 ft, 
694 for 1,000 altitude above the outlet. Thus, because 
the influence altitude Eq. 6a, the flood peak multiplied 
factor consequence increased median altitude from 
200 1,000 ft. All three the other s-factors and associated 
constants run off the drawing before attaining altitude 
600 ft. Their trends magnify the influence altitude drain- 
age area above the outlet greater degree 1,000 than has 
been noted for minor floods. 

(c) The pondage factor, used the authors shows but slight 
variation between and 5%, attaining value only 1.067 
1.169 50%. other words, the regulation afforded 
ponds covering fully half the drainage area would decrease the 
peak 1/1.169 85.6% its magnitude with pondage. 
the fractional power were 0.4 instead 0.04, the corresponding 
values the a-factor would become 1.90, 3.90, and 1/3.90 
25.7%. Thisis more nearly accord with the writer’s experience, 
although 0.2 might serve somewhat better, with values 1.38, 
2.19, and 1/2.19 45.8%. 


The proper scope the authors’ equations and formulas this paper lies 
mainly not entirely below 1,000 for both and s—drainage area and median 
altitude above the outlet, respectively. Similarly, the width formula derived 
Colonel Pettis was intended serve primarily for drainage areas ranging 
from 1,000 10,000 miles; occasionally served fairly well outside those 
limits, but more often wandered beyond control. 

For many years the writer has striven for simplification expression within 
the hydrologic field. One the prime essentials seemed formula 
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capable comparing the floods either the Missisippi River the Amazon 
River with the peaks torrential floods their headwaters, regardless 
areal extent. Investigations covering several years led measure ap- 
proval and acceptance the hydrologic and engineering profession, demon- 
strated wide use the square root drainage area convenient first 
approximation the influence areal extent for the entire range above to7 
miles. the same conclusion was reached the result later 
investigation conducted committee the Boston Society Civil Engi- 
neers. This report laid some the foundations for later development the 
field unit hydrographs, has been tacitly acknowledged the originators 
that For the usual form drainage area, with length con- 
siderably greater than average width, introduction the expression may 
regarded equivalent adopting the width the governing dimension 
with either coefficient power greater than unity; might construed 
being the same using the dimension length the denominator, with 
either coefficient power somewhat less than These developments 
have been presented considerable detail some the references the 
paper, well other publications deserving 

spite the time-consuming mazes complex computations involved 
the authors’ procedure and formulas, not mention the occasional excep- 
tions divergences opinion freely expressed the writer constructive 
suggestions, the entire investigation, together with its lucid presentation, 
both profitable and enlightening the profession. The concepts basin 
characteristics for not only flood phenomena but also rainfall and runoff rela- 
tions generally are not new. The late Emil Kuichling, Am. Soc. E., the 
late Talbot, Past-President and Hon. Am. Soc. E., Major Myers, 
and many other authors formulas made special provision for accomodating 
the various factors throughout their probable range. The writer’s modifica- 
tion the Myers formula included extension the range both above and 
below that originally suggested Major Myers. The use the Modified 
Myers formula simple and elementary deter its adoption those 
who presumably seek the more abstruse and complicated forms. matter 
how involved any preferred procedures, the results are readily translated into 


terms Myers percentage ratings, demonstrated Table 
Committee Floods, Journal, Boston Soc. Civ. Engrs., September, 1930, pp. 297 and 


Applied Hydraulics,’’ Calvin Davis, Editor, McGraw-Hill Book Co., Inc., New 
York and London, Ist Ed., 1942, Sections pp. 1-132. 


Stream-Flow with Reference Flood-Formulas,” Jarvis, Transactions, Am. 
Geophysical Union, 1937, pp. 409-419. 


Floods the United Jarvis, ibid., 1939, pp. 157-166. 


Annual Yield from Ungaged Drainage-Basins and Effect Land-Use Treatment 
Surface Jarvis, ibid., 1941, pp. 875-893. 


River Discharge from Measurements Jarvis, ibid., 1942, 
pp. 


from the Field Hydrology,”’ Jarvis, ibid., 1943, pp. 651-693. 
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FAILURE THE TACOMA NARROWS BRIDGE 


REPORT THE SPECIAL COMMITTEE 
THE 
BOARD DIRECTION 


Discussion 
HoMER HADLEY 


Homer Assoc. Am. Soc. additional com- 
ment may well made the depth-to-span ratio the late Tacoma Narrows 
Bridge 350) which conjunction with its narrow breadth made struc- 
ture extreme slenderness and flexibility. number places the report 
the Special Committee the Board Direction and its appendixes refer 


this the minds many people seemingly constitutes 


damning indictment the bridge and itself the full and sufficient cause 
the disaster. Therefore, appears desirable direct attention certain 
other structures which not conform this criterion but stand denial it. 
Such criterion, properly limited and restricted, may well accepted 
every one—in its field application. Great the accomplished fact. The 
bridges that have been built and stand and serve their purpose fully and well 
should given the most thoughtful and respectful consideration. Reality 
them, and reality greater than theory—greater than reams mathe- 
matical calculations, greater than laboratory experiments, matter how 
technical and elaborate. Because this is, however, warrant for arbi- 
trarily taking one aspect reality and declaring that the essence the 
whole. 

Bridge Table Report (a) none other than the George Washington 
Bridge New York, Various data concerning are given and its 
depth-to-span ratio recorded footnote corrects the tabulated 
dead weight and the number-of-track values for present single-deck con- 


report was published December, 1943, Proceedings. Discussion this report has 
appeared Proceedings, follows: April, 1944, Francis Witmer. 


Regional Structural Engr., Portland Cement Assn., Seattle, Wash. 
Received the Secretary May 1944. 
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dition” but fails correct the depth-to-span ratio for this condition. Foot- 
note (*), Table states that “depth” Col. stiffening girder.” 
Accepting this the value used “the present single-deck condition” 
the depth-to-span ratio found about 1300, since the stiffening girders 
are more than the top chords the future trusses which will built when 


the lower deck added. the greatest possible combination longi- 


tudinal stiffening members fascia girders, ratio 
still 350, the same that the ill-fated Tacoma Narrows Bridge. 

Undoubtedly minor movements occur the deck the George Washing- 
ton Bridge but far the writer has ever heard these movements are very 
minor, spite the open exposure wind and storm. Why should this 
condition be, fairly low essential for stability? The great weight 
the deck and the four 36-in. cables one explanation offered. Another 
course the fact that there very little projection this deck for hori- 
zontal wind impinge on—little cause for wind-stream disturbances across 
its upper surface; but whatever the cause the stability the fact remains that 
here demonstrated the inadequacy any broad and sweeping generalization 
that would seek establish certain criteria for judging the 
ciency design. 

Another example the Ronosfoss Bridge Norway with center span 
movement this bridge due wind was reported, and carries 
traffic satisfactorily.” 

very remarkable suspension bridge—the design which truly fin 
example engineering skill and judgment its adaptation the site and 
the need for extreme economy—is the Royal Gorge Bridge 
crosses the canyon the Arkansas River, 1,053 above the tracks the 
Denver and Rio Grande Western Railway. Its main span, 880 ft, not much 
less than that the Peace River Bridge the Alcan Highway British 
Columbia, Canada, completed 1943. floor system, which suspended 
from the cables 10-ft. panel points steel suspender rods, consists 
15-in. floorbeams, which turn support seven lines 8-in. I-beams and 
two lines 8-in. channel This construction yields maximum 


the use wind cables, stability obtained. Wind cables are evidently 


Still another limitation that the basic shape and form the deck 
itself and the resistance offered the passage wind around it. There 
should question about the importance and reality this fact any 


Engineering News-Record, December 16, 1943, pp. 74. 
June 1930, 922. 
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one’s mind. The analogous condition the flow water around bridge piers 
and the effect pier shape freely recognized The effect shape 
shown the airplane wing. Time was when the biplane was the dominant 
form, with its upper and lower trussed together struts and wires. 
Where the biplane today? practically obsolete and has been replaced 
the monoplane not because the biplane was structurally inadequate but 
because the monoplane much better suited its purpose. 

the experimental work with the 100-ft-long model the Tacoma Nar- 
rows Bridge jointly sponsored the University Washington, the Wash- 
ington Toll Bridge Authority, and the Public Roads Administration, 
Professor Farquharson found that the scale-replica model the open steelwork 
the deck that failed (that is, 8-ft side girders, floor beams, stringers, wind 
trusses, etc., without floor slabs, sidewalks, diagonal stays midspan) was 
extremely sensitive and flexible winds low velocity (see Fig. 12(a)). 


Heavy Paper 
Wind 
Heavy Paper 


(a) APPROACHING DESTRUCTION (b) STABLE WIND 
30-MILE WIND 


developed more pronounced movements than the prototype and 30-mile 
wind, also laboratory scale, whipped violently that the test had 
terminated. Following this the same identical frame was given paper encase- 
ment and was provided with semicircular nosing both sides, shown 
Fig. When this model was exposed the 30-mile wind, there was 
movement. hung inert. Vigorous undulatory and torsional motion was 
imparted the model and, soon hands were removed, the motion damped 
out and again hung still. all winds 125 miles per the behavior 
was the same with the exception that the 125-mile wind the imparted motion 


continued. The did not change. The moment inertia the deck 


did not change. The cables did not change. What had changed, however, 
was the flow the wind stream around the deck; and that action very 
wonderful change was effected wind loads. Tests along these particular 
lines have not been extensive but sufficient work has been done indicate the 
revolutionary possibilities inherent the situation shape and form can 
substituted for and great weight. 

Report (a) states there doubt that sufficient knowledge and 
experience exist permit the safe design suspension bridge any prac- 
ticable span.” What this “practicable one thing current 
convention must adhered to; another thing can revised and 


Channel Obstructions,” David Yarnell, Technical Bulletin No. 442, 
ovember, 1934 
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simplified. one thing with lavish economy; another thing lavish- 
ness must end. Certainly one can gainsay the fact that long-span suspen- 
sion bridges entail heavy costs. Neither can gainsaid that there are 
numerous places the United States where, due navigation requirements, 
depth water, limited traffic, etc., narrow long-span bridges are badly needed 
but are not “practicable” under existing standards. The Tacoma Narrows 
site was such place. was effort meet the very real and urgent 
necessities the project lowest possible cost that its design was developed. 
Unfortunately the roadway was placed below the top the girders with 
the resultant turmoil and turbulence the wind stream which has revealed 
itself every bridge with I-shaped stiffening members. If, however, with 
necessary readjustments framing, the roadway had been placed the top 
the girders, the section had been closed with plywood the bottom the 
girders, and fairing the form plywood nosing had been added the sides, 
appears very probable that the bridge which failed disastrously would have 
been standing today, giving 100% satisfactory service without the slightest 
change its width its 

Every one wiser after the event. What apparent today was far from 
evident yesterday. The dramatic failure the Tacoma Narrows Bridge has 
profound message for all engineers, but imposes necessity adhere 
rigidly past practices. Its lesson other than this: Whatever you 
must right with Nature woe unto you! 

bit high school Latin absurdly clings the writer’s memory when 
more important things are forgotten—two lines from Ovid concerning Phaé- 
thon, the bold youth who essayed drive the sun chariot his father, Apollo, 
across the sky. The wild horses the chariot broke from his control and 
dashed headlong their course, careening and down, scorching earth and 
heaven. length Jove was compelled strike Phaéthon down with 
thunderbolt and fell earth. The water nymphs gathered his remains, 
buried him, and placed stone his head which they inscribed these lines: 


situs est Phaéthon, currus auriga paterni. 
Quem non tenuit, magnis tamen excidit 


Freely translated they are: 


lies Phaéthon, driver his father’s chariot. 
Although was unsuccessful, failed great undertaking.” 


the first Tacoma Narrows Bridge let remembered that its failure 
was failure great undertaking. 
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Assoc. Am. Soc. clarify this discussion, the 
method framework analysis presented the author summarized 
has follows: 
redundant unknowns are divided into two linear dis- 


| 

placements the joints; and rotations the joints, for which the moments 

hen the joints can substituted the case relatively simple structures. 

2.—Group (a), Item represented the relative variations lengths 


(strains) the sides statically determinate triangulated system which 
and “maps” the actual structure according method well known surveying. 
and 3.—Group Item can calculated relatively simple manner 
function group (a) when the framework the type; that is, when 
the computer progresses along the framework definite direction without 
intersecting hinges (free rotation) fully restrained multiple supports (no 
rotation), member can encountered twice. this special case, moments 
can considered directly, instead rotations. The procedure becomes 
increasingly complicated when multiply-connected frames are considered and 
becomes necessary introduce third group redundant unknowns, which 
must added the moments and which includes the rotations some 
the joints. 
4.—The equilibrium moments the joints obtained without consider- 


ing the equilibrium the longitudinal and transverse forces. These forces 
can defined terms the moments the joints. Then the structure can 
replaced equivalent hinged structure and the conditions stability 
and equilibrium forces corresponding that structure give the equations 
necessary define the unknowns group (a), Item 


paper Yves Nubar, Am. Soc. E., was published April, 1944, Proceedings. 
Structures Engr., Consolidated Vultee Aircraft Corp., Allentown, Pa. 
Received the Secretary May 1944. 
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profitable comparison can made between this paper and the simple 
and direct analysis presented David Am. Soc. E., who 


(1) Divisions the unknowns into two groups the basis both methods. 

(2) For group (a), Item Mr. Wilson uses rotations and Mr. Nubar uses 
strains. The two groups are similar, and only practical considerations can 
indicate the superiority Mr. Wilson limits himself structures 
which axial deformations are negligible, which the “less general case” 
treated Mr. Nubar; but the treatment the more general case suc- 
cessive approximations well known; the conventional method for the 
who examines trusses, beams continuous supports, and frameworks. 

(3) For group (b), Item Mr. Wilson uses rotations moments, and 
applies slope deflection and moment distribution methods. Mr. Nubar uses 
moments wherever possible and applies the method fixed (conjugate) points. 
All these methods are well known and the special conditions each problem 
or, more often, the habits the analyst, dictate the choice between these 
methods. However, for multiply-connected circuits, the moment distribution 
method has the definite advantage being well developed, available with 
numerous examples any textbook. contrary, the fixed points method, 
developed many years ago straightforward manner for open circuits only, 
far from being perfected (as shown the author’s statements) for use 
analyzing multiply-connected circuits. The author indicates, substance, 
that some rotations must introduced auxiliary unknowns and that the 
system equations containing these unknowns must then solved. 
known and condition that must applied any statically indeterminate 
structure. method new and useful essentially the extent which 
presents definite and straightforward manner choosing unknowns and 
solving systems with high degree indeterminacy. 

(4) The determination group (a), Item when group defined 
terms group (a), identical the Nubar and the Wilson methods. 


From the foregoing, appears that the Nubar method follows closely the 
pattern conventional methods. remains decided whether the 
method introduces element newness and simplification over conventional 
procedures. 

indicated, the division unknowns two groups conventional. 
Generally the unknowns group (a) are: The horizontal displacements for 
multi-story frames; the vertical displacements for Vierendeel beams; and the 
axial loads (or strains) for triangulated trusses. 

The author replaces unknowns that result from choice adapted the nature 
each problem, system arbitrary unknowns which are the same all 
cases. This attitude can compared that analyst who would always 


Rigid Frames Superposition,” David Wilson, Proceedings, Am. Soc. 
February, 1944, 129. 

Method Successive Approximation for the Analysis Statically Undetermined 
Calisev, Publications, International Assn. for Bridge Structural Eng., Vol. IV, 1936, 
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use rectangular coordinates, even other sets coordinates (polar, cylindrical, 
spherical, etc.) are more adapted the problem. Such approach utterly 
unsatisfactory; and sometimes the main difficulty the solution problem 
the choice proper system coordinates. the present case, since 
always necessary calculate the angular displacements the straight 
lines connecting the ends each member (called unit deflections the 
author), appears logical, wherever possible, use the latter principal un- 
knowns. For instance, Fig. advantage results from the use 
unknown, against the use either the horizontal displacement the 
common rotation EA, Defining this rotation, Eqs. 36, 37, 
and become: 


= = + 61-2) = 2.777 6 


The simplification obvious and clearly demonstrates, means the ex- 
ample given the author himself, that the method very little practical 
interest for frameworks, even unconventional type. 

For conventional framework, there doubt that the linear displace- 
ments the various stories, defining the rotations glance (even better) 
the rotations each line columns between stories, are the logical unknowns 
used (Fig. 14). axial strains are considered, they are added the 
rotations form the first group unknowns. 


His illogical choice the unknowns for rectangular frames has forced the 
author use, all cases, what called compatibility” which 
relate axial strains angular displacements. These equations are well known 
under other names. are the equivalent the Williot-Mohr con- 
struction. The condition expressed terms moments (Eq. the 
used for the calculation deflections curve beams, the sole dif- 
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ference being that the quantity replaced the 


quantity generally replaced the integrated ex- 


pression: 


Eq. 10a well known, English" well literatures under the 
equivalent form (see Fig. 15): 


Topology purpose the 
section entitled “Topology Frameworks” appears 


appears several times the paper. However, this 
proposition self-evident without any preliminary proof: The each 
have two degrees freedom the plane; the entire framework joint dis- 
placement has degrees freedom, plus degrees the linear dis- 


placements support, minus degrees for the rigid (nondeformable) 
members. Thus: 


52, (in which index and not exponent) may vary from 
the total number members. appears that this entire section 
extraneous the subject the paper. 

Boundary Conditions—Since moments and axial strains are used 
unknowns, and, perforce, they define the stresses and loads the structure 
completely, obvious that the boundary can expressed 
terms such elements, and general Eqs. not clarify the problem for 
practical applications. 

End Moments.—The method used determine succession the moments 
structure the “open” type the well-known method ‘‘conjugate 
points” (English terminology), (French terminology), “fixed 
points” (German terminology). 

The development Eqs. can presented simpler and more 
metric manner: Using the author’s notation, and with denoting any member, 
there exist relations the type 


being the joint the end the member r), relation the same 
holds for the nth member. 
402. 
des Fontviolant, Vol. II, 1927, Paris, 148. 
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This fact obvious from the general relations linearity between moments 
and rotations. The coefficients the nth member can related 
the other coefficients. The equation equilibrium joint is: 


which couple applied the joint. Expressing the rotation any 
member the joint terms the loads and couples, 


Wy = + 6 Ek; (2 Mrs MrR = (55) 


All rotations are equal. Replacing terms means Eq. 
53, the value can expressed follows: 


M 6E k, 4 + 2 ms mr + Des 
Wy C, + 2 C, 4 2 
Substituting Eq. 54: 


E k, 6, 2 alia r D, k, 


When m,, and are all equal zero, there may rotation 
under the action moments This shows that the coefficient 
equal zero (except for completely restrained joint for which and 
2). Thus 


Also 


Eqs. and are absolutely symmetrical and are equivalent Eqs. and 
developed from Eqs. and mere extension the method fixed 
points for continuous beams, adapted “open” structures. that form 
was presented Prof. early 1936. 

For multi-connected circuits, necessary introduce the rotations 
some the joints, which cannot eliminated directly. Very rapidly, the 
method becomes extremely complicated, due the number simultaneous 
unknowns, and loses its value, least its present stage. Moreover, the 
author limits himself simple suggestion the method, without attempting 
develop into practical procedure. For conventional truss, the number 
rotations introduced the integer equal immediately superior 
one third the number the circuits. Fig. 16(a), the “opened” joints 
are represented hinges. For conventional Vierendeel structure (Figs. 
16(6) and 16(c)), the corresponding ratio one half. The last result also 
applies one-span, multi-story portal, and two-story, multi-span 
portal (Fig. 16(d)). For multi-span, multi-story portal, the number 


for Bridge Structural Eng., Vol. IV, 1936, Zurich, 509. 
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approximately two thirds one half the order multiplicity the 
portal. For instance, the portal represented Fig. 16(e) 42d order. 
can reduced the 6th order shown, and then the “open” type 
shown Fig. 16(f). Thus total twenty-eight rotations are needed for 
framework the 42d order, and four rotations are required for framework 
the 6th order. obvious that such cases only iteration methods are 
practical interest, and these methods may advantageous use all the 
rotations joints symmetric manner, avoiding the separation any 
that may considered fundamental unknowns order calculate the 
remainder terms the segregated unknowns. 


(a) 


Fie. 


When the author states (see heading, End Moments”) 
that proposition will established that the end moments the members 
rigid frame can found linear functions the 5,_, elements and 
the linear movements the supports,” simply states that solve system 
linear equations the following procedure can followed: 

Solve the equations terms the given constants, some the unknowns 
being considered fundamental unknowns. The remaining unknowns will 
expressed linear functions the fundamental unknowns and the con- 
stants. Thus, the author’s statement, presented new proposition, does 
not reflect any new light the problem. Moreover, the proposition in- 
complete, since the end moments are functions, not only the 5-elements, but 
also the rotations that must introduced multi-connected circuits. 
However, for structures, that procedure, which could designated 
Ritter’s method, great interest because: (1) The number fundamental 
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variables small; and (2) the secondary variables can expressed 
relatively simple manner terms the fundamental variables. 

similar methods were available for multi-connected circuits, great 
step would made toward the general direct analysis structures that 
type. appears that such methods not yet exist. 

Determination the Axial Stresses and Linear Movements the Supports.— 
indicated Mr. the conventional procedure consider appro- 
priate free bodies and determine the transverse forces the members. 
The conditions stability and equilibrium the totally 
hinged structure, the conditions equilibrium 
several free bodies that are statically independent 
one another, define all the conditions required find 
the fundamental unknowns. For instance, the sim- 
ple case the frame Fig. 17, the hinged structure 
unstable. Thus, the equilibrium two free bodies 
must examined: The condition equilibrium body 
defines the load acting body Then the equilibrium body can 
expressed, defining the load acting the vertical member. 

order solve the same problem, the author introduces fictitious diag- 
onal and stipulates that the load that diagonal must equal zero, which 
appears useless complication. Actually, the author, Fig. 13, does not 
use that diagonal, and his method, somewhat artificial, does not give the 
required equation directly. possible find the result directly. For 


example, projecting the forces shown Fig. 18, order eliminate the 
loads along and 
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which strictly equivalent Eq. after the subsequent elimination 
from the latter equation. 
Conclusions.—The following conclusions deserve emphasis: 


1.—The author does not present substantially new method analysis 
redundant structures; merely gives somewhat different presentation 
existing methods. The presentation shows tendency generalize the field 
application some methods. the case presented, this approach 
erroneous, since absolutely general methods are available (least work). The 
trend future research must toward the development special methods 
well adapted each category problems, since absolutely general methods 
are impracticable highly redundant structures, because long numerical 
computations. 

2.—The author introduces conditions that create undue complications 
the simple and well-known processes, such his introduction 
conjugate frames and his treatment the subject Frame- 
also presents well-known concepts with changed and special 
terminology, such compatibility equations, boundary conditions, the incom- 
plete proposition relative the determination end moments linear 
functions the 6,_,-elements, and the fundamental relation the fixed- 
point method. 

3.—Precisely the point where the method seems suggest new develop- 
ments (in the introduction minimum additional rotations), the author 
merely states that these quantities must introduced and then eliminated. 

incorrect interpretation given the number the principal 
unknowns. Actually, the unknowns defining the joint displacements must 
added the unknowns defining the joint rotations multi-connected frame- 
works. 

5.—The author has repeated summary two places his paper—first 
the “Introduction” and again the “Application the General Case.” 
The second summary much more intelligible than the first but, when com- 
pared with the outline given Mr. appears that both summaries are 
practically identical step step, far principles are concerned. 
practical standpoint, however, there doubt that Mr. Nubar’s paper will 
prove much less valuable the engineer who desirous finding simple 
and rapidly applicable method. 


Correction for Transactions: April, 1944, Proceedings, page 447, line 28, 
change “‘they will will page 447, line 38, change 
and will also be” will also page 447, line 46, change 
“all them” “‘all coefficients and, page 449, line 12, change 
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Discussion 


ALLEN CARY AND ADE JASKAR 


twelve cases discussed this paper 
constitute important problems continually confronting the highway engineer 
wherever road construction and road maintenance are undertaken. The writer 
was particularly interested the references geophysical methods test 
applied their solution. The possible value such exploration methods 
the geologist and engineer alike should emphasized. 

For some years the Division Physical Research and Tests the Public 
Roads Administration has been engaged research the use geophysical 
methods test highway engineering. Conclusions resulting from these 
studies and the work several the state highway departments, well 
other agencies the United States and elsewhere, indicate that the methods 
test mentioned the author—the earth-resistivity and the refraction 
seismic tests—may have considerable economic value the solution many 
the problems listed Mr. Specifically, Problem (1) and Prob- 
lems (5) (10), inclusive, appear susceptible solution partial 
solution geophysical tests. The following paragraphs are submitted 
amplification the information concerning geophysical tests presented 
the author the belief that such methods subsurface exploration may prove 
increasingly useful the geologist further demonstrating the useful- 
ness geology highway engineering. 

Problem (1).—The location and delineation deposits road-surfacing 
materials earth-resistivity tests have been described various investi- 
gators during the ten years since 1934 letter the 


paper Marshall Huntting was published December, 1943, Proceedings. Dis- 
cussion this paper has appeared Proceedings, follows: February, 1944, Crosby; March, 
1944, Carl Brown; April, 1944, Kellogg, Berlen Moneymaker, Bean, Bleck, 
Philip Keene, and Jacob Feld; and May, 1944, Marshall, Woods, and 

Belcher. 
Associate Civ. Engr., Public Roads Administration, Federal Works Agency, Washington, 
Received the Secretary April 25, 1944. 
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writer, Henry Porter, Am. Soc. E., reported rather favorably such 
tests the State Highway Department Texas. Investigations conducted 
the Public Roads Administration have further established the usefulness 
this test such work. 

Problem (5).—Both the refraction seismic test and the earth-resistivity test 
have been found applicable the classification excavation materials 
(53). The resistivity test has been applied rather successfully the study 
highway grading projects the Missouri State Highway Department, the 
test data obtained being considerable aid predicting the character_of the 
materials encountered (54). 

Problem (6).—Although little work has been done, yet, determining 
the usefulness geophysical tests connection with landslide areas, would 
seem highly probable that such tests could appreciable value disclosing 
potential slide areas. The refraction seismic test has been used preliminary 
investigations determine the presence impervious shales beneath talus 


slopes which constitute potential sliding surfaces upon which the overlying 


porous material may move (55). The earth-resistivity test should prac- 
tical value locating impervious strata shale clay beneath less consoli- 
dated and more porous soils. Where such conditions exist, road structures 


may benched into the impervious layer instead being built entirely upon 


the porous materials subject slipping along the lubricated surface the 
underlying impervious stratum. the case large fill structures highway 
construction routine resistivity tests made over the fill immediately after 
construction and intervals thereafter may prove value disclosing varia- 
tions moisture content. this way abnormal increases moisture may 
detected within the fill where increasing hydraulic pressures may cause 
slip-out requiring costly reconstruction work and causing considerable incon- 
venience established flow traffic. data are hand substantiate 
such conjecture, however. 

Problem (7).—Both methods test may the study bridge 
foundation materials, the seismic test being more positive its identification 
solid rock formations areas where either solid rock gravel formations 
may expected. The low cost involved and the rapidity with which these 
tests can executed make them particular value preliminary surveys, 
often eliminating the need for the time-consuming and costly methods ordi- 
narily used direct explorations the substrata. conditions warrant, 
the direct methods exploration may used the final location selected. 

Problem (8).—The earth-resistivity traverse constant-depth test should 
aid the geologist locating buried river channels throughout 
area, particularly where sand and gravel are concentrated within such channels 
and where rapid erosion might conceivably occur during flood periods. The 
determination the presence, depth, and lateral extent such deposits might 
well influence decisions the desirability particular bridge road 
location. 

Problem (9).—Earth-resistivity tests have been used successfully 
nection with explorations involving tunnel site hydroelectric develop- 
ment Canada (56)(57). The use this test explorations highway 
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tunnel sites likely largely influenced the depth and character the 
overburden. However, the presence faulting and zones shattered and 
highly altered weathered rock where water may prove troublesome tunnel 
construction should rather easily determined the resistivity test. The 
refraction seismic test has been used several instances the Public Roads 
Administration surveys made connection with the location proposed 
highway tunnels. The data obtained from these surveys have not been pub- 
lished but are contained the files. such studies the broken and weathered 
rocks have seismic wave velocities considerably different from the characteristic 
velocity the solid rock sought for the portal construction. Predictions 
the need for lining for particular tunnel may made, based upon 
intermediate seismic wave velocities indicating fractured, jointed rock mass 
through which the tunnel pass. 

Problem (10).—There have been numerous instances which the resistivity 
test has been used successfully the search for underground water supplies. 
References (58)(59)(60) are examples such studies. times the water 
may located more less directly reason low resistivity resulting 
from various salts carried solution, but more often located through its 
association with some stratigraphic structural formation such clean sand 
and gravel rock formation conducive the entrapment and storage 
water. 


help both engineers and geologists recognize more frequently those 
situations which geology does enter into highway engineering.” How ex- 
tended complete the discussion should open question—it could range 
from the extreme minimum single statement obvious fact maximum 
volumes which all known possible examples were cited. Whether 
not the length and style this paper sufficient for the purpose, therefore, 
largely matter personal opinion and hardly subject technical discussion. 

The author has chosen subdivide the topic into twelve “problems” 
might more appropriate term). inference these 
comprise the entire field cooperative endeavor. Actually there are numer- 
ous other classes problems which are equally more important—such 
(a) selecting the most desirable two more alternative routes for highway, 
either locally regionally, from the viewpoint geologic conditions, (b) locat- 
ing alining highway avoid undesirable local conditions, (c) peat-bog 
problems, and (d) grade-fill problems—to cite from experience alone. 

Mr. Huntting has cited incidents exemplify these “problems.” Natur- 
ally these examples are those with which personally familiar which have 
impressed Any one working the dual field geology and highway 
engineering could supplement these other incidents, perhaps even more 
striking typical. The number and variety such incidents would de- 
termined geographic and temporal range and variety experience, but 
greater lesser degree Mr. Huntting’s examples could duplicated anywhere, 


Geologist and Head, Areal and Eng. Geology Div., State Geological Survey, Urbana, 
Received the Secretary April 28, 1944. 
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they may considered sufficient illustrate the classes engineering 
problems discussed. 

Within the limits the paper, Mr. Huntting has done probably well 
any one could. Because the brevity the paper, many generalized state- 
ments have been made which, like all generalities, may not universally true 
and which exceptions limited application could made. Excluding such 
exceptions, experience Illinois would corroborate all the statements ap- 
plicable within the state, general way, except one: has been 
observed that silt, not clay, the major constituent subgrade material most 
affected frost heave, and the writer does not subscribe the opinion (see 
heading, that clay acts the medium through 
which water rises from the water table capillarity until reaches the frost 
line, where freezes and causes irregular upheaval the road surface.” 
The presence some clay the silt may add the effective capillarity and 
the ice-layer expansion which the principal cause frost heave, but 
material that has large percentage clay seems favor neither capillarity 
nor ice-layer expansion. Theoretically, because the smaller grain size 
clay, capillary water should rise higher clay than silt, but other factors 
apparently operate offset the grain size. 


ALLEN Assoc. Am. Soc. E., anp Jun. 
Am. Soc. contributes the solution airport engineering 
problems much those highway engineering. The author has pre- 
sented twelve possible problems which geologist’s services might profit- 
able highway design. Most these have been major importance 
airport design for several years. The urgency the military airport program, 
the heavy wheel loadings, and the tremendous amounts construction ma- 
terials involved have forced the designing engineers obtain large quantity 
accurate data minimum length time. many, but not enough cases, 
the services full-time geologists have been requisitioned. 

the preparation profiles, procurement truly representative 
soil samples, exploration for fill materials and aggregates, and providing 
drainage data, the geologist has been like the man who charged one dollar 
for fixing the typewriter few moments. Objection elicited the statement 
that charged only one cent for the fixing, but ninety-nine cents for knowing 
how. geologist has been able minimize extensive explorations 
involving test pits, auger borings, etc., his professional knowledge probable 
extents uniformity sediments (both horizontally and vertically) and his 
understanding ground-water behavior. 

addition advising the design the runways, the geologist has been 
inestimable aid water supply and building foundation problems. 

The author has ably discussed the various applications geology 
highway engineering that not necessary enumerate further examples. 
However, the relation between the engineer and the geologist does 
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elaboration. far the best single publication available this subject was 
written Legget 1939 (12). 

satisfying note that the constructive service the geologist 
more appreciated now than any time the past and bids for further recogni- 
tion the future. review civil engineering the past, and its status 
the present, demonstrates that the future trend inevitable. the earliest 
years modern civil engineering, when most projects were small moderate 
scope, the engineer was able cope with his problems being his own geologist 
considering large factor safety (more correctly, factor ignorance) 
his work. When modern society began demand larger and more elaborate 
structures, became necessary for the civil engineer himself specialize. 
Today civil engineering graduate may have such subtitle structural, 
sanitary, railroad, highway, irrigation, commercial, municipal, hydraulic, 
materials engineer. These special fields may even subdivided further. 
Many years general experience are required before man can qualify for 
the title civil engineer. 

Upon resuming responsibility for the design major projects, the civil 
engineer busy with general aspects that has neither the time, nor, 
possibly, the training details, necessary obtain the geological data for 
design. must either delegate this responsibility specialists gamble. 
The latter course has been taken too frequently demonstrated failures 
the past (12). Too often the principal use consulting geolcgists 
determining the causes failure important engineering structures. 

stated previously, the constructive service the geologist becoming 
more appreciated, and bids for higher recognition the future. This prophecy 
supported not only the results mistakes that have been made the 
past, but also the fact that the most favorable sites for engineering projects are 
gradually being The ingenuity the civil engineer being taxed 
the utmost the design large concrete and earth-fill dams, long-span 
bridges, large airports (often built the poorest foundations), intricate high- 
way systems, elaborate water-supply systems, and many other projects. Thus, 
the geologist needed more than ever. 

Where find the experienced engineering geologist, who can fulfil such 
responsibilities, and devote full-time service the engineer, important 
consideration, but one beyond the scope this discussion. 


(12) “Geology and Engineering,” Legget, McGraw-Hill Book Co., 
Inc., New York and London, 

(50) Practical Applications Resistivity Measurements Highway 
Problems,” Kurtenacker, Transactions, A.I.M.E., Vol. 110, 
1934, pp. 49-59. 

(51) “Prospecting for Road Metals Geophysics,” Stanley Wilcox, 
Engineering News-Record, February 21, 1935, Vol. 114, 271. 

“Sand and Gravel Prospecting the Earth Resistivity Method,” 
Stanley Wilcox, Geophysics, Vol. No. January, 1944, pp. 36-46. 
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(53) Exploration Earth Resistivity and Seismic Methods,” 
Shepard, Public Roads, Vol. 16, No. June, 1935, pp. 57-67. 

(54) “The Earth Resistivity Method Applied the Prediction Materials 
sented the Twenty-Fifth Mississippi Valley Conference State 
Highway Depts., Chicago, February, 1933 (unpublished). 

(55) “Seismic Investigations Made Connection with the Relocation the 
Knife Edge Section the North and South Highway, Mesa Verde 

National Park, reported Williams, Public Roads 
Administration, Regional Office, San Francisco, Calif., November, 1938 
(unpublished). 

(56) and Geophysics Engineering,” Dolmage, The Canadian 
Engineer (Toronto), Vol. 67, 1931, 19. 

(57) Study Predicts Rock Conditions Tunnel Site,” 
Carpenter and Leonardon, Engineering News-Record, Septem- 
ber 1930, Vol. 105, 364. 

(58) “The Application Electrical Resistivity Prospecting Ground Water 
Problems,” Tattam, Colorado School Mines Quarterly, 
Vol. 32, No. January, 1937, pp. 117-138. 

(59) Prospecting for Underground Waters Desert Areas,” 
Lee, Information Circular No. 6899, Bureau Mines, 
August, 1936. 

(60) Prospecting for Water,” Loehnberg and Lowenstein, 
Mining Magazine (London), Vol. 55, No. 1936, pp. 143-152. 
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DISCUSSIONS 


RAINBOW ARCH BRIDGE OVER 
NIAGARA GORGE 


SYMPOSIUM 


Discussion 


GEORGE STERN 


the literature practical design arches long span contained this 
Symposium. The design arch bridge 950-ft span represents struc- 
tural problem the first magnitude. The use high-strength structural 
steel with high allowable unit stresses leads structural members small cross 
section and comparatively small rigidity the structure. 

arch rib long span considerable deformation will occur and design 
the usual so-called theory” will involve appreciable errors, was 
correctly stated the authors the first Symposium paper. Many these 
potential errors will avoided the use the approximate method sug- 
gested Symposium. obtain correct results, however, the deformation 
the arch under the design loading must considered setting the 
equations equilibrium—that is, the “deflection theory” must used. 

Although the so-called exact method, outlined the first Symposium paper 
obviously can used determine the final moments trial-and-error 
procedure, direct analysis when not too laborious would preferable. The 
equilibrium equations required for such direct analysis may set 


The bending moment M(z) point two-hinged arch can ex- 
pressed 


Symposium was published October, 1943, Proceedings. Discussion this paper has 
appeared Proceedings, follows: December, 1943, Egidio Genova, and Charles Mackintosh; 
January, 1944, Goodrich, and Kennard Thomson; February, 1944, Spofford, 
Mensch, and Ammann; March, 1944, Louis Balog, William Rapp, and Leon Beskin; April, 
1944, Neil Van Eenam, and Gronquist; and May, 1944, 

Stress Analyst, Aircraft Design Dept., Otis Elevator Co., New York, 

Received the Secretary April 29, 1944. 


Deflection Theory Two-Hinged Arch Ribs,” Stern, thesis No. 516, Dept. Civ. Eng., 


Univ., New York, résumé theoretical studies made the Technical Univ. Vienna, 
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which the true horizontal thrust the arch and the bending 
moment the same point simple beam the same span and cross 
section the rib. The vertical displacement point considered 
positive when downward. The second derivative Eq. 


dx? 


the other hand, the following well-known relationships apply for small 
element curved bar subjected axial force, bending moment, and 
temperature change: 


and 


which the angle formed the end tangents the element ds, 

the change length the neutral axis, and the angular displacement 

the element ds. The remaining notation that used the Symposium. 
From the geometrical relations, 


and 


the following expressions for the vertical and horizontal displacements, and 


dp, are derived: 


and 


differentiating Eq. 102a twice, 


and, combining Eq. 103 with 100, the following expression can obtained: 


the abscissa Eq. 104 may rewritten: 


general, the second term the right side Eq. 105, representing the effect 


n 
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the normal stresses due the axial thrust, has small influence the results. 
The second derivative this equation, substituted into Eq. 99, gives 


Eq. 106 expresses the vertical displacement two-hinged arch rib 
function the horizontal thrust differentiating Eq. 105, the second 
term dependent was neglected. the latter part this discussion, 
however, the final solutions will also given for the case which that term 
not disregarded. 

Another relationship between horizontal thrust and vertical displacement 
can derived from the fact that the relative displacement the two hinged 
supports must zero. Hence, from Eqs. 102, after some transformation, 


substituting Eq. 100a and introducing arbitrary cross-sectional area 
Eq. 107 can transformed mathematically into 


L 
0 


Eqs. 106 and 108 represent two formulas with two unknowns, the hori- 
zontal thrust and the vertical displacement soon the general solu- 
tion the differential Eq. 106 obtained, and can computed. 

Eq. 106 nonhomogeneous differential equation the fourth order; the 
derivation its general solutions far too laborious given detail. 
The procedure followed was first solve the corresponding homogeneous 
differential equation. the case arch rib constant cross section 
(that is, I(x) cos constant), the characteristic values and char- 
acteristic functions Eq. 106 can easily computed since they repre- 
sent simple sine curves. Thus the vertical displacement and the moment can 
expressed the form Fourier series. the other hand, the moment 
inertia the rib varies, the characteristic functions are longer true sine 
curves. obtain these some additional preparatory computa- 
tion must made. The formulas that can derived express these sinu- 
soidal curves are not too elaborate, but their complexity increases the more 
closely the variation the moment inertia considered. 

Once the characteristic values and functions the corresponding 
homogeneous differential equation are computed, the general solution the 
nonhomogeneous differential equation can expressed uniform and con- 
vergent infinite series these previously mentioned functions. This method 
solving series was first used Bleich his deflection theory 
suspension Just the theory arches shows good analogy 
Berechnung verankerter Hana Bleich, Springer, Vienna, 
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that suspension bridges, also the respective deflection theories corre- 
spond rather well similar methods attack are used. the case two 
hinged arch ribs, the aforementioned series converges rapidly that, for all 
practical purposes, more than satisfactory use only few terms the 
series. Thus, the general solution Eq. 106 may written 


which the factor given 


The parameters and are defined the integrals, 


0 

and 

stated previously, Eq. 109 can now substituted Eq. 108. Thus, 
integrating term term and transformation, expression can finally 
derived which contains the horizontal thrust the only unknown value: 

which 
and 

Eqs. 109 and 112, the effect the normal stresses due the thrust 
disregarded. was mentioned before that this effect the vertical displace- 


ment very small. Thus, most cases, Eqs. 109 and 112 will give satis- 
factory results. 

However, the effect the normal stresses not neglected, the general 
solution the differential equation the two-hinged arch problem may 
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which the terms representing the effect the 


normal stresses, average values Au, and are used for area, moment 
inertia, and the angle Using Eq. 114, the following final formula can 
developed: 


2 tan 
(a; + b; H) b; Cc Sis “a 2 | 


The derivation Eqs. 109, 112, 114, and 115 shows that they can applied 
any symmetrical two-hinged arch, irrespective the shape the arch axis, 
the variation cross section the rib, and the loading condition in- 
vestigated. 

The horizontal thrust can now calculated from Eq. 112 115 the 
trial-and-error method. Once its value known, the bending moment, shear, 
and unit stress may easily computed the conventional methods, and the 
vertical displacement can immediately obtained from Eq. 109 
writer informed that Friedrich Bleich, Am. Soc. E., has developed, 
but not yet published, deflection theory two-hinged and hingeless arches. 
uses method similar that outlined herein, but, neglecting the effect 
the normal stresses due the thrust the vertical displacement, equations 
are derived from which can calculated directly.) 

solution for the final moments the preliminary two-hinged Niagara 
arch gave results which agree with those obtained the authors the first 
Symposium paper. Using the method outlined herein, the writer computed 
60% increase the moment the quarter point due the deformation 
the arch rib, compared with the 63% increase obtained the so-called 
exact procedure the first Symposium paper. 

Although few engineers will have the opportunity design arch bridge 
the magnitude the Rainbow Arch Bridge, the record its design and 
construction will prove very valuable future designers similar projects. 
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DISCUSSIONS 


HYDRAULIC MODEL INVESTIGATION 
MITER GATE OPERATION 


Discussion 


EDWARD SOUCEK 


Epwarp Assoc. Am. Soc. adoption the 
test program ably described Mr. study was made the 
methods which had been used previously design miter gate moving 
machines. The usual basis for the hydraulic load differential head the 
leaves acting direction opposing the motion. The magnitude the head 
commonly fixed judgment—in some cases supplemented observations 
existing installations. The reliability comparative studies for use 
design limited principally uncertainty allowing for differences gate 
length, submergence depth, and kinematics the gate motion; sill 
clearances are also believed affect the hydraulic loads but probably 
lesser extent. 

Some designers compute the varying differential head upon the leaves 
consideration the rates change the areas separated the gates and 
the discharge past the gates. This method gives results that are much too 
small because the loads are principally due water-level variations which are 
transient character and local extent—surges induced the gate movement 
—and therefore not related any simple manner the adjacent water-surface 
areas. 

During the search for design data, paper that dealt with loads moving 
miter gates was found (written German). appeared very promising 
because contained many equations. Upon translation, however, all sig- 
nificant equations were found contain that elusive term—the unknown and 
varying differential head—which had evaluated before the equations 
could used design. Owing the lack reliable design data, use 
model tests was, the author states, the only alternative the use ex- 
tremely liberal safety factors for the design considerable number large 
machines. 


paper Maurice Amster, Assoc. Am. Soc. E., was published March, 
Proceedings. 


Senior Civ. Engr., Engr. Office, Wilmington, 
Received the Secretary May 22, 1944. 
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This situation not only gives Mr. Amster’s paper “scarcity which 
insures that will widely used, but also requires that most careful con- 
sideration given the limitations stated the author and any other 
conditions that might affect the applicability the test results the design 
gate moving general. Operation the model and interpretation 
the test data accordance with Froude’s law are believed correct since the 
phenomena associated with the hydraulic loads are predominantly inertial 
character. Elasticity the gates, effect flow through the lock culverts, and 
reduction effective chamber length caused the nearness large ship 
might mentioned factors not represented the model additional those 
listed the author. Their effects are relatively small but should not always 
considered negligible. 

Under “Application Test Results Other Miter Gates,” the author 
outlines method for application the test data the design dissimilar 
prototypes. The writer believes that the problem not simple this 
section the paper implies. gate moving machines are not geometrically 
similar, the gate operating time not sufficient measure the kinematic 
characteristics the gate movement. the model gates were operated 
uniform angular velocity, rather than with slow movement near the limits 
travel, the hydraulic load for the same operating time would increased 
many times. For this reason, would preferable consider only the 
first degrees the movement for the critical opening operation rather 
than the entire opening period some cases. 

However, even dismissing the questions kinematic characteristics (other 
than operating time), leaf angle when fully open, sill clearance, and other 
factors which were not investigated, not clear that the proposed adjust- 
ment for operating time generally valid. the author explains, the ratio 
prototype leaf length model leaf length establishes the scale ratio. Then 
from Fig. torque curve can selected that correct for particular 
prototype operating time (determinable from the scale ratio and the model 
operating time). may, but probably does not, correspond that which 
the prototype will operate. Fig. 6(a) shows the effect varied operating time 
only for gate with leaf length approximately equal the submerged depth. 
The effect operating time upon peak torque indicated the curves shown 
Table does not follow from these data that Fig. 6(a) valid basis 


MINUTES 


Description 
1.75 2.00 2.25 2.50 2.75 
Torque, millions foot-pounds........ 22.3 20.0 16.5 14.5 12.5 
Percentage time 0.0 14.4 28.5 43.0 57.5 
Percentage torque 0.0 10.5 26.0 35.0 44.0 


for adjustment torque curves terms operating time when the leaf 
length much greater than the submerged depth—a frequent case barge 
locks. Pending the presentation further data conclusive analysis, 
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would appear that the relationships suggested should used only when the 
leaf length approximately equals the submerged depth. 

The treatment the forces due density differences covers the ideal case, 
static condition with the plane separating fresh and salt water known 
position. reality, determination the elevation this plane quite 
difficult. may interest point out that series tests the existing 
Panama locks, corroborated small transparent model, indicated that the 
excess closing force produced difference density cannot eliminated 
increasing the elevation the culvert outlets. Consideration this problem 
detail beyond the scope this discussion. may suggested, however, 
that designers miter gate moving machines should question the 
elevated outlets reducing the loads upon the machines. 

the Panama Canal Third Locks Project, the terminal gates which are 
most subject the density effect will less deeply submerged than most 
the other gates the machines the terminal gates will have substantial 
overload capacity. all gates, the water the lock culverts 
produces transient opening force which can used reduce the critical 
opening load, either swinging the gates near the time maximum 
overtravel closing the culvert valves and the favorable water 
level differential. course falling tide eliminates the latter possibility 
the terminal gates. 

The writer’s remarks are intended solely clarify certain design problems 
and emphasize the limitations upon general applicability the test data. 
Many significant conclusions can drawn from the paper addition those 
modestly claimed the author. His efforts would justified the paper 
did more than call attention the large reduction operating forces and 
peak power loads which can attained reducing the operating speed the 
gates and nonsynchronous operation the leaves. Since the gate operating 
time relatively small part the time required perform lockage, these 
methods frequently can adopted with scarcely noticeable adverse effects 
upon over-all lockage efficiency. 
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DISCUSSIONS 


STORAGE AND THE UNIT HYDROGRAPH 


Discussion 
FRANKLIN SNYDER 


FRANKLIN Assoc. Am. Soc. dealing with 
unit hydrograph and flood routing procedures have been quite numerous 
recent years and the author complimented for his efforts clarify the 
inherent relationship between these hydrologic tools. However, the extent 
which that objective has been accomplished modified the introduction 
number confusing and unsupported statements about certain hydrologic 
phenomena. 

the “Introduction,” the author, referring article,’ states that 
example, determining seven-period unitgraph from flood resulting from 
three periods runoff, possible determine thirty-six unit hydrographs 
which, applied the three items runoff, will reproduce the composite hydro- 
graph Theoretically, the unit hydrograph procedure were exact, 
there would one unit hydrograph that would reproduce the composite 
hydrograph exactly. For illustration, one can assume unit hydrograph and 
compute composite hydrograph. There exists, therefore, unit hydrograph 
which will reproduce the composite hydrograph exactly. actual practice 
there unit hydrograph that will reproduce observed composite hydro- 
graphs exactly. However, the composite hydrograph not the result too 
complex storm, numerous methods are available for determining unit hydro- 
graph that will reproduce the observed composite hydrograph with satisfactory 
degree accuracy. example’ cited the author illustrates this theory, 
the most probable distribution-graph values obtained 
for composite storm not reproduce exactly the observed hydrograph. 

Fig. shows the results study the effect storage that equally 
function inflow and outflow wave form with initially very sharp 


paper Clark was published November, 1943, Proceedings. Discussion 
this paper has appeared Proceedings, follows: February, 1944, James Sweet, and Otto Meyer; 
March, 1944, Sherman, Gordon Williams, and Ray Linsley, Jr.; and May, 1944, 
Brater, Crawford, Robert Kennedy, and Victor Cochrane. 

Hydr. Engr., War Dept., Office Chf. Engrs., Washington, 

Received the Secretary May 15, 1944. 

Procedure Applied Flood Control Reservoirs with Reference the Muskingum 
Flood Control Project, Appendix IV,” The Engineer School, Fort Belvoir, Va., 1936. 
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peak. The author has made common mistake” obtaining outflow 
hydrograph different from the inflow hydrograph after assuming that the 
reach storage equally function inflow and outflow 0.5, that 
storage varies the sum inflow plus outflow). 

The fact that the outflow hydrograph should identical with the inflow 
hydrograph except for time shift readily proved from the general storage 
equation and evident from consideration the Muskingum formulas, 
6d, 7a, 7b, and 7c. Another result obtained when certain basic require- 
ments fundamental most routing procedures are (1) The routing 
interval, must small enough define, satisfactorily, variations rate 
flow; and (2) the reach should such length that the time travel 
approximately equal the example under discussion 2-hr value 
time travel, was assumed the author for imaginary reach 
river. The value that was used derive the hydrographs Fig. 
not given, but, according requirement (2), hours. Then, with 
end the period equals inflow the beginning. Requirement (1) was 
violated since value equal hours too large define, adequately, 
the hydrograph Fig. However, will always equal from Eq. 

the same basis the outflow curve labeled 0.5 Fig. error 
and should identical the original curve except for time shift hours. 
The negative values for the curve, 0.5, and for the other outflow curves, 
also, are definitely due the selection units and the limitations the routing 
procedure, and couple these values with physical phenomena either real 
imaginary seems unwise. 

Another erroneous result often obtained when routing hydrograph such 
that for Danville Fig. 12, through reach with value considerably 
greater than reasonable value The slow prolonged rise and rapid fall 
coupled with the assumptions the Muskingum routing procedure result 
peak outflow rate higher than the inflow rate even though value con- 
siderably less than 0.5 assumed. Walter Assoc. Am. 
Soc. E., has suggested adjustment the Muskingum procedure 
compensate for the inability the weighted flow represent storage the 
reach adequately when either inflow outflow changing rapidly. 

justify the basic assumptions followed developing the Muskingum 
routing procedure, the value should made approximately equal the 
value this not done (and seldom possible practical applica- 
indicating the portion the weighted flow that derived from the inflow. 

The attempt prove that shifting hydrograph time and routing 
through small value storage gives hydrograph quite similar that 


Direct Method Flood Routing,” Wisler and Brater, Transactions, Am. Soc. 
Vol. (1942), 1519. 
Flood Routing,” Edward Rutter, Quintin Graves, and Franklin Snyder, ibid., Vol. 104 
278. 
Storage and Unit Studies,” Walter Langbein, Transactions, Am. 
physical Union, Pt. II, 1940, pp. 620-62 
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obtained with the Muskingum routing procedure using value 0.5 has 
resulted failure make clear important objective the basic idea. 
Most flood routing procedures actually provide time advantage for fore- 
casting purposes—that is, the inflow data must extended into the future just 
far desired obtain values outflow. The procedure proposed 
the author—namely, shifting the inflow hydrograph time and routing 
through predetermined amount reservoir-type storage—does give time 
advantage equal the amount the time shift. 

The writer has been accumulating data connection with such routing 
assumption and, would expected from the ease with which routing 
procedures can verified, has had little difficulty obtaining satisfactory 
results when routing observed hydrographs for various streams the eastern 
United States. 

analyzing these records procedure convenient for routing the transposed 
hydrograph through certain amount reservoir-type storage was developed. 
stated the author, where storage, and outflow, are uniquely defined 
the elevation the outlet such reservoir 0.0 the Muskingum 
procedure), the ratio successive increments storage increments 


discharge, has the dimension time and represented this discussion 
the symbol storage acre-feet and cubic feet per second, 
units days and equals (very nearly) 0.5 Starting with volu- 
metric expression the storage equation for some satisfactory routing period, 


the mean inflow, minus the mean equals the change storage. 


For units cubic feet per second for discharge, acre-feet for storage, and 
days for 


With the usual assumption straight-line variation inflow and outflow 
during the routing period, and, since then equals with the 


subscript indicating the value the variable the beginning routing 
period: 


Substituting for 0.5 and solving for gives 
which 
(22) 


This solution the storage equation lends itself readily tabular 
putation and the use variable values plotting the variation 
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which recession coefficient the outflow hydrograph equal 
1 

Fig. shows graphically some the time relationships utilized analyzing 

observed records. Local inflow has been separated from the outflow hydro- 

graph that the latter represents the total effect the transposition, 

and the reservoir-type storage, the hydrograph inflow the upper 


Transposed 
” 
a 
=16 Hrs 


Outflow, Equals 
Routed Through Reservoir 


Discharge 


Time, Days 


end the reach. The total time shift between center mass inflow and 
outflow, equal the sum and that is, 


The value also equal the value the Muskingum routing 
procedure provided the fundamental routing requirements previously dis- 
cussed are not violated. However, the inflow hydrograph, Fig. 
routed using the Muskingum coefficients with hours, 0.40, and 
hours (too long reach), the time lag between center mass inflow 
and outflow, will hours instead hours. quite similar 
hydrograph obtained this procedure, using 10, 11, and 

The value any particular case available from inspection the 

hydrograph and the following relationship useful estimating values 
From Fig. 17, tpr equals plus which the total time lag the 
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crest discharge, the crest lag due the transposition, and the lag 
the crest discharge due reservoir-type storage effect; therefore: 


T; = Tr tor + tye (25) 


The time between center mass the inflow and outflow hydrographs, 
readily obtained transferring the hydrograph shapes with carbon 
paper sheet stiff paper, cutting out the hydrographs, and locating the 
centers area. 

the beginning the writer had hoped that values and could 
assumed constant throughout wide range flow conditions, assumption 
that theoretically incorrect but one that made most routing procedures. 
variation can handled without much difficulty but when both 
and vary that their sum also varies, the problem becomes more difficult. 
might expected, observed hydrographs for many streams could 
duplicated assuming constant values the time and storage factors. 
However, when stream with wide flood plain such the Missouri River 
goes overbank, there variation several hundred per cent the factors. 

Under the heading, ‘‘Valley Storage,” the author states that reservoir 
pool above dam may not reduce flood Most engineers will 
subscribe readily this generality, and its acceptance indicated the fact 
that considerable attention given this problem the design and operation 
reservoirs for hydroelectric projects. The author also states that, 
the pool level the dam held constant gate control, the peak 
outflow may larger than the peak The writer could also agree 
with this statement the words “natural peak” discharge the dam site 
were substituted for the words “peak inflow’’; but will take considerable 
more evidence than the Muskingum routing hypothesis convince him that 
the peak outflow can greater than the peak inflow. matter fact, 
because backwater storage above level reservoir pool, the level the 
pool the dam would have lowered somewhat make the peak outflow 
flood large the peak inflow. The peak inflow some cases larger 
than the natural peak the dam site would have been prior construction 
the project because the elimination some valley storage 
because the change timing the concentration the runoff from the 
areas draining directly into the reservoir. 

The author’s extension storage theory ground-water flow quite 
interesting. The use unit hydrograph for distributing ground-water 
runoff may quite practicable, but the use storage values determined from 
recession curves for deriving the rising part the runoff hydrograph subject 
more criticism the case ground-water runoff than the case surface 
runoff, and the procedure open question even for surface runoff. 

Other than the use constant factor for dividing total runoff into surface 
and ground-water runoff, the author’s procedure probably would satisfactory 
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for many purposes. However, the writer takes exception the inference 


q 


that, current practice treating ground-water discharge 
unsatisfactory. Reference made discussing this and other 
phases runoff phenomena which both direct runoff and ground-water 
runoff for the Schuylkill River above Pottstown, Pa. (drainage area, 1,147 
miles) were computed storm storm for year record using only observed 
values precipitation and temperature. This article presents graphical 
comparison the computed ground-water discharge and the observed hydro- 
The accuracy the ground-water phase the computations proved 
more than sufficient for any need the writer flood forecasting hydrologic 
studies. 

The author’s use the quantity, and the time concentration the 
development unit hydrographs confusing. Under the heading, 
Hydrograph Index Storage,” stated that, the time most 
rapid relative rate decrease discharge, ‘“The falling legs large number 
hydrographs flood discharge for streams major drainage basins from 
New Hampshire North Carolina exhibit ratio decrease discharge 
total discharge which, for value 0.5, indicates value closely 
approximating the time elapsing between termination runoff producing 
rainfall and the point most rapid relative decrease With 
0.0 the values for these same recessions determined Eq. 12) 
would equal one half the value time concentration 
the stated that the time-area concentration curve should 
routed through reservoir storage “‘equal times the outflow when the 
time indicated Eq. 12b” the time most rapid relative rate decrease 
discharge when 0.0; that However, the author states 
(see heading, ‘‘Determination the Instantaneous Hydrograph for Streams 
Without Flow amount storage which the time-area con- 
centration curve should modified will approximate 10% 25% the 
time concentration the example for the Appomattox River 
hours) whereas the average value for the James River Lick Run 
determined from the data Fig. was 0.5 hours, hours). 
With somewhat similar procedure, which about the most common and 
oldest known rational method determining synthetic hydrographs, the writer 
has obtained satisfactory unit hydrographs from time-area curves routed 
through storage values equal the unit hydrograph defined time 
from center mass effective rainfall peak the unit hydrograph. 
Values are almost always greater than one half the time concen- 
tration defined the author. 

Howard Turner, Am. E., and Allen using 
similar procedure, obtained fairly good reproduction observed hydrographs 
for the Connecticut River Sunderland, Mass., routing time-area 


Conception Runoff-Phenomena,” Franklin Snyder, Transactions, Am. Geophysical 
Union, Pt. IV, 1939, pp. 725-738. 

Flood Howard Turner and Allen Burdoin, Journal, Boston 
Civ, Engrs., Vol. XXVIII, No. July, 1941. 


1 | 
7 


June, 1944 SNYDER STORAGE 941 


inflow hydrograph through reservoir storage for which equals hours 

The principal purpose the preceding discussion emphasize the fact 
that there are still two variables the development synthetic unit 
storage factor (expressed percentage the author). The use 
slope factor has been introduced the author aid determining the 
time factor, but the storage factor has not been correlated directly with any 
physical features the drainage basin. 

From the writer’s experience attempting correlate the two variable 
factors with physical not too optimistic about the author’s use 
slope factor determine the time variable except localized areas. The 
river slope adopted the author believed more representative slope 
factor than the average over-all slope usually obtained dividing the total 
fall the length. However, the data presented the paper indicate range 


values which from 0.6 for the Smith River Bassetts 


2.0 for the Appomattox River near Petersburg. value average slope, 
defined the author, 11.6 per mile for length miles was deter- 
mined from available reports and used for the Smith River Bassetts instead 
the value 16.6 given the paper.) Although the writer believes the ex- 
ponent the slope factor should smaller than 0.5, variation the exponent 


the expression would have little effect the relative range 


values which vary from 0.6 2.0 for the Virginia drainage basins dis- 
cussed the paper. 

The suggestion the author that the determination the time and storage 
factors for streams without records should based actual measurements 
should emphasized. time factor such can obtained for any 
river station where discharge principally function stage, without installing 
any equipment, observing the time crest stage and comparing with 
field observations official records the rainfall characteristics the storm. 
storage factor not easily determined but not subject great 
range value the time factor and the time factor provides partial index 
the actual value the storage factor. 


1938, nit-Graphs,” Franklin Snyder, Transactions, Am. Geophysical Union, Pt. 


cal 
the 

ber 

ely 
ing 

ites 

on- 

144 
Ts). 

and 
iter 
phs 
ysical 


AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 1852 


DISCUSSIONS 


SURGES PANAMA CANAL REPRODUCED 
MODEL 


Discussion 


POSEY 


this type are all too infrequent engineering literature. When the under- 
lying theory generally accepted, there tendency adjust the design 
fit the results the model tests and omit subsequent tests the prototype. 
Only when the model theory the subject controversy the laborious 
and difficult job securing prototype data for comparisons likely under- 
taken. Many hydraulic engineers would consider satisfactory model-prototype 
correspondence, this particular problem, almost foregone conclusion 
and prototype tests unnecessary. The writer disagrees with this viewpoint, 
and believes that such comparisons should undertaken whenever practicable. 

Although substantial agreement with respect wave velocities evident 
from inspection Figs. and statistical study comparative water- 
surface elevations might give more definite basis for the authors’ conclusion 
(see heading, Model and that agreement 
considered very Plotting the model and prototype test 
data the form shown Fig. should make easier judge the degree 
correspondence. 

The authors’ comparison the results synchronization and 
synchronization, shown Fig. certainly special interest. seems safe 
assume that these are the extreme cases and that surges resulting from 
operation intermediate periods will not exceed the values obtained for 
synchronized operation. 

paper Edwards and Edward Soucek was published January, 1944, Pro- 


ceedings. Discussion paper has Proceedings, follows: April, 1944, Harold 
Weggel; and May, 1944, Harry Leypoldt. 

Research Engr., Iowa Inst. Hydr. Research; Associate Prof., Hydraulics and Structural 
State Univ. Iowa, Iowa City, Iowa. 

Received the Secretary May 29, 1944. 
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DISCUSSIONS 


MILITARY AIRFIELDS 
SYMPOSIUM 


Discussion 


CASAGRANDE, AND HARRY COTTON 


this Symposium the logical development both construction and design 
problems, together with the solutions permitted under the strict limitations 
imposed the expanding, ever-changing, and accelerating activities war- 
time. Colonel Stratton’s interesting relating subgrade support and 
the various methods stabilization compaction, the addition base 
courses suitably rolled, recalled some the writer’s early construction experi- 
ence. addition the usual provisions for both surface and subdrainage, 
was sometimes necessary insure against both lateral and upward flow 
compressible and plastic materials. One the most satisfactory methods was 
the placement crushed slag stone trenches for compaction under pile- 
hammer impact. The virtual interlocking the angular fragments accom- 
plished this manner seemed make such footings react very similarly 
plain concrete slabs. Unquestionably the inertia such heavy masses among 
the footings accounts part for their remarkable stability under varying 
applications live load and required dead load. 

The writer vividly recalls the announcement the division design 
federal bureau concerning proposed thickened edges pavement strips 
replace either the prevailing uniform thickness, even the heavily crowned 
road pavement, with edges concrete thinned the amount rise the 
crown. less startling was the introduction dowels extend across all 
construction and expansion joints, thus preventing the uneven settlement that 
had proved destructive concrete under heavy metal tires. 


Symposium was published January, 1944, Proceedings. this Symposium 
has appeared Proceedings, follows: February, 1944, James Newman, Jr.; March, 1944, 
Thomas Stanton; April, 1944, Howland, and David Jenkins; and May, 1944, William 
Rudolph, Raymond Irwin, Greulich, Hibbert Hill, Jacob Feld, and Robert Horonjeff. 

Senior Civ. Engr., Intelligence Div., Office, Chf. Engrs., Washington, 

Received the Secretary April 10, 1944. 
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Colonel Stratton’s frequent references traffic distribution, repeated 
stresses, and fatigue pavement materials indicate his familiarity with, and 
his practical application of, the basic principles disclosed the. experiments 
the late Van Ornum, Am. and others. How 
else could any one justify the departures from previous practices described 
under the heading, Paving: Allowable Working Stresses for Con- 
crete The use safety factor 1.75, corresponding 57% 
the 28-day flexural strength the pavement, even the use, for limited 
operations, safety factor only 1.40, representing 71% the stress ordi- 
narily required cause failure with single application loading destruc- 
tion, entirely logical view the increasing strength after days and the 
wide dispersion traffic, well the wing lift accompanying the landing 
planes. Who could have foreseen that the vibratory strains the taxiways 
and the warming-up aprons, rather than the stresses due landing impact, 
would govern both current and future design? 

Summarizing the conclusions that may drawn from the Van Ornum 
series tests, application 50% less the loading that would normally 
cause failure may repeated indefinitely without appreciably decreasing the 
resistance reinforced concrete columns, beams, slabs; but, for 75% 
loading greater, only few such applications should tolerated because 
the definite tendency toward fatigue materials, accompanied measur- 
able reduction ultimate strength. 

The observed variation among individual test specimens em- 
ployed the writer research conducted 1905 and 1906 (38), and also 
among the more numerous samples subjected the more elaborate tests 
the Van Ornum series, would occasionally becloud the results. For example, 
loading 85% would cause failure with single application, and the next 
specimen would withstand similar loading 100 more times. 

the light all available tests date, difficult improve the 
range indicated Colonel Stratton’s treatment allowable working stresses— 
57% the ordinary 28-day strength resistance, for limited operations, 
71%. 

Although there are valid objections the use concrete pavements 
thick in., because the attendant temperature stresses and also 
because the difficulty securing the required compaction and density, there 
appears practicable limit the base-course thickness stone slag, 
duly compacted. the writer’s experience, has been demonstrated re- 
peatedly that base course rubble quarry stone, with interstices mostly 
filled with fragments, would support lightly reinforced pavement 3-in. 
thickness plus whatever penetration was attained keying into the stone 
base, with better results and less cost than would plain concrete double 
treble the pavement thickness, and without base course. 

Mr. Hathaway deals with noteworthy accomplishments broad field 
that seems have awaited well-organized, vigorous, clear-sighted, and ana- 
lytic approach. Some other attempts have fallen far short the goal because 


18% Numerals parentheses, thus: (37), refer corresponding items the Bibliography the end 
the Symposium, and the end discussion this issue. 
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the lack required basic data. However, additional material has been 
supplied recent years through the cooperative efforts various federal 
departmental representatives including, notably, Mr. Hathaway and Merrill 
Bernard and Grover, Members, Am. Soc. E., their administrative 
supervisors and technical staffs, legislative committees, and, finally, the great 
mass enlightened citizenry who favored liberal support hydrometeoro- 
logical related research. addition this financial and scientific sup- 


port, the development fundamental procedures exemplified the 
Symposium certainly required the wholehearted assistance specially qualified 
members the technical staffs various field offices under the Chief Engi- 
neers, well those men immediately under Mr. Hathaway’s direction, 
tacitly implied expressed under 

attempts were made include the such technical 
papers and other publications have contributed vitally the background 
and the fundamental knowledge airfield engineering (including meteoro- 
logical, infiltration, runoff, and related phenomena), the list would too long 


and cumbersome for practical purposes. 

The writer realizes that any constructive suggestions for either amplifica- 

tion clarification concepts may mistaken for criticism, matter how 

helpfully intended. However, passengers maritime vessels other 

transportation units, who fail communicate vital information that may 
beyond the range regular lookouts (an open port too near the water line, 

signs shallows impending collision—anything endangering the common 

interest) thereby partly forfeit their claims safe passage. The writer 

deeply concerned about the field hydrologic and related scientific research, 

le, particularly regarding, 

(a) Thorough utilization existing basic data; 

Expansion procurement programs broaden their scope; 

(c) Recognition patterns among natural phenomena, capable afford- 

ing guidance supplying reasonably close approximations missing 


basic data; and, finally, 
(d) Simplification procedures whereby the divergent factors may 
made yield integrated, logical, and conclusive result, capable 


withstanding repeated practical tests consistently and satis- 
factorily. 


The writer’s first glance Fig. disclosed what appeared break 
the continuity the series between the curves for and 120 min. Obvi- 
ously, 60-min curve (intermediate yet not quite midway between the adjacent 
graphs) should have been included Fig. would appear straight 
line connecting the origin coordinates with the abscissa value 4.5 the 
top margin. the very nature this relation, numerical values ordinates 
and abscissas all intersections the grid for the 60-min curve will equal. 
The fact that the hourly rainfall intensities are displayed along the vertical 
axis was not overlooked. 

additional curve that might well drawn Fig. concerns the 24-hr 
Assuming that the only rainfall during such period occurs within 
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one 
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end 
| 


946 JARVIS MILITARY AIRFIELDS Discussions 


the 240-min (4-hr) interval which the curve represents, the average hourly 
4.5 
line would connect the origin coordinates with the point the top margin 
corresponding abscissa 0.75, representing the extreme limit toward the 
vertical axis. Allowing for the gradual increase precipitation prior the 
beginning the 4-hr period high intensity and for corresponding decrease 
later, well for the associated showers which combine make the 24-hr 
maximum for each designated frequency within the series, the final position 
the graph would connect the origin with the value approximately 1.0 
along the top margin. This result logical because 24-in. rainfall depths have 
been attained rarely, ever, within the scope officially accredited records 
the United States. The most outstanding example found the 
writer’s investigation (first reported 1930), covering period some five 
years but dealing with several times many stations were included the 


intensity would in. Therefore, the new curve straight 


final tabulation (39), was the record Taylor, Tex., 23.1 in., associated with 


maximum hourly rate 4.25 and monthly maximum 25.9 in. Similar 
maxima recorded for Austin, Tex., another tabulation more recent date 
are 19.0, 4.38, and 20.8 in., respectively (40). 

maximum observed values are desired for some the world’s outstand- 
ingly high rainfall stations, one must include 24-hr depths 21.9 and 28.7 in. 
for remote stations Hawaii and Jamaica, British West Indies, respectively, 
and 25.0 in. for Madras, India, and 40.0 in. for Cherra Punji, India, associated 
with hourly maxima somewhat less than 4.0 in. and with daily and monthly 
maxima 42.8 and 360.1 in., respectively (39)(40). 

Similarly, for the Baguio Station, high among the mountains Luzon, 
Philippine Islands, the maximum hourly rate was nearly 4.0 in., whereas the 
maxima for daily and monthly periods were recorded and published 46.0 
and 133 in., respectively. Obviously, the curve for century-frequency maxi- 


mum intensities that includes all known values should connect point the 


top margin with abscissa nearly equal 2.0 in. per with the origin; would 
therefore lie slightly the right the curve Fig. 14, representing the maxi- 
mum 4-hr intensities within the United States. The last-described curve, 
dealing with world maxima, should preferably excluded, since would dis- 
turb the orderly march the series developed for the United States. 

far the writer was able determine, much the concept relating 
rainfall patterns was new when was published 1930 (39). Two illustra- 
tions (39a)(39b), table (39c), and textual references pertaining all three 
incline toward the establishment nearly straight-line relationship for the 
series. For many the stations listed (39d), comprising 820 the principal 
station-records throughout the world summarized form, the writer reported 
that the mean annual rainfall divided the average number rainy days per 
year would represent the maximum rainfall depth expected within 
3-min period (39e). Obviously, this curve would lie somewhat the right 
the one Fig. and might well serve envelope the plotted values 
below the graph for 5-min duration. Inasmuch the majority the 
tions listed the writer (39d) for world-wide distribution had data for 
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periods shorter than one day, the intermediate values between the 3-min period 
and the 24-hr period had determined interpolation. For this purpose, 
the writer determined trial that curves for the 15-min and 60-min durations 
practically trisected the space between the 3-min and 24-hr maxima. For 
that reason, principally, the 30-min rainfall intensities were purposely omitted, 
being less significant members the series than were those retained. Other 
members the plotted series rainfall data (39a) were maxima for both 
monthly and yearly rainfall depths. Differences lengths records for the 
respective periods, from day less months and years, were inclined 
disturb the true relationships. was not uncommon find only decade 
rainfall intensity records for day less, whereas the monthly and 
yearly values published covered half century more, thus calling for 
some adjustment reduce such data basis comparable frequencies. 
Under the heading, Mr. Hathaway states: 


study rainfall intensity-frequency data for large number 
precipitation stations indicates that there are fairly consistent relationships 
between the average intensity rainfall for period and the average 
rates comparable frequency for shorter intervals, regardless the geo- 
graphical location the stations frequency 1-hr 


This terse statement far-reaching importance and somewhat startling 
first reading, links hourly and shorter intensities regardless geographical 
location and frequency each series. proves almost axiomatic the 
light further study. fact, part the concept rainfall patterns 


which the writer sought express 1930 (39) and earlier. effort 
define the essential similarities inherent precipitation intensities total 
depths for varying periods time from min year, the writer expressed 
all these quantities percentages the mean annual depth rainfall. Thus, 
the maxima for 5-min period ranged from 0.17 in. Juneau, Alaska, 
0.91 in. Denver, Colo., 1.05 in. Tampa, Fla., 1.15 in. Terre 
Haute, Ind., representing the maximum officially reported and observed depth 
the United States. Reduced percentages annual mean rainfall depths 
for the respective stations, the maxima become 0.2%, 6.4%, 2.1%, and 3.0%. 
Some other stations have established claims for somewhat higher records but 
have not fully qualified them official data. Undoubtedly still higher values 
have occurred for tropical storms, although official observations are lacking. 
Maps.—For several years after the publication technical re- 
ports, including the volume entitled Rainfall Eastern United States” 
(41), the Miami Conservancy District, many hydrologists considered that 
the Great Plains and Western States were not susceptible the same treat- 
ment. The reason for this opinion was related the erratic and unpredictable 
rainfall occurrence, amount, intensity, and areal extent, well the 
dominant influence altitude over latitude. the face such professional 
opinion, the writer undertook extend series isohyetal lines representing 
not only daily but also monthly maxima throughout the United States 
(399), and then completed the maximum monthly isohyetals for such parts 
the various continents were covered officially published records. Mr. 
Yarnell (15) adopted scales that precluded the portrayal local effects, such 
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are associated with topographic forms, rain shadows, normal paths air- 
mass movements, and the orographic barriers surmounted. over- 
come this difficulty part, Mr. Yarnell introduced special maps more 
generous scale, covering the Pacific Coastal States, with regard 24-hr rainfall 
depths for various frequencies, from 100 years. Lack data had been 
noted the Pacific Coastal strip for some the shorter periods—namely, 
16, and hr—and implication for still shorter periods. 

The writer’s experience the construction and drainage arsenal and 
airfield facilities lead him offer the following observations: 


(A) matter how liberally one may interpret the series isohyetals 
the small maps Fig. 13, there urgent need for either supplemental maps 
the Pacific Coastal region scales permitting the portrayal the pro- 
nounced local variations, special note caution concerning the practical 
certainty that the indicated limit 1.25 in. per has been and will greatly 
exceeded within 25-yr periods record. 


TABLE RAINFALL RECORDS THE PACIFIC 
REGION 


ANNUAL 
Maxi- 


Washington: 
Tatoosh 


North Head.... 
Oregon: 
Portland... 
Roseburg 
California: 
Monumental.... 
Crescent 


Point Reyes.... 
San Francisco... 
Sacramento 


Los Angeles. 

San Diego 

Campo 2,543 

Squirrel 5,280 
(2.0) 


ono 


Parentheses enclose unofficial estimated quantities. All units Northwest. 


(B) study the selected station-records Table (which were 
intentionally represent the upper bracket observed intensities, 
most the actual maxima high mountainous regions are unrecorded) 
indicates that the influence latitude subordinated that altitude 
bined with windward exposure. 


| 
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Item Altitude 
Station (ft) 
a eee 480 (2.3) (6.6) 31.3 108.4 120.3 153.8 
300 (2.0) (6.0) 24.6 79.6 124.4 141.0 
300 (2.4) (7.0) 42.3 79.2 129.5 157.1 
‘b Seattle......... 246 0.80 3.52 15.3 20.0 33.8 56.4 
109 0.79 3.79 18.9 22.0 39.9 57.3 
211 0.86 4.29 17.8 39.9 55.5 78.6 
(2.2) (6.5) 37.2 63.6 95.9 132.1 
(2.5) (7.0) 29.3 82.5 109.3 134.4 
575 (2.7) (8.0) 50.2 97.3 130.5 167.3 
1.31 7.70 20.1 30.8 42.8 67.2 
510 1.14 4.10 12.8 21.1 33.3 44.3 
2,700 (3.0) (9.0) 43.8 153.5 
125 (2.6) (8.0) 31.9 107.6 
; 15 eS Se 64 1.30 6.87 19.5 2 64.5 
Red Bluff....... 307 1.72 6.10 20.7 49.0 
490 1.40 3.10 11.4 32.4 
206 1.07 4.70 24.4 38.8 
295 1.36 2.90 9.5 21.5 
8,060 (2.8) (8.5) 39.8 94.3 
7.36 15.8 40.3 
3.20 9.3 27.6 
(8.6) 20.4 35.3 
(6.0) 29.0 125.5 
tic 


(C) The seaward slopes the coastal ranges extending from the Olympic 
Peninsula Washington San Diego, Calif. (and probably even farther) and 
from the western slopes higher ranges farther inland (such the Cascade 
Range and the Sierra Nevada) are capable developing fair approach the 
series intensities plotted around the Gulf Coast. These intensities approach 
about in. per outstanding stations, although they are ordinarily limited 
roughly half that depth for available periods hourly record. 

(D) Table the rainfall patterns disclosed neighboring stations are 
used for practical purposes estimate appropriate values the missing items 
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basic data. All these are enclosed parentheses indicate their unofficial 


character. These estimated values may subject uncertainties approxi- 

mately 20%. 

Both limiting and periodic intensities for selected rainfall stations seem 

vary with the square root the number minutes the respective time 
intervals considered, has been stated official publications (42)(43). Thus, 
the envelope conceivable intensities—perhaps far beyond the attainable 
limit for most the outstanding stations high rainfall—would have depths 

Minutes Inches 

43,200 (month days)............. 207.4 

525,600 (year 365 days).............. 725 


Some the foregoing depths have been attained even exceeded—for 
example, the daily, monthly, and yearly depths Cherra Punji and the daily 
and monthly depths Baguio. Some have been either attained closely 
approached several stations India and the tropical islands the Pacific 
Ocean. During seven the fourteen years available record for Cherra 
Punji (1905 1919), the maximum daily rainfall depths were excess 
in.; and three these exceeded in., follows: 33.00, 31.05, and 32.35 in. 
during 1910, 1914, and 1915, respectively. Desert regions may attain much 
greater proportion the short-period rainfall intensities, say, for less than 
hour, than should expected for days months. Otherwise, they would 


cease classified deserts, inclining more toward the humid character. 
Patterns for most stations will follow more less closely some proportion 
the foregoing envelope depths. 

Runoff the interest simplified procedure and computa- 
tions, Mr. Hathaway has wisely made number assumptions, some which 
are admittedly only approximations cover variable conditions and behavior. 
(Either compensating cumulative errors may result for given problem.) 
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Among these are the assumed constancy rainfall, infiltration, and resulting 
rate supply for overland flow during specific storm. Inasmuch the 
assumed rectangular form pluviograph seldom ever attained and 
alternate form isosceles triangle the same time base but with its altitude 
double the height the rectangle would enclose equal area and would 
measurably alter the proportions for infiltration and supply volumes, this 
alternate form pluviograph may bear investigation. Although would still 
approximation, would occupy mean instead extreme position. 

Mr. Horton’s equations for overland flow, presented Mr. Hathaway, 
well single page textual material Mr. Yarnell’s notable bulletin 
(15), employ the so-called for runoff, which 
expresses runoff cubic feet per second, representing some percentage 
decimal fraction (C), the ratio runoff rainfall for the critical period, times 
the product rainfall intensity inches per hour and the number acres 
the drainage area. The infiltration capacity likewise expressed inches 
per hour; as, for example, 0.6 in. per adopted Mr. Hathaway for turfed 
areas. Similarly, the allowable pond storage 2,000 per acre the 
example used Mr. Hathaway may expressed depth 0.55 in., since 

For the example solved Mr. Hathaway according the rigorous pro- 
cedure recently developed with discharge the inlets the drainage system 
given 2.1 in., from rainfall intensity 3.6 in. per hr, equal 
0.583, 58.3%. Evidently, the surface storage outside ponds was 
0.35 in., which added 0.60 0.55 2.10 3.60. Only that portion the 
surface storage (10% the total rainfall) subject estimation; could 
not seriously error; neglected entirely, the error would only 10%. 

Perhaps just plain coincidence, but interesting note that the 
X-ratio value derived graphically Mr. Hathaway 59% (compare with 
58.3%), representing the ratio drain-inlet capacity peak rate runoff, 
the difference being absorbed pond storage. 

The impressive array well-chosen tables basic data and related com- 
putations, together with the charts and diagrams which the illustrative 
solutions typical problems are readily traceable, and the clear-cut, logical 
treatment the entire procedure the accompanying text Mr. Hathaway 
leave but little supplied discussion. gratifying know that the 
design procedures have already been incorporated the Engineering Manual 
the Construction Division, Office the Chief Engineers (23). 

The susceptibility the entire analysis and rigorously derived results 
comparison and expression terms the so-called “rational formula,” 


acre-inch 


demonstrated herein, should insure their incorporation professional practice 
much more widely than would have been possible without such popularly 
established device. This formula results closing survey circuit well within 
the tolerance limits, when results derived, capable derivation both 
methods, are found substantial agreement. Mr. Hathaway’s paper 
represents noteworthy accomplishment and progress the field hydro- 
meteorology. 
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splendid method for rational design drain- 
age capacity contained Mr. Hathaway’s paper. The engineer given 
measuring stick which may correlate the effect the storm deemed 
critical for particular installation with the airport topography and the drain 
capacities necessary for the required use the airfield. Such collection 
concise data and workable information appreciated the drainage engineer. 

Proper use subgrade drainage may often simplify surface drainage facili- 
ties. Where extensive flat areas are involved, better over-all solution the 
problem may sometimes obtained combination several methods. 
The following analysis and comments give some reasons for serious considera- 
tion subgrade drainage facilities. 

The 1-hr infiltration value subtracted from the 1-hr rainfall quantity iden- 
tifies the standard supply curve and forms the basis for succeeding computa- 
tions. the example quoted Mr. Hathaway’s computation storm drain 
capacities assuming supplemental ponding, the design storm chosen presumes 
2.2 in. rainfall per hour. The 1-hr infiltration capacity for turfed areas 
assumed 0.6 in. per and represents 36.7% the rainfall quantity. 
Surely this important enough warrant further investigation and more 
accurate data infiltration capacities. The pertinent factors are many and 
their interrelationship apparently vague. Problems presenting greater diffi- 
culties have been solved previously through research. 

36.7% the rainfall quantity may disposed through infiltration, 
why not either dispose the water quickly through subdrainage facilities and 
thereby reduce the amount ponding, or, lowering the adjacent water 
table, create ponding area below the surface and away from traffic. 

The possibilities subdrainage pipe seem wide and valuable. Quick dis- 
posal filtered water will reduce overland flow and minimize surface erosion. 
The subdrainage system creates storage area below the plane traffic, 
thereby reducing the number required inlets with their attendant safety 
hazards. 

Lowering the water table reduces the damaging action frost subgrades 
under runways and aprons and, through better stabilization, can reduce the 
cost surface treatment. proper subdrainage system will make available 
for landing outside runways areas which might never attain usable condi- 
tion the entire disposal depended inlets and storm sewers. 

view the many facets the subdrainage problem, hoped that the 
same clear, usable type analysis described Mr. Hathaway for storm sewer 
disposal systems may evolved for the design and efficient use subdrainage 
pipe. 


engineer who any way connected with 
the design construction airfields, whether for military use civilian use, 
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should find this important and timely Symposium indispensable. contains 
many sound principles utmost value the designer airfields, presented 
simple, straightforward, and stimulating manner. 

connection with airfield design and construction problems, the Corps 
Engineers has proceeded quietly and systematically designing, constructing, 
and testing, meet adequately the demands the airfield war construction 
program. Although greater and greater wheel loads have been predicted, com- 
paratively little information has been published methods proposed meet 
such airfield problem. 

Base-Course Materials—In discussing base courses for con- 
crete pavements, Colonel Stratton uses the term “high-quality base materials” 
without sufficient definition. Because common misunderstanding 
what constitutes base and because the belief 
many engineers that clay should added granular materials for use 
base course (simply because the satisfactory performance has been observed 
well-graded sand gravel mixtures containing clay when used surface 
courses exposed directly light wheel loads without the use substantial 
pavements), should stated that, general, clean, well-graded granular 
base materials with minimum say, not over passing 
No. 200 mesh, make the most desirable high-quality bases for both rigid and 
flexible pavements. 

Only weight finer than 0.02 well-graded granular material 
used base course northern areas may cause troublesome frost action— 
breaking concrete due irregular heaving winter and insufficient support 
due softening the spring. 

southern regions, natural clayey gravel sand-clay mixtures may give 
satisfactory results placed properly compacted, homogeneous condition. 
However, the deliberate addition clays clean, well-graded granular ma- 
terials for use base courses should discouraged, not only because the 
difficulty, expense, and length mixing process, but also because the inherent 
increased sensitivity moisture control for proper compaction, well 
because the increased potential susceptibility subsequent moisture 
changes. 

Design Concrete conditions have restricted the use 
steel which turn has caused the revision concrete pavement sections. 
More than one good engineer pleaded that the reinforcing steel left pave- 
ments, and especially those which had placed directly poor subgrades 
without the benefit intermediate base courses. 

There question about the necessity for thickening the edges concrete 
pavement done the revised sections the pavement placed directly 
poor subgrade. However, would very desirable make additional 
investigations the necessity for thickening any edges when the pavement 
placed high-quality, well-compacted base course, considering the 
minimization differential settlement, together with the loading frequency and 
the natural traffic distribution airfields. not probable that the effect 
disturbing the base course shaping for thickened edges counteracts 
benefits obtained from the thickened edges? 
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Control Frost discussing pavement design for frost action, 
Colonel Stratton states that the most generally accepted method insuring 
pavement base-course stability over potential frost-action soil the pro- 
vision sufficient thickness insulating material. 

The writer would question the mere use over 
potential frost-heaving soil insure stability unless the term “insulating 
material” were defined granular material containing more than small 
percentage material finer than 0.02 mm. Without such restriction, one 
can well imagine granular material with, say, 10% finer than 0.02 which, 
were possible keep dry, would act insulating material but which, 
because its high percentage and high potential capillarity, might 
cause serious damage from frost heave used cover another potential frost- 
action soil contact with water supply. Therefore, using soils insu- 
lating materials, necessary place limit, not necessarily the grain- 
size curve itself, but the allowable amount soil finer than certain size, 
unless the soil can entirely separated from water supply either using 
intermediate impervious course selecting the soil that the limit 
rise will not fall within the zone frost penetration. 

Subsurface the past much time and money have been wasted 
installing subsurface drainage systems which, Colonel Stratton states, 
could frequently have been eliminated providing adequate and properly 
designed surface drainage system. 

subsurface drainage system referred herein system installed 
primarily for the purpose removing water from pavement base courses 
subgrades—some form improvised drain,” system open-jointed 
pipe, assemblage well points. 

When studying subsurface drainage, expedient classify soils 
one the three following groups: (a) Well-draining soils such clean sands and 
gravels; (b) Poorly-draining soils like the slightly plastic organic and inorganic 
silts; and (c) Impervious soils—for example, homogeneous clays. (These re- 
marks drainage are based lectures Arthur Casagrande, Am. Soc. 
E., the “Control Soils Airport Construction,” delivered the 
Graduate School Engineering, Harvard University, Cambridge, Mass.) 

Colonel Stratton states that where the ground water above the 
frost line causes soft, boggy conditions areas where soil stability re- 
quired, subsurface drainage should noted, however, that 
(a) lowering the water below the frost line may may not prevent frost heave, 
depending the height capillary rise, and (b) may much easier and 
more practicable some cases lower the water table using circum- 
ferential ditches providing drainage through underlying impervious 
stratum when such favorable condition exists. 

For clean well-draining granular materials, subsurface drainage can often 
eliminated altogether, least reduced minimum providing circum- 
ferential ditches intercept ground water from surrounding areas and 
sealing off the surfaces reduce the infiltration water into underlying courses 
other words, the water should kept from entering the base 
courses, but, some does penetrate, the perviousness these courses should 
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advantageous removing the water. the well-draining materials contain 
appreciable amount fines, then, prevent possible serious damage from 
frost heave, will necessary keep the water table not only below the 
frost line but distance below the frost line greater than the height capil- 
lary rise for the particular soil. 

strict sense, the statement that poorly-draining soils are “drained” 
with subsurface drainage system incorrect. The height capillary rise 
such soils usually great, and long the water table lies within the 
height capillary rise the soil will remain saturated capillary water. the 
frost zone comes within this zone capillary water, there may danger 
serious frost heave. Therefore, there advantage attempting use 
subsurface drains eliminate danger frost heave poorly-draining soils, 
unless the drains can placed such depths that the piezometric water table 
will lowered distance below the frost line greater than the height 
capillary saturation. For some soils, this impracticable impossible. 
However, benefits can derived from using subsurface drains poorly- 
draining soils. their use soft soils can consolidated and thereby their 
bearing value can increased. 

using subsurface drains for consolidation purposes the following should 
emphasized: 


(a) Since any benefits consolidation are obtained from capillary forces 
(that is, long the soil remains completely saturated), the deeper the drains 
placed within the zone capillary saturation, the greater will the 
force producing consolidation. 

(b) The closer the drains can placed for given depth, the less time 
required consolidate the soil between the drains. 

(c) function properly the system must designed that water can 
removed fast collects. may possible accomplish this 
gravity system certain favorable locations, but, more often, will 
sary use pumpsinsumps. Pumps are also needed when well points are used. 

(d) Drain pipes.for such system may very small since the actual 
water moving through them may more than (The 
lying soil remains saturated capillary water that the greater part 
from that source will only that which squeezed out consolidation) 
Since the soil nature drains poorly, only small quantity water can 
expected from infiltration and ground-water movement. 

(e) The ideal situation would one which the poorly-draining soil 
underlain pervious layer, which case water could removed from tht 
pervious layer means well points circumferential ditches. The 
subsurface drainage system homogeneous clay extremely slow that 

little value for improving the condition the clay. 


Construction designs may vitiated poor 
and should stressed that high-quality construction control and 
are essential many types construction. 

The value subgrades and bases more and more 
excellent pavement may deteriorate inferior base and 
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whereas second-grade pavement properly designed and constructed 
foundation may give excellent service. engineered foundation requires 
proper control from design through completion. requires careful selection 
and placing materials, proper moisture control, and adequate compaction— 
the real value which cannot overstressed. Many the defects settle- 
ment, swelling, and shrinking, often blamed the pavements, can 
traced the foundations and could have been minimized proper compaction 
control. 

Remarks Training Soils teaching Army courses 
“The Control Soils Military Construction,” the writer stresses particularly 
the value identifying soils correctly because thereby the possibility making 
major mistake may eliminated. After soil correctly identified, the 
properly trained soils engineer and his inspectors should know how improve 
that type soil, how judge its drainage characteristics, how decide the 
degree moisture control necessary for efficient compaction, and how 
determine what method should used for compaction. 

prewar days, too little attention was given the proper training 
field men and inspectors for all types highway and airfield construction. 
Frequently proper control and inspection was sacrificed for impressive looking 
plans. the postwar period bring extensive highway and airfield con- 
struction with ample available man power, anticipated, then surely some- 
thing will have done assure not only the proper utilization experi- 
enced personnel returning civilian work, but also the provision ways and 
means for the training additional men. The time consider these problems 
has already arrived. 


facilities—roads, bridges, railways, and waterways—has been age-old source 
problems the civil engineer, and now the construction airports, with 
scope undreamed only few years ago, may well making important new 
chapter civil engineering. rapid developments take place, reliable 
published information often lags far behind. This Symposium fortunate 
exception. 

The following discussion deals primarily with the first paper the Sym- 
posium. Colonel Stratton has succeeded remarkably well summarizing all 
pertinent information and essential considerations this subject, with par- 
ticular emphasis the many important contributions which the Army Engi- 
neers have made the design and construction airfields for heavy plane loads. 
Although restricted scope military airfields, should emphasized that 
most the Symposium also applies the commercial airfields the future. 
fact, almost all the extensive research the design and construction 
airfields, conducted under the direction Colonel Stratton, will inesti- 
mable value for the postwar development commercial air facilities serve the 
air liners the immediate future. Therefore, hoped that this paper will 
followed others describing the results the large-scale traffic tests which 
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have been undertaken the Engineer Department, the results which 
will indispensable every engineer designing pavements for heavy planes. 

Design Flexible design procedure adopted the Office 
the Chief Engineers based the empirical method developed the 
California State Highway Department. stated Colonel Stratton, the 
design curves Fig. have been confirmed large-scale traffic tests. 

may well explain why has not yet been possible develop 
successful theoretical approach the design flexible pavements, based the 
theory elasticity. Perhaps the most important reason that the coarse- 
grained materials, which pavement and base courses are generally composed, 
not follow Hooke’s law. Therefore, the distribution stresses the 
subgrade will show considerably greater concentrations than that computed 
the theory elasticity. addition, most flexible pavements not have 
enough tensile strength resist the considerable tensile stresses which develop, 
according the theory elasticity, for large ratios the modulus elasticity 
pavement-plus-base over that the subgrade. Even tension cracks 
develop, the computed stresses the subgrade would the unsafe side be- 
cause the deformation the pavement tension much greater than that 
compression; other words, for practical purposes the pavement cannot take 
tension. These remarks are not intended discourage further attempts 
theoretical studies the design flexible pavements, but merely emphasize 
difficulties which are often overlooked and which can lead conclusions very 
much the unsafe side. 

Stabilized Versus Compacted confusion still prevails the 
the construction base courses. cite two extremes: (1) Excellent concrete 
pavements have been built Long Island hydraulically placed, clean, beach 
sand, without even attempt compacting the sand; whereas (2) different 
place, well-graded, fairly clean, glacial gravel (the best base material next 
crushed stone) was adding 15% material passing No. 200 mesh, 
thereby changing into frost-heaving material. Even where there frost 
action, the addition plastic binder gravel has often resulted type 
base which, condition capillary saturation (evaporation from its surface 
prevented the overlying pavement), was much less resistant deformation 
than the coarse material without the binder. Not only costly mix 
binder with coarse material, but difficult make the mixture uniform. 
Even small excess clay binder may cause local destruction pavement, 
which would not happen the admixture such binder were omitted. 

Instead paying for the admixture fines, time and money should 
directed toward very thorough compaction the base-course materials. 
though has been emphasized, far too little attention paid proper com- 
paction. fact, some engineers really believe that uniform sand cannot 
compacted; that always remains loose. (Perhaps this belief due the 
ambiguity the term which also used designate grains with 
cohesion.) The writer foresees the development compaction equipment 
which means vibration will produce degrees compaction coarse- 
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grained soils unattainable with any methods now used, and little 
additional cost. 

The importance thorough compaction best illustrated the fact that 
well-compacted sand has about three times more resistance deformation 
and bearing capacity than the same sand its ordinary loose condition (not 
bulked). For well-graded gravel this ratio may high fifteen. The 
large increase resistance deformation gravel surfaces exposed constant 
traffic well known. Such surfaces may not satisfactory unless traffic 
compacted. Does this not prove that possible compact such materials 
far better than ordinary construction methods? What can done 
traffic can accomplished vibration equipment especially designed for the 
compaction sandy and gravelly soils. 

Considering that the cost compaction only fraction the cost the 
material itself, particularly where long hauls are necessary, only logical 
improve the material compaction the utmost limit. Manufacturers 
construction equipment would well pay attention this problem and 
its partial solution already attained before the beginning World War (44). 

Design Concrete Stratton advances pertinent argu- 
ments for restricting the thickness concrete pavements about in., and for 
providing the necessary support poor subgrades appropriate increase 
the thickness the gravel base. Success such design will depend the 
uniformity subgrade support. When dealing with nonuniform support such 
results from irregular subgrade conditions, soils subject swelling and 
shrinkage, nonuniform frost action, spite good base times there may 
little support, and corners may break under the design load, even less load. 
For such unfavorable subgrade conditions, the use thicker slabs with steel 
reinforcement cannot avoided. Warping stresses might controlled 
sufficiently small size the slabs, say, for pavement in. thick. 
The wire mesh reinforcement will preferably used two layers, spaced about 
from the top and bottom Some additional reinforcement the 
top the pavement should provided the corners. such design the 
steel would actually participate carrying stresses, instead acting merely 
tie broken slabs together. 

thick pavements even greater attention will have paid the design 
the joints for very heavy loads. The most desirable solution would 
make all joints alike, each one permitting slight amount expansion. 
trapezoidal tongue and groove design could used with premolded asphaltic 
boards, about 3/16 in. thick, containing enough air voids (for example, form 
the total expansion possible for 100 pavement, assuming 12-ft spacing 
the joints, would about This allowance ample view the uni- 
form distribution this expansion and the greater assurance that all joints will 
remain operative, compared with the standard, wide expansion joints that 
have been known become partly inoperative with time. Heavy, short steel 
dowels should provided, painted with asphalt and having thin compressible 
disk glued each end, allow for the necessary slight amount movement. 
vertical steel pins are used hold the premolded boards place (piercing the 
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trapezoidal tongue and being driven into base subgrade), the section inside 
the tongue may require sleeve cover allow movement between this part 
the pin and the concrete. 

uniform joint construction, distributing expansion uniformly over all 
joints, would aid greatly solving many undesirable features the concrete 
pavements the past. 

Insulation for Frost the feasible methods for protection 
against frost action the use insulating material beneath the pavement, 
such thickness prevent frost penetration into the subgrade. Where 
cheap timber available and gravel scarce, timber deserves serious considera- 
tion, because thickness in. will offer the same protection in. 
gravel. 

Another material which the writer has seriously considered for this purpose, 
ever since became acquainted with 1929, porous concrete made from 
mixture cement mortar and foam. This Danish development has been 
used Europe and South America for heat insulation and other construction. 
Fracture surfaces have the appearance slag. can manufactured 
different unit weights and insulating characteristics, and has the important 
advantage that can mixed the field and poured like concrete. This, 
some other cheaply made material full air voids and having compressive 
strength least 300 per in. (which will not soak water into the air 
voids when submerged), poured either place pre-cast slabs, would solve 
the problem frost action cases where the necessary thickness gravel 
would prohibitive cost. areas permanently frozen ground such 
insulating materials could serve protect the frozen ground under pavements 
from melting. 

Wheel Loads Airplanes.—Airplane designers should realize that placing 
the entire weight the heaviest planes only two wheels, two pairs 
dual wheels with little space between the dual tires, creates the necessity very 
expensive airfield pavements and limits the number airfields which such 
planes can land safely. For the heaviest planes the load should spread over 
four wheels, with the dual tires spaced least apart. For further increase 
plane load 150 tons, the writer foresees the necessity spreading the load 
over two four-wheel carriages, resembling perhaps the wheel carriages railroad 
cars. With such eight-wheel landing gear, 150-ton plane would have the 
enormous advantage that emergency could land ordinary 
mercial airfield without danger breaking through the pavement and serious 
consequences for plane and passengers. Civil engineers specializing airfield 
construction should consulted connection with the design the landing 
gear air liners. 

Closing Remarks.—The writer takes this opportunity express his apprecia- 
tion the broad application soil mechanics the design earthworks, 
foundations, water-front construction, and now airfields, which the Army 
Engineers have followed during the past decade, and which Colonel Stratton 
has played leading The Army Engineers have contributed more the 
advancement applied soil mechanics than any other governmental private 
organization. 
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Harry Am. Soc. Corps Engineers, with 
which the two authors (Colonel Stratton and Mr. Hathaway) the Symposium 
under discussion are connected, has contributed much the various problems 
airport design and construction—especially base and subgrade design, 
proper backfill for subdrains, and runoff criteria for the design storm drains. 

Backfill for Subdrains.—The design base-coarse drainage Fig. 8(b) 
Colonel Stratton’s paper shows the perforations placed which contrary 
good practice especially proper backfill (filter material) used. This sug- 
gests that discussion order backfill, the proper design the infiltration 
area (perforations), and the amount infiltration area required. 

The chief enemy subdrains clogging either the pervious trench 
backfill the pipe itself 
been solved investiga- MATERIAL 
tion undertaken the 


Standard A.S.T.M. Percent. 
Waterways Experiment Station 
Vicksburg, Miss. (45). The 100 
remove ground water from 
water bearing zone. Such zones No. 
No. 


fine sands, coarse sands, and 
gravel. prevent the en- 
trance the fine soil particles into the trench backfill, the investigation shows 
that the backfill should not coarser than that Table 10. 

clean backfill having grain-size curve between the curve grading 
Table and that medium sand (No. sieve) will satisfactory. This 
radical departure from the usual state highway specification for material 
mostly between in. and Some states allow 10% passing No. 
sieve and some allow material large The idea that coarse backfill 
needed with pipe subdrain probably originated with the requirements for 
drain (without pipe) which large voids are necessary for lateral 
flow. The vertical flow capacity medium sand ample for all ordinary sub- 
drain conditions. The pipe provides lateral flow. 

Such fine grained backfill not only necessary prevent clogging but, 
many localities, more economical. However, necessary that the drain 
pipe designed eliminate the entrance the backfill and yet allow for the 
entrance water. 

Because most the backfill material smaller grain size 
than the perforation diameter, drops through the holes gravity action 
when the perforations are the upper half the pipe. This fact was demon- 
strated the Vicksburg tests well all other similar tests previously 
conducted. Therefore, instead the perforations being up, shown 
Fig. they should down, the upper part the lower half the pipe. 
Also, prevent the backfill from entering, due side pressure, the upper 
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row holes should not higher than 22.5° below the horizontal axis the 
pipe shown Table 11. 

There seems undue concern among engineers regarding the loss 
water once has been collected inside the pipe. There are few instances 


TABLE 11.—RECOMMENDED PERFORATION DESIGN 


Water Table 


Tamped 

Material 
(a) FOR BASE DRAINAGE FOR 
AND LANDSLIDES FOR SUBDRAINAGE 


Number of: 


Section (a) 


Perforations per linear foot, averaged per section pipe (minimum)....... 


Section (b) 


subdrainage work where such concern justified. For the usual condition 
lowering water table saturated soil (section (b), Table 11), the zone under 
the pipe already saturated and therefore loss water will occur. How- 
ever, for landslide and terrace bank drains, collect seepage that may 
flowing top thin impervious zone, the bottom the pipe may lie 


Runway Transition Strip Shoulder 


Bituminous Bound Rock (1'' Size) 
Pervious But Stable Under Traffic 


Compacted 
Stabilized Foundation 
Subgrade 


Size) 


Depth 
Variable 


Tamped Impervious Material 
Minimum Outside Diameter 
+6" 


pervious soil and precaution should then taken prevent loss water 
methods such those shown section (a), Table 11. For these cases and for 
base drainage, the quantity water collected small and two rows 
perforations are more than adequate. When drain pipes are installed adjacent 
pavements impervious soil zones for base drainage and also for intercepting 
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drains, impervious trench bottom, shown section (a), Table 11, and 
Fig. 38, advisable insure complete removal free water collected the 
trench. The backfill ‘‘filter over the top and sides should 
tamped mixture sand and gravel not coarser than the grain-size curve 
material Table but having permeability sufficient for runoff. 

Required Infiltration Area (Perforation hydraulics subdrain- 
age not exact science. The engineer may guided research 
ground-water flow, but the answer largely based judgment and experience. 
(Typical such research are the runoff tests the Engineering Experiment 
Stations Iowa State College Ames and the University Minnesota 
Minneapolis. Also the Darcy formula may used: 


section area flow.) The ground-water flow seldom sufficient fill 6-in. 
pipe within 1,000-ft run. Even for the most severe seepage condition en- 
countered highway and railway cuts, 6-in. pipe has sufficient capacity the 
first 400 500 the run. change pipe size the run (for example, 
in. in.) required carry the accumulated interception but not because 
increased ground-water flow the drain per linear foot. Therefore, the 
larger drain pipes need more infiltration area (perforations) than required 
6-in. pipe. 

Little known about the infiltration capacity perforations pipe 
surrounded with sand sand and gravel. One such test with metal pipe 
gave approximate answer 0.009 per sec per in. perforations 
holes) for condition which little head was assumed above the 
top the pipe. this basis, 6-in. pipe, 1.0% slope with perforations 
specified section (b), Table 11, could fill itself less than ft. When 
consideration given the fact that the area perforations run 
specified section (b), Table 11, equal the end area 6-in. pipe, 
self-evident that the infiltration area many times larger than necessary; 
other words, large factor safety provided. When too many rows 
perforations are specified pipe, requiring the top row near the 
spring line (horizontal axis), solids are likely enter. prevent such en- 
trance one the requirements well-designed subdrain pipe. 

Impact.—Colonel Stratton’s discussion loading criteria leads the reader 
conclude that impact not factor the design runway pavements and 
may even construed infer that impact does not occur. The subject 
impact has probably received more discussion than any other one item 
connection with airport design. The justification for designing runways for 
wheel loads without impact seems based mostly the observation that 
runway pavements show less damage than apron pavements used for warming 
airplane engines. The writer the opinion that the apron failure not due 
increased load, but mostly the failure the foundation, which becomes 


puddled vibration, creating condition similar that beach sand when 
tapped with the foot. 
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can assumed that there impact the instant perfect landing 
but, this date, published data indicate definitely the existence impact for 
the average landing: Evans, the Goodyear Tire and Rubber Com- 
pany, has stated (46): 


connection with this work [landing tests determine the efficiency 
tires landings variety soils and used two planes, 
one the class and one the class. They were equipped 
with recording accelerometer whereby the actual impact factors developed 
during landing could determined. result this work, arrived 
the following conclusions: Ordinarily our typical landings made hard 
surfaces developed impact factors ranging from 1.25 1.75 Anything 
less than 1.25 would rated easy soft landing and one above 1.75 
2.0 would called hard landing. The hardest landing which 
ever experienced all this series tests was The design factor 
both these planes was, recall, 6.5 that still had quite ways 
before approaching the limit impact values which the planes were de- 
signed withstand. this connection, may interest point out 
that all types airplane tires, operated recommended values load 
and inflation will not ‘bottom’; that is, deflect down contact with the rim 
until the impact load reaches 3.5 4.0 times the static load xx.” 


contrast the Goodyear impact tests with light planes, the Douglas 
Aircraft Company conducted tests (47) DC-4 weighing 63,000 lb. The 
normal landing impact factor was 1.43 and the extreme landing impact factor 
was 3.5. Although there was wide difference the weight planes used, 
the impact factors developed both tests are close agreement. 

basis for the design the primary airplane structure, the Civil Air 
Regulations (48) require, addition factor safety, depending the 
part considered, limit load factor for level landings and three-point landings. 
For planes heavy 2,000 gross weight this factor 4.33, whereas for 
planes heavier than 2,000 lb, 


which the gross weight pounds. For 74,000-lb plane 
static wheel load) the factor becomes 2.91, which added factor safety. 

summing this impact question, should remembered that the 
ground reaction the sum total the forces causing tire deflection, shock 
absorber travel, and stresses the chassis. Also, whenever the tire deflection 
the shock absorber travel greater than that under static load, impact 
must exist. 

The California Bearing Ratio method now used extensively for the design 
airport pavements and pavement bases. This method empirical, based 
California highway experience, and doubt contains impact factor for 
highway traffic, the existence which was definitely proved tests con- 
ducted the Iowa Engineering Experiment Station. these tests trucks 
moving fast miles per were used, and transmitted loads were scale 
weighed. The resulting for the design highway structures, 
including culverts, was that impact factors from 1.5 
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fore, would seem preferable state that impact allowed for the 
fornia method rather than that impact used, infer that normal 
landings create impact. 

The time factor another item that has received considerable discussion. 
The contention some engineers seems that load applied and removed 
suddenly does not produce the same transmission effect the pavement founda- 
tion buried structure would static load slow moving load. 
Whether this true, what extent the transmitted load may lessened 
supersudden load, compared sudden load, not known. Ob- 
servation has shown that mud pumping occurs from slab deflection under fast 
moving trucks well under slow moving trucks. Also the deflection 
large flexible culvert under heavy railroad loading with earth cover 
has been observed. this case distinct deflection occurs with the passing 
each the front and rear group wheels and with deflection recovery be- 
tween the application these loads, thus indicating the speed with which 
moving loads are transmitted through earth. 

Surface Hathaway’s criteria for surface drainage based 
entirely catch basins inlets. Two methods have been adopted for 
the collection surface water—the catch basin (or inlet) and the surface inter- 
cepting drain. combination has been used. Experience 
has shown that both methods have their advantages and disadvantages. 

The intercepting was devised because the apparent 
difficulty concentrating surface water one point without the use 
pronounced gutters ditches, which are objectionable airport surfaces. 
The idea intercept the surface water without changing the uniform sheet 
flow piled-up flow; and, when properly functioning, these drains just 
that. The two chief difficulties have been silting the top layer and dis- 
placement the top layer traffic. These difficulties are not the fault the 
idea, but are due improper use and construction. Silting will occur when 
the surface water carries silts and fine sands; and, therefore, when used 
indicated Fig. 38, all turf areas should slope away from the drain. Regard- 
ing displacement, the tendency use too little bitumen. The correct 
amount can determined best making tests with the particular stone 
used. Contrary the opinion asphalt technicians, the claim some- 
times made that enough bitumen used stabilize the surface, porosity 
lost. This problem needs and deserves much more research work. 

Catch basins (or inlets) function more satisfactorily sumped areas than 
any other type. Those ordinary size are not efficient collecting water 
from valley surfaces where warping limited airports. This fact 
recognized those who use extra long grate-top inlets across valleys runway 
pavements. These valleys are generally located outer part the pave- 
ment shown Fig. 8(c) collect the runoff from the runway and thus 
prevent erosion and softening the shoulder. Catch basins are even less 
efficient collecting surface flow turf valleys. combination the 
catch basin and surface intercepting drain system shown Fig. would 
insure much better interception. For this use, the intercepting drain backfill 
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the surface should not coarser than the grain-size curve the filter 
material Table 10. newly graded fields, silting may, and probably will, 
occur, but, using the recommended backfill, the silt and fine sands, carried 
the surface water, will stopped approximately the upper in. 
This can bladed off and replaced with clean filter material. Under ordinary 
sod-growing conditions, such replacement will not required more than once 
twice. 


Surface Flow 


Perforated Pipe 


Variable 


Catchbasin Inlet 

Minimum Outside 

(a) Diameter 


(b) SECTION A-A 


Fie. AND INTERCEPTING DRAIN TuRF VALLEYS 


Combined Surface and Subdrain (or Base Drain) combined 
system, shown Fig. 38, has received criticism the ground that moisture 
may introduced into the subgrade, water may forced into the pervious 
base, and when the intercepting drain becomes full water. For storms 
greater than that for which the drain pipe designed (frequencies year, 
years, years, etc.), the trench backfill may become full water, but theo- 
retically only once the design frequency period and then only for few 
minutes. considering this question, one must remember that water follows 
the line least resistance and will not enter the pervious base unless under 
head—a condition infrequent negligible. for the introduction 
moisture into the subgrade capillary action, the design provides for com- 
plete removal free water the trench, and the remaining moisture will 
probably more than normally contained the subgrade soils. 

For economic considerations backfill, etc., the accumulated water the 
intercepting drain pipe usually diverted into main sewer when pipe larger 
than in. required. 

The realization the need for pervious bases and base drainage increas- 
ing. The proper place for base drain the edge the pavement. Also, 
prevent damage runway shoulders, the proper place intercept surface 
water the edge the pavement. combined drain that performs both 
functions usually more economical than separate system. However, 
other engineering problems, the conditions must govern the decision 
whether the combined system separate systems are used. 

Runoff Computations—The computation runoff from various types 
sizes drainage areas has always been unsatisfactory task. two engi- 
neers working independently will arrive the same answer. This fact 
probably due more the values selected for the variable factors than the 
formula used. 


Sand Gravel Backfill 
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Most runoff formulas contain “runoff factor” which expressed 
percentage rainfall. The Horton formula (Eq. applies this factor 
the rainfall less detention. This seems step the right direction 
self-evident that for like surface conditions (soil, vegetation, slope, etc.), 
the percentage runoff would greater for 2-in. rainfall than for 1-in. 
rainfall. the other hand, the engineer confronted with the task assign- 
ing value for hourly detention which there seems information. 

Mr. Hathaway has done commendable job simplifying the Horton 
formula computations the use charts, but the solution still more 
laborious than most other formulas. However, after sifting out the es- 
sential steps, Mr. Hathaway’s solution less difficult than appears first. 

The writer has not made sufficient number computations support 
definite statement obtained the Horton formula compared 
with runoffs obtained other formulas, such the Burkli-Ziegler and the 
rational formulas. many cases the answers are close agreement provided 
the same runoff factors are used. 
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AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 1852 


DISCUSSIONS 


SHEAR AND BOND STRESSES REINFORCED 
CONCRETE 


Discussion 


GRINTER, AND ANDERS BULL 


Am. Soc. authors approach their problem 
with the logic straightforward analysis. They accept the usual assumptions 
the analysis reinforced concrete beams and derive the standard formulas 
introducing such generalizations will permit them study the influence, 
the cross section, thrust well bending moment. Their thesis 
that the thrust must considered for the design the cross section resist 
shear and bond stresses properly. 

Perhaps without intention they have chosen example that about the 
strongest reasonable case for supporting that thesis. beam having 
prismatic ratio (length over depth) less than three contrasted common 
ratio building design from ten twenty. Moreover, this deep beam must 
resist great thrust produced uniform lateral pressure the earth over the 
22-ft height the structure. Figs. and establish certain major dis- 
crepancies between the theoretical bond and diagonal tension stresses (not shear 
stress labeled Fig. for the structure studied. discrepancies, 
probably little consequence, can anticipated for building frames, but this 
criticism the paper, which gives only casual attention building 
frames and appears intended study heavy foundation structures. 
excellent brief presentation interesting thesis. 

The writer has ardently supported mathematical analysis about the only 
device available the designer for the production balanced structure; but 
the designer must stop and analyze the real meaning such calculated stresses 


reach logical decision where the refinement proposed the 


authors valuable and where merely theoretical 


paper Stanley Benscoter and Samuel Logan was published March, 1944, 
Discussion this paper has appeared Proceedings, follows: May, 1944, Shanley, 
Aleck, Dean Peabody, Jr., and William Wilbur. 
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few assumptions used the authors may reviewed profitably because 
similar assumptions are always necessary for the analysis any structure. 
The analysis moments and thrusts follows from assumptions the loads 
shown vertical uniform pressures the top and bottom the box (Fig. 4). 
Such pressure lines are never uniform but are actually curved. Likewise, the 
side pressures vary; they increase toward the bottom and are also likely 
unsymmetrical with respect the vertical axis. These factors seldom can 
brought into design because they are not determined theory but local 
conditions. 

Based upon conventional load assumptions, designer computes set 
moments using stiffness factors, all likelihood computed from J-values 

3 

with the transformed section (see Eqs. 10), what transformed section shall 
use? Evidently, since each member frame usually has one more points 
contraflexure, the transformed section could either the usual one 
upside down section. reality, over large part the length the beam 
the effective section certainly section uncracked concrete, since crack 
covers only fraction aninch. Actually, the beam has constantly varying 
moment inertia, the degree variation being unknown. desperation, 
engineers can excused for abandoning the idea using the transformed 
section and adopting for the moment inertia determining stiffness 
factors; but doing they must not forget that this crude assumption has 
had its inevitable influence computations for bending moments. 

Then there the problem the corner action which cannot solved 
common methods analysis. Perhaps the designer takes the theoretical span 
the beam the clear span; perhaps chooses the span center-to-center 
columns; perhaps with corner fillets uses less than the clear spacing be- 
tween columns the span. Whichever chooses, the result expressed 
bending moment diagram only approximation—a result judgment 
well mathematical analysis. 

Such crudities analysis are even more obvious when the member contains 
both bending moment and thrust. The fit structure within confining 
material may shift line thrust small distance that would 
consequence represented the vertical reaction under horizontal beam; 
but, represented the location longitudinal thrust acting vertical 
member, may have enormous influence upon the stress the member 
because shift even in. the eccentricity longitudinal thrust could 
disastrous. 

Finally, Figs. and reveal the theoretical fact that section moduli, 
stresses, bending moments, bond and shearing stresses all change sharply, 
whereas common sense dictates that such changes are actually gradual because 
stress builds slowly and not suddenly. This inaccuracy the result the 
designer’s inability treat the structure mathematically series small 
particles. must adopt gross procedures such the treatment 
beam cross section single element. 
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necessary analyze the suggestions offered this paper with these 
facts mind order decide just where the proposed equations should 


used. The specification formula, 120 240, not statement 


about any stress that actually exists the structure. just rule for con- 
trolling size cross section, more justified experiment and observation 
its satisfaction practical cases than its theoretical derivation. other 
words, whether the calculated shearing stress high low immaterial 
the specification writer takes care that fact his choice allowable 
working stress shear. 

The example offered the authors seems exaggerated case, and the 
effect thrust bond and shear beams columns common length- 
depth ratio may not usually sufficient importance consider. most 
the inaccuracies that have been mentioned this discussion exist the 
analysis given reinforced concrete structure, unlikely that the designer 
can hope know the stress given point any closer than 25%. That 
being so, can obtain little benefit from the theory presented the paper 
unless the influence the thrust bond diagonal tension proved 
about this magnitude. 

worth noting Fig. that the theoretical shear reaches 180 per 
in. cross section where the shear formula gives 120 per in. and the 
authors’ theory shows the critical shear causing diagonal tension only 
per in. Presumably, shear reinforced concrete not very serious 
cannot produce diagonal tension. However, well recognize that 
the formulas the paper, which show very low shear stresses, refer only 
shears the diagonal tension areas the beam; elsewhere, larger shears than 


those given may exist. Even though such increased shears may 


not serious, the authors have called attention change stress pattern 
that designers have probably been overlooking. serious fault over- 
look possible stress conditions. They have habit suddenly becoming 
critical. For example, the authors have not advised designers look for the 
possibility that negative thrust tension) can upset their computations 
producing more serious diagonal tension and bond condition than those 


represented the formula Nevertheless, the theory presented 


the authors can applied the case negative thrust. 


reinforced concrete beams may designed more conveniently the section 
modulus concept than the methods common use. This may true for 
beams plain bending, where the section moduli and are referred the 
centroidal axis which, that case, coincides with the neutral axis. When 


Civ. Engr., Board Transportation, New York, 
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thrust also acting, however, the fact that the two axes longer coincide 
leads complications. 

Formulas published the 1921, although similar Eqs. 12, are 
applicable both types problems. bending moment and thrust P’, 
applied distance from the center line tensile steel, may replaced, 
according well-known principle, eccentric force P’, acting 


distance from the tensile steel; that is: 


When thrust acting, the term drops out ©), which case 
Pe=M. 

Treating the moment, the equilibrium equations for rotation about 
the respective centers flange strain yield, for linear stress distribution, the 
following formulas which denotes the critical stresses: 


Designating the fractions Eqs. and respectively: 
and 


represent the section moduli the cross section about the neutral axis. 
The average shearing and bond stresses are found be: 


and 


There are numerical from which and can determined 
veniently, using and for arguments. When the tables are not available, 


the problem solved first determining the value which the only 
unknown Eqs. easily done means Eq. 36, which resolved 


Stresses Reinforced Concrete Arches and Other Structures,” Anders Bull, 
(Detroit, Mich.), December, 1921, 242. 
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into the following third-degree equation: 


For plain bending and Eq. reduced to: 


Evidently the hypothetical example chosen the authors has not been 
correctly treated, since the external and internal moments not balance. 
For cross section in. from the face the support (Fig. 3), the moment 
due roof load 791 in.-kips, and the moment due thrust 
16.2 713 in.-kips—a total external moment 1,504 in.-kips. Since the 
center steel 4.4 in. from the surface, the stresses indicated the diagram 
could develop only the following internal moments: 


Tension steel.......1.203 417 18.4 92,310 

Total internal moment (inch-pounds)........... 721,580 


—an amount which suffices for balancing the moment due thrust only. 
Furthermore, with thrust kips applied in. above the bottom the 
slab, there could point inflection Fig. the roof load moment 


cal meaning, the thrust remains the same height throughout the span. 

For comparison, Fig. shows the stresses obtained from Eqs. The 
thrust applied 6.7 in. from the bottom the slab, being the resultant 
3.5 kips acting the center the cross section and 7.5 
kips acting the bottom. 

novel feature the paper the distribution the shearing stresses 
across transverse section. The findings are interesting from theoretical 
point view, but their practical value may questioned, the average 
shearing stress much more easily determined. The assumptions which 
reinforced concrete design based are far from consistent regard tensile 
stresses concrete. Such stresses are considered when diagonal but are dis- 
regarded when longitudinal, although the latter stresses obviously contribute 
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the bending strength. The aim simplify the computations. For the 
same reason any highly refined treatment would seem somewhat futile. After 
all, the allowable stresses governing the design are determined the cut-and- 
try method. The possibility future excavation alongside subway structure 
like the one Fig. must also considered may cause the thrust 
removed temporarily, eliminating its relieving effect the shearing stresses. 


16852 per In. Center Line Steel v=107 per In. 


n 


Kips 
n=10 


per 


Axis_ 


Axis 


The significance the bulging shape the shear distribution curve 
easily seen investigating another cross section, located little closer the 
center span. The new stress distribution line (inclined dashed line Fig. 
then found intersect the old line the centroidal axis, showing that the 
compressive stresses below the axis are reduced whereas those above are 
increased. These changes conform the shape the curve, whose first deri- 


vative 


proportional the change compressive stress any given height 
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AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 1852 


DISCUSSIONS 


ANALYSIS RIGID FRAMES SUPERPOSITION 


Discussion 


FISCHER, LEON BLOG, ARTHUR MCGEE, 
AND JAROSLAV POLIVKA 


presented the author certainly another valuable method for the analysis 
rigid frames. The illustrative examples presented the paper show that 


2.00 0.00 
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0.08 
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40% 40% 


27.3% 
0.00ER 


2.00 4.00 1.09 0.82 
-0.19 0.58 0.09 
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paper David Wilson was published February, 1944, Proceedings. Discussion 
Goldberg. 
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the method general and can used for complex frames. 1938 the 
writer found the method superposition very practicable for calculating 
moments structures several bays wide and several stories high. 

Example Three-Column method moment 
distribution presented Professor Wilson general. some cases con- 
siderable time can saved taking advantage permissible short cuts. 
For example, the moment distribution Fig. 3(6) can shortened considerably 
first balancing joint and then carrying over one half the balancing 
moment the adjacent joints, before the adjacent joints are balanced. This 
method correct theoretically, and the writer uses all his moment dis- 
tribution analyses. Fig. shows this method balancing moments beginning 
joint After carrying over one half the balancing moment the ad- 
jacent joints, the next step balance joint carry over one half 
the balancing moment joints adjacent joint and then balance joint 


69.2% 


9.00ER, 
0.34 0.76 
0.02 
0.03 


53% 


53% 


1.52 0.06 

0.00ER, 

0.76 

0.34 0.02 


15.—Two-Srory Bent Loaps 


comparison with the formal procedure the paper the amount labor 
found reduced considerably. 

Example Bent with Horizontal symmetrical 
structure, such Fig. the K-factors for the horizontal members can 


parentheses, thus: (20), refer corresponding items the Bibliography (see Appendix 
the paper), and the end discussion this issue. 
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modified slightly assuming that there point inflection, hinge, 
the center, which makes the K-factor equal Thus, the 
K-factors for members and are equal 1.5 1.5 2.25, and 
balancing moment need carried over the adjacent joints for these members. 
Fig. 15(a) gives the moment distribution for the final moments shown 
Fig. and Fig. shows the moment distribution for the final moments 
shown Fig. 9(c). With this slight modification, considerable labor can 
saved the method theoretically correct. 

Example Trapezoidal Bent.—Fig. 16(a) represents single-story trape- 
zoidal bent similar Fig. but with bottom member and horizontal 
load 1,000 applied joint (compare analysis this bent Amirikian, 


Since for member equal (Fig. 16), the fixed-end moment Map 


(Fig. The first joint balanced joint and one half the balanced 
moments are carried the adjacent joints. Next balance joint and carry 
one half the balanced moments its adjacent joints, balance joint and 
carry one half the balanced moments its adjacent joints, and balance 
joint carrying one half the balanced moments its adjacent joints. 
Repeat the process the order joints and The final moments 
are shown brackets Fig. thus accounting for steps and 
step the author’s procedure (see Fig. 16(c)): 


also, Fig. 16(d), 
and 


The values and Eqs. are shown Figs. 16(e) and 16(f), re- 
Fig. 16(e), 


29.51 


and, Fig. 16(f), 
40.11 
step substituting the value from step Fig. gives the 
following final moments foot-pounds: 2,006 Map; 2,218 


5A. Truss with Inclined problem, 
course, differs somewhat from Example the paper (compare analysis 
Mr. Amirikian (1b)). Fig. represents five-bay Vierendeel truss with 
load each lower panel point. The value considered unity 


Correction for Transactions: indicate final sketch 
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throughout, and the values J/L are shown. Figs. and 17(c) 
show the left portion the truss the distorted condition with the R-values 
expressed terms and using Eq. 21. For this particular symmetrical 


0.86 

0.41 


1 2 
(c) 

(e) 


truss assume hinge the center the members AA’, BB’, and CC’ and 
multiply each K-value 1.5 that, with this slight modification, balancing 
moment required carried over the adjacent joint for these members. 

The fixed-end moments for the various members are shown Fig. 
18(a). For the moment distribution method the sequence is: Joints 


. 


| 
1.14 0.43 3.86 
7.72 
0.64 
2.28 1.26 
0.39 0.31 
031 
4.58 2.29 
167 
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and For member the carry-over balanced moment from joint 
joint one half the balance moment joint with the sign changed, 
because then necessary only apply the moment distribution method 
the left part the truss when the load and truss are symmetrical. The final 
moments terms are shown brackets. 


K=2 


FI 


400 200 200 200 200 400 


TRUSS WITH INCLINED CHORDS 


Fig. 18(b) gives the fixed-end moments terms and the sequence 
balancing is: Joints and thus accounting for steps and 
Fig. 18(c), for step 


and, Fig. 18(d), 


(41a) 
and, Fig. 18(d), from summation moments 


The values and Eqs. are shown acting Figs. 18(e), 18(f), 
and 18(h). Figs. 18(e) and 18(f), respectively, with and 
Figs. 18(g) and 18(h), 
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and 
Substituting Eqs. 40, 
and 
13.46 58.24 4,800................. (43b) 


Solving Eqs. determinants, 78.12 and 64.36. 
The numerical values the moments Figs. and 18(6) may now 
calculated (step 5), and the following moments foot-pounds the ends 

The foregoing intended lend emphasis the writer’s sincere opinion 
that the method superposition presented Professor Wilson 
unusually valuable tool for the solution rigid frames. 


importance using the method superposition 
when analyzing rigid frames emphasized this paper; but the paper has 
additional value. the several examples, the author has explained various 
simplifications and clarifications previously unknown the writer. Although 
superposition was not required for solving Example there considerable 
merit the technique its solution. separating the R-term Eqs. 
and writing separately Eqs. and then expressing the R-term 
other members terms that one member, considerable simplification 
the statement the values for the end moments occurs. Through stating 
necessary evaluate and quantitatively. Values and radians 
are generally low order and lead errors the final moment values unless 
carried throughout the solution with many significant figures. Another 
simplification that the value occurs both sides Eqs. that, the 
author’s method, need not considered separately. Only the product 
must evaluated and substituted the values the end moments which 
mum error the moment values sought. The technique just described 
better than any known the writer for the solution rigid frames either 
the slope deflection the moment distribution methods. 

Superposition method solving, steps, certain problems which 
applicable. When solving problems dealing with rigid frames involving 
asymmetry any kind, one-step solutions moment distribution are not 
possible. Two steps are required: Moment distribution evaluate mo- 
ments due pure rotation the joints; and (b) moment distribution obtain 
moments due movements occurring transversely the axis each member 
that can deflected. The moments evaluated each step are added alge- 
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braically each joint derive the required moment. Step may solved 
any good system moment distribution. Step often solved the 
unbalanced shear method (2). Hardy Cross, Am. Soc. E., does not 
step anything, but unbalanced shears are actually the tools required 
when correcting for sidesway his method. 

Professor Wilson does not use the unbalanced shear method step 
solving Example Fig. 3(b), distributes fixed-end moments which are 
distributed moments. The moments thus found are added algebraically 
those Fig. 3(a) get the working moments. This consummates very 
direct solution Example 

analyzing two stories, necessary use set R-values: appertains 
the lower story and the upper tier. After that variation from Ex- 
ample the technique solution parallels that for Example 

For the solution Example the author evaluates moments terms 
EA. Here, inspection will not yield R-values for the vertical members 
terms the R-value selected vertical member illustrated Example 
this example was evident that the R-values for vertical members joined 
the same horizontal member are inversely proportional the length the 
vertical members. Fig. 5(a), the joint between points and pre- 
cludes evaluation deflection angles inspection. Otherwise the method 
solution closely follows that for Examples and 

For solving the Vierendeel truss Example Professor Wilson able 
profit expressing moments terms Otherwise, the solution self- 
explanatory. 

Equally valuable the exposition the relation between the values 
for the members quadrangle. The knowledge this relation permits very 
direct attack upon the type problem considered the ensuing comment. 
The meaning Fig. may made somewhat clearer drawing 
horizontal line through toward the right meet perpendicular drawn 
downward from C’. 


For the solution the frame Example the R-relation Eqs. and 
was used. Fig. not quadrangle, but the fixity the legs and 
and are applicable the solution Example The distributed moments 
are determined using manner closely resembling the solution 
Example 

Example Fig. solved easily using the quadrangle concept. The 
lower story has legs which are fixed, and fictitious member may 
assumed. The second story full quadrangle. Otherwise the solution 
closely follows that Example 
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Fig. 10(6) Example the author uses the quadrangle. The dashed 
line corresponds the basic quadrangle. Moments are all ex- 
pressed terms and and thereafter the equations statics 
dominate the solution. 

Finally, the writer believes that Professor Wilson has made distinct con- 
tribution the analysis rigid frames homogeneous material. Should 
any the members variable shape, the resultant variation end moment 
inertias would preclude the simplification which results Eqs. using 
the concept relative K-values. However, each Eq. would merely contain 
the K-value for that member and would lead some additional numerical 
work. the frame not homogeneous material, modular ratios such 
those for reinforced concrete may used establish all the equations 
terms the selected member. This also would involve some extra 
numerical work, but the extra work would common other known methods. 


value for the analysis rigid frames. The examples the paper illustrate 
the method procedure workable manner; but, since Examples 
and can solved shorter and easier methods, they not demon- 
strate the value the method. the method justice and show its 
great value solving difficult problems the writer wishes show its applica- 
tion two-bay, saw-toothed structure with tapering members. 


19.—Two-Bay, WITH TAPERING MEMBERS 


Fig. shows such frame loaded unsymmetrically with vertical and hori- 
zontal loads. The tapering members have been given the same moments 
inertia the ends that one set beam constants will suffice for all members. 
These constants are taken from the tables compiled Walter Ruppel, Assoc. 
Am. Soc. application these constants the several 
members given Fig. (22). 
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The notation used this discussion varies slightly from the author’s but 
should cause confusion. The symbol used herein for the linear dis- 
placement one end member, and represents angle equal 
Also the writer considers end moment the moment 
member rather than resisting moment acting the member. This, 


u=0.250 v=0.427 
(a) p=0.292 
(c) 


(d) 


Relative Relative A=1.71M Relative Relative 


course, reverses the signs the end moments from those used Professor 
Wilson. 

Table gives the effect applying positive (clockwise) moment each 
joint the frame, without sidesway. The table may computed moment 


TABLE Moments APPLIED JOINTS 
(Not Corrected for Sidesway) 


+1,000 
moment applied 


distribution. The writer, however, prefers the traverse method (22). Fora 
moment joint the traverse shown Fig. 21. The moments the 
diagram are multiplied obtain the results for joint Table 
corresponding traverse must run for joint The intermediate joints, 
and are easily calculated using data taken from the two 
already run. Table will apply any combination fixed-end moments 


that may act the frame. 
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When the frame permitted deflect indicated Fig. 22(a), the 
relative displacements the members AB, BC, and are shown Table 
the triangle abc, side drawn normal member AB, side normal 
member BC, and, likewise, side normal member CD. These displace- 


o 


ments are —0.8d, and respectively. The equations for fixed-end 


inertia members, are follows: 


Eld 
and 


Table shows the fixed-end moments for members AB, BC, CD, and DH, 
Fig. 19. Member like member BC, member like member CD, and 
member like member BA. 

There are three conditions deformation for this frame. The first con- 
dition shown Fig. 22(a). The coefficients calculated Table 
have been multiplied 1,000 arbitrarily place the decimal point more 

1,000 
AB, BC, and are placed the diagram their respective locations and 
the moments are balanced the use Table Thus, joint 204.12 
128.75 75.37, clockwise moment which must balanced throughout 


convenient location. The fixed-end moment coefficients for members 


the frame. Therefore, times the values Table for joint will 


added algebraically the respective end moments the members. Like- 
wise, joint —75.40 31.81 —43.59, the unbalanced moment, and 


lis- 
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128.75 
10.85 
1.20 


149.18 


(a) 
0.44 
0.05 
270 
84.79 


049 
+123.69 
14.08 +14.08 061 
+1232 


+15:72 


184 


Discussions 
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AB, BC, CD, DH, Fic. 


TaBLE 2.—RELATIVE DISPLACEMENTS AND MEMBERS 

: 


ie 


Member Displacements Moments 


| | ; 
| 
6.04 
9.49 
0.20 
6.04 
28.39 
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—43.59 
1,000 
Table will spread this unbalanced moment throughout the frame. The same 
procedure followed for the unbalanced +54.32 moment joint The 
end moments for each member are now added and the first condition defor- 

mation finished. 


200 
500 


times the moments shown for +1,000 ft-lb moment joint 


(c) 


1 
54.84d3 26.50 
(i) 
Fig. 


Fig. 22(b) shows the second condition deformation, and the coefficients 


computed before. The third condition solved Fig. 


Fig. and the remainder this discussion, the expressions 


have been contracted and respectively. 


82.25d 127.39d 


Try 


rr 


rr 
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From Fig. 23(a) evident that 


Similarly, Fig. 23(c), 

and, Fig. 23(d), 

Substituting Eqs. 46c, and 46d 46a gives 


Fig. 23(f), taking moments about point 


Fig. 23(g), taking moments about point and substituting the value 
from 


Substituting Eqs. 46d and 47c Eq. 47a gives 


From Fig. 23(h), 


Fig. taking moments about point 


Fig. 23(j), taking moments about point and substituting the value 
from Eq. gives the equation, 


Substituting Eqs. 46b and 48c Eq. 48a gives 


Solving Eqs. 46e, 47d, and 48d simultaneously gives —15.512, 

The moments Figs. 22(a), 22(b), and 22(c) multiplied their respective 
d-values are shown Table final values the moments the ends 
all members are given the bottom line Table 


same 
The 
efor- 
9.494, ] 
| - 
q 
q 
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Additional bays can added the frame Fig. without increasing 
the problem unduly. Each additional bay would add one more condition 
deformation. other words, such three-bay structure would require four 
conditions deformation; four-bay structure would call for five; 


MEMBER 
Factors 


The solution this type problem superposition may rather long, 
but the individual steps are logical and easily followed and the end difficult 
problem has been solved. 


JAROSLAV Am. Soc. principle superposition 
generally recommended and used cases which the direct analysis either 
involves the solution greater number simultaneous equations requires 
relatively tedious determination certain constants the structure. For 
example, using the slope-deflection method, the examples presented the 
paper require the solution the following number simultaneous equations: 


Example No. Number equations Example No. Number equations 


Examples and there are only four equations, course, since the frames 
are symmetrical. 

most cases several load systems must investigated, and for this reason 
preference should given method analysis which certain constant 
values characteristics the given structure, common and applicable any 
type loading, can determined. The method analysis presented the 
author possesses this advantage, since the computation the internal forces 
(bending moments, shears, and thrusts) produced lateral displacements 
the joints treated independently transverse load. 

Since accuracy greater than usually not required for design 
poses (except for checking the more complicated statically indeterminate struc- 
tures), graphical combined methods are useful for types framed structures 
which the method superposition applicable. One the most promising 
methods that fixed (characteristic) points, which has not yet been fully 
utilized the United States. The solution general and can applied 


Research Associate, Univ. California, Cons. Engr., Kaiser Co., Berkeley, Calif. 
Received the Secretary May 1944. 
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hinged and elastically yielding supports. The latter condition (yielding sup- 
ports) often more serious engineering practice than either the extreme 
conditions fully restrained hinged supports. 

The designer can adapt graphical, semigraphical, algebraic processes 
his purpose. Graphical devices are quicker; they are easier check; and they 
are accurate the precision the draftsman who prepares them. 

graphical solution points, applied Example demonstrated 
reference Fig. 24. The determination the fixed (characteristic) points, 
and each structural member shown Fig. 24(a) using the values 
the elastic weights, and the angular displacements, the 
joints. The distances the characteristic points, and can checked 
calculated the following relationships: (full restraint point A); 


The angular displacement has the significance the elastic weight sup- 


port the member point and equals because restraint point 


The successive angular displacements can expressed terms 
distances between the opposite fixed points and the supports; for example, 

Fig. shows the bending moments the bent, subjected central 
load member BC, without considering the sidesway, found the fixed 
points and the crosslines for the loaded span. The crosslines are determined 
values and Ke, derived from fixed-end moments from formulas for 
each type transverse loading. 

take the lateral displacement the column into account, the shears 
and are calculated from the moment diagram, and the sum, 
opposite sign, applied the temporary support, point The effect 
this temporary reaction obtained allowing arbitrary translations 
the columns and calculating the shear each column. The shear 
the total elastic weight column AB, restrained point and elastically 
supported point and the semiaxis the ellipse elasticity this 
column. The elastic centroid and are found graphically shown Fig. 
24(c). Simple relationships permit the computation moments points 


6 Ap 


Mpa = ~ — (49) 
and 


The moment then distributed over all remaining structural members. 
Similarly, the moments due and which are produced 


. 
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Ug=3.33! Third Points v2=48 aime 85! 3 4= 2.86! 


20! L4=20! 4 


(a) CHARACTERISTIC POINTS 


284 M) 
(b) MOMENT DIAGRAM 


=104 


Moments Due 
—-— Moments Due 
Total Moments 


(c) DETERMINATION ELASTIC CENTROID 


PL> 
Mce= —1268 
— 3 
= 
H 1 H 2 
ie 
4 
| 
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are calculated. Finally, all moments produced equal lateral dis- 
placement the columns are summed and the corresponding shear along 
member BCE calculated. 


The moment distributed the other joints using the final carry-over 
factors which are determined fixed points and angular displacements (at the 


joint C). The carry-over factor from point point Cac 


etc. Similarly, the values and Mpc (for column DC) 


and (for column FE) are obtained and carried over the other 


joints (Fig. 24(c)). For the purpose superposition, only the ratio between 
the temporary reaction point (Fig. and the total shear along member 
ECB (or the total thrust along this member) are considered. Therefore, the 
common denominator 710 can eliminated (Fig. 24(c)). The temporary re- 


action point (Fig. and the total shear produced 


The additional moments produced sidesway are calculated multiply- 
ing the values Fig. 24(c) the ratio 0.63934. 

Professor Wilson uses the moment distribution method (except Example 
which the slope-deflection method used) for the determination the 
superimposing effects. The relative accuracy depends the number cycles 
and usually completely sufficient for rational design shown the values 
obtained the author compared with more accurate results the writer’s 
computations the method fixed (characteristic) points: 


hence, the ratio sought 


Example Moments 


Mcp 933 (—935.76); Mpc 557 (—558.67); 
1,612 (—1,612.93); Mpc 799 (—808.81); 
2,203 (2,212.61); 1,109 (—1,112.61); and 

8.71 (8.79); 14.58 (14.65); 2.70 (2.82); 
Mep 4.06 (3.92); Mpc 6.30 (6.31); 7.75 
12.15 (12.02); 6.76 (—6.74); 

2.66 (—2.64); Mra 9.49 (—9.38); and 


| 
112 
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Example Moments 

228 (230.93); 428 (429.18); 640 (642.25); and 

569 (568.95). 

Mcp 2,250 (2,256.11); Mpc 550 (539.99); and 
3,800 (3,831.97). 


the foregoing, the writer’s values are parentheses. The results are ac- 
curate, the differences being less than 4%. Example differences were 
found. 


(1) “Analysis Rigid Amirikian, Government Printing 
Office, Washington, (a) (6) Fig. 104, 145. 
(2) “Analysis Continuous Frames Distributing Fixed-End Moments,” 
Hardy Cross, Transactions, Am. Soc. E., Vol. (1932), pp. 1-156. 
Fischer, Concrete, November, 1938, 
(21) Transactions, Am. Soc. E., Vol. (1927), 170. 
(22) “Relative Flexure Factors for Analyzing Continuous Structures,” 
Ralph Stewart, ibid., Vol. 104 (1939), 521. 
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SIMULTANEOUS CONSOLIDATION 
CONTIGUOUS LAYERS UNLIKE 
COMPRESSIBLE SOILS 


Discussion 


BARBER, AND JACOB FELD 


1 

Epwarp timely analysis the consolidation 

layers foundation soils having different consolidation charac- 

teristics offered this paper. The mathematically correct approach pro- 
vides opportunity evaluate the results furnished the less laborious but 

approximate methods. approximation that provides reasonably simple 

solution for the consolidation any number contiguous layers and takes 

into account influencing factor not considered the author’s method 
average coefficients—namely, the relative position the consolidating soil 
layers and the drainage faces—is outlined this discussion. This approxima- 
tion produces results that agree more closely with the exact theoretical analysis 
those the method average coefficients. 

considering the average gradient the incompressible layer, shown 
Fig. 1(a), required carry away the water squeezed out the consoli- 
dating layer, any degree consolidation, approximate formula 
obtained for the curves Fig. follows: 


have the coefficient permeability the drainage layer, large enough 
provide unrestricted The results furnished Eq. are good 
for the combination large values and values 


paper Hamilton Gray, Jun. Am. Soc. E., was published February, 1944, 
Proceedings. 


Associate Highway Engr., Public Roads Administration, Washington, 
Received the Secretary May 10, 1944. 
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The dimensionless coefficients, and used Professor Gray’s 
formulas may simplified expressing them terms the strain per unit 


can obtained directly from the dial readings consolidation test. Actual 
thicknesses and permeabilities may used instead their reduced values and 
the void ratio becomes unnecessary. 


Thus: 

ky 8) 

(38c) 

2 


Reference now made the comparison between the results furnished 
the method average coefficients and those the exact analysis. 
appears that the consolidation 


TABLE several layers should averaged the 
LAYERS ILLUSTRATIVE basis the settlements rather than the 
EXAMPLE thicknesses, shown Eq. 35. such 


case, the agreement between the exact 


Total agreement the two methods discusse 


the author may not always close 

that indicated Fig.8. illustrate, 

assume two contiguous consolidating layers with one drainage face and prop- 
erties given Table 

Using the method average soil coefficients outlined Eq. and 


4.7.5 
(a) CASE (b) CASE 


Free Drainage 


Free Drainage 


Drainage Drainage 


Using Professor Gray’s exact method, there are two cases (see Fig. 10): 
Case impervious face next layer and case impervious face next 


Consolidation, 


of 


Degree 


: 
4 
4 
q 
=0.4 k=10 
| 


June, 1944 BARBER CONSOLIDATION SOILS 995 


layer II. The constants are: 


Constant Case Case 
0.5 


The time-consolidation curves resulting from these two methods calcu- 
lation are shown Fig. 11. Also, Fig. are shown the curves obtained 
means the approximation introduced the writer (which modifies the 


Case Exact 


0.2 


Consolidation, 


0.6 


Degree 


1.0 
0.5 


method average coefficients take into account the effect the relative 
position the boundaries and thus agrees more closely with the exact method). 
The effect this approximation illustrated the example for three 


contiguous layers outlined Table The time computed the method 
average coefficients. Thus, using the relationship, 


| 
« 
this, 
Case Approximate 
Time, 
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the elapsed time expressed 


which becomes for However, the actual layers are replaced 
their H/k (Col. Table and, instead using the 


TABLE THREE LAYERS 


GIVEN 
SETTLEMENT 
Layer 41) For 
Final Thick- Equiva- Layer 
ness cient lent ent part (Total Hp) 
No. ment layer permea- thick- 
(1) (2) (3) (4) (6) (7) (8) (9) (11) 


This for drainage above and below. For impervious lower boundary make this value 0.50. 


Degree Consolidation 


0.2 0.4 0.6 1.0 Free Drainage 


Part Thickness, 

> 


1.0 0.8 0.6 0.4 0.2 Free Drainage 
Hydrostatic Excess Pressure, 


Fig. ANALYSIS TIME-CONSOLIDATION COMBINED LAYERS 


same percentage load transfer each layer, the part the load transfer 
each part the total used shown Fig. and Table 


Technical Memorandum No. Bureau Reclamation, Dept. the Interior, 1939, Fig. 34. 


> 
di 
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which values AU, are taken from Table the difference between the 


TABLE TIME-CONSOLIDATION COMBINED LAYERS 
(Degree Load Transfer Layer, Up) 


Part total 
thickness 


0.30 


0.40 


0.50 


0.60 0.70 0.80 0.90 


0.10 0.20 
0.00 0.000 0.000 
0.05 0.028 0.037 
0.10 0.042 0.062 
0.15 0.046 0.081 
0.20 0.048 0.088 
0.25 0.049 0.095 
0.30 0.049 0.096 
0.35 0.050 0.098 
0.40 0.050 0.100 
0.45 0.050 0.100 
0.50 0.050 0.100 
0.55 0.050 0.100 
0.60 0.050 0.100 
0.65 0.050 0.102 
0.70 0.051 0.104 
0.75 0.051 0.105 
0.80 0.052 0.112 
0.85 0.054 0.119 
0.90 0.058 0.138 
0.95 0.072 0.163 
1.00 0.100 0.200 


DEGREES CONSOLIDATION 


f@ 
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
0 0.049 0.100 0.154 0.217 0.29 0.38 0.50 0.66 0.95 yi 
0.2 0.027 0.073 0.126 0.186 0.26 0.35 0.46 0.63 0.92 
0.4 0.016 0.056 0.106 0.164 0.24 0.33 0.44 0.60 0.90 
0.6 0.012 0.042 0.092 0.148 0.22 0.31 0.42 0.58 0.88 
0.8 0.010 0.036 0.079 0.134 0.20 0.29 0.41 0.57 0.86 
1.0 0.008 0.031 0.071 0.126 0.20 0.29 0.40 0.56 0.85 
1.5 0.006 0.024 0.058 0.107 0.17 0.26 0.38 0.54 0.83 
2.0 0.005 0.019 0.050 0.095 0.16 0.24 0.36 0.52 0.81 
3.0 0.004 0.016 0.041 0.082 0.14 0.22 0.34 0.50 0.79 
5.0 0.003 0.013 0.034 0.069 0.12 0.20 0.32 0.48 0.77 
10.0 0.003 0.011 0.028 0.060 0.11 0.18 0.30 0.46 0.75 
« 0.002 0.009 0.024 0.048 0.09 0.16 0.28 0.44 0.73 
aa 0.002 0.008 0.018 0.031 0.05 0.07 0.10 0.14 0.21 


Pressure drainage face divided pressure impervious face Values time factor for two 


drainage faces. 


U,-values the limits the layer. Thus, for layer (see Col. 10, Table 5), 
when 0.60 the value AU, (for 0.30 0.80) equal 0.440 


0.23 


0.50 


Pp H, ( ) 
q 
0.000 0.000 0.000 0.000 0.000 0.000 0.000 
0.044 0.046 0.047 0.048 0.048 0.049 0.049 
0.103 0.113 0.123 0.129 0.134 0.140 0.145 
0.121 0.138 0.149 0.160 0.171 0.180 0.191 
0.147 0.188 0.224 0.258 0.294 0.328 0.365 
0.149 0.194 0.236 0.278 0.320 0.360 0.405 
0.150 0.200 0.250 0.300 0.350 0.400 0.450 
0.179 0.262 0.351 0.440 0.529 0.620 0.709 
0.197 0.287 0.377 0.471 0.566 0.660 0.755 
0.256 0.354 0.453 0.552 0.652 0.751 0.851 
0.300 0.400 0.500 0.600 0.700 0.800 0.900 
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The writer has assumed that the increase pressure applied the con- 
solidating layers uniform throughout their total thickness. include the 
effect trapezoidal distribution the initial excess pressure, values 
may taken from Table which based the following interpolation 


which constant and the subscripts are values 


Am. Soc. clearly indicating how the limita- 
tion actual field conditions affects the theoretical solutions, the author has 
provided valuable method for determining the range values, not only 
magnitude but also rate occurrence. This procedure will not mislead the 
engineer who looks such tools for guidance the determination expected 
consolidation soil fills and under foundations. All theoretical solutions are 
based assumed end and surface conditions. Unless such ideal conditions 
occur, the formulas can only approximate true values and probably only give 
range within which results may expected. Only when prophecy 
consolidation, both extent and rate, made the basis theoretical 
conditions alone, correlation found between actual occurrence and ex- 
pected values, and one loses faith all quantitative soil studies. 

Sudden clear-cut changes soil layers are not often encountered nature 
—even rock formations. There usually conglomerate, mixed skin 
separating unlike soil layers. Such condition difficult evaluate 
theoretical solutions, and thus subjects this nature fall the class the 
practical art rather than that the theoretical science. 

The same problem occurs building road and airport surfaces sub- 
base materials. When layer pervious material placed directly clay 
subgrade, apparently the ideal condition contact soil layers occurs. 
Actually, experience indicates that the two materials mix and form interven- 
ing skin, which similar neither the parent layers. accelerated field 
demonstration this phenomenon (probably not new) occurred during the 
construction some twenty-five miles pavements extremely fine and 
impervious Seneca County (New York) clays 1942. examination all 
existing roads that area, and experience other projects built the pre- 
vious year, showed that the fine clay worked through the run-of-bank gravel 
sub-base, through penetration macadam, and even through bituminous mix 
surfacing. one year, under rather light traffic, the clay had come through 

total thickness in. gravel sub-base, in. crushed dolomite base, and 
penetration surfacing. 

However, inspections, after almost two years use, the roads built 
under the writer’s supervision showed vertical intrusion the clay into the 
road construction. these cases, fine bank sand containing leachings from 


Barber, Public Roads, March, 1937, 18, Eqs. 53, 54, and 55. 
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limestone deposits had been placed about in. deep the excavated roadbed. 
Over the sand in. in. run-of-bank gravel, in. crushed stone base, 
and in. tar mix surfacing were then placed. Similarly, some test roads 
built directly top 2-in. 3-in. weathered clay “topsoil” containing 
some grass roots were also satisfactory. Apparently the introduction 
fine sand screen, rolled into the clay and forming conglomerate skin, changes 
all the consolidation characteristics normally encountered the clay. 

consolidation forecast filled materials must also take into account the 
possible altering the exposed surface with the formation crust which will 
entirely change the expected free drainage surface conditions. the construc- 
tion large airport, grading operations included the removal glacial till 
hills and the filling valleys depths almost ft. The expected shrink- 
age and consolidation did not occur. When the paving was placed, about 
year later, was found necessary remove about the fills bring the 
surface subgrade elevations. This top crust was very hard and high 
density, but was floating almost liquid material which had not con- 
solidated. The result could prophesied the author’s methods had 
been assumed that there was practically surface evaporation. However, 
that assumption not made, erroneous result found from the use 
any consolidation formula. 

placing fills over completed subway and similar structures densely 
populated areas, the most important requirement rapid consolidation. 
Otherwise, there not sufficient time replace the subsurface utility struc- 
tures such manner that failures will not occur after the streets are repaved. 
Too little has been done experimenting with artificially treated backfill 
materials, when ideal porous filling materials are not available. The necessity 
for such treatment can determined from the methods developed the 
author, and proper assumptions the relative depths treated and untreated 
fills can also obtained. one personal experience, the use river silt 
taken from subaqueous tunnel operation was found unsatisfactory because 
the slow drainage characteristics. However, when the tunnel operation 
ran into some difficulties and appreciable amount hay was used for 
up” the heading the tunnel, the silt material became excellent 
backfill material and consolidated very rapidly. Once was consolidated 
there was little danger future volume changes because the restrictions 
subsurface structures, foundations, and impervious street paving. course, 
this idea ancient origin, dating back biblical, earlier, days when 
bricks were manufactured the addition straw clay. 

The author congratulated for the presentation workable formulas 
which can checked the field and which will give consistent and true results 
if, and only if, consistent and true field conditions are evaluated (as 
carefully explains) before using these formulas. 


Correction for Transactions: February, 1944, Proceedings, page 149, the 
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DISCUSSIONS 


ERECTION HYDRAULIC TURBINES 
AND GENERATORS 


Discussion 


PAUL GISIGER, AND FRANK MORGAN 


the main points observed the erection Kaplan and Francis wheel 
hydroelectric generating units. The writer would like add one recommenda- 
tion this competently presented paper—namely, that those who are actual 
charge erection make and keep permanent record all elevations, clear- 
ances, tolerances, etc., established the time erection. This record should 
checked and necessary supplemented immediately before the unit goes 


into regular operation. 

The writer has several times regretted the lack such record. 
fairly obvious that the levels and tolerances established during erection 
hydroelectric unit are not permanent. the best concrete subjected 
slight changes volume and shape due the factors shrinkage and 
creep, apart from temperature variations and elastic deformations. The dis- 
locations caused these factors are generally equal, not higher, order 
magnitude than the tolerances within which unit can erected. The result- 
ing changes alinement are most cases inconsequential and may un- 
noticed for years, but occasionally they are large enough require correction 
and then the question arises how much change has taken place since the 
time when the unit was erected. The horizontal alinement generator 
stator may have changed, the runner blades Kaplan turbine may have 
left trace some sector the circumference the discharge ring, bearing 
may show signs pressure concentrated one side, shaft may found 
slightly out plumb. each case several causes may responsible 
for the change. Without accurate record for reference difficult, not 
impossible, judge the nature and causes the trouble and apply suitable 
remedies. 


paper Andrew Komora, Am. Soc. E., was published March, 1944, 


Proceedings. 
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There are normally least three parties concerned with the erection 
hydroelectric unit—the purchaser, the turbine manufacturer, and the generator 
manufacturer. something does get out plumb alinement, natural 
for the equipment manufacturers question the structure which this equip- 
ment supported. All the interested parties then need complete record 
the condition the unit the beginning its useful life. 

From this viewpoint cannot recommended too strongly that com- 
plete record kept each unit, giving all key measurements and elevations 
such as, for instance, diameters the discharge ring least four directions, 
elevations marked points the speed ring, elevations the points sup- 
port the main bearing brackets, and elevations the stator base plates. 
addition these points parts the generating unit itself, will helpful 
maintain record observation points that permit following any slight 
change the structural supports. Such points may located, for instance, 
the two main diameters the lowest accessible points the pit liner and 
corresponding locations higher the concrete foundation the generator 
stators. The cost installing such observation points bench marks the 
time erection insignificant when compared the value the information 
they may reveal during the life the unit. 


longfelt need has been met the publica- 
tion this interesting paper. The writer wishes add his comments the 
hope that Mr. Komora will clarify few points. For example, under 
Placing and Erecting Turbine Embedded Parts,” the author states: 
level transit used [in setting equipment] depends the engineer’s 
preference.” not true that there choice since the transit not 
sufficiently accurate instrument for this purpose? The use high-grade 
level quite essential insure accurate results. 

The author mentions difficulties experienced with the I-beam type bearing 
bracket because improper shimming, etc., the side arms. This has been 
serious problem and much trouble has been experienced with guide bearings 
because distortion the bearing retainer resulted from increased deflection 
under load. This defect has been corrected later designs hinging the 
side arms. would interesting learn how the bearing clearance deter- 
mined, with the machine under load recommended the author. 

Mr. Komora mentions the use four plumb lines when plumbing the shaft. 
Why are four lines necessary? should possible obtain all the data with 
two plumb lines 90° apart—reducing the amount work considerably. 

The author states that customary make the coupling bolts mils 
oversize. What are the particular benefits this procedure? This practice 
requires considerable extra time and labor during subsequent dismantling work. 
Furthermore, such tight fitting bolts tend distort the faces the flanges 
and make difficult the accurate determination the stress, the bolts, 
during the process tightening. light tapping fit for coupling bolts 


Maintenance Engr., Pennsylvania Water and Power Co. and Safe Harbor Water Power Corp., 
Baltimore, Md. 
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considered adequate. the bolts are properly stressed when tightened, the 
friction between the coupling faces should sufficient transmit the torque. 
one particular case the bolt holes one flange shaft coupling were 
inadvertently reamed oversize. The unit, large one, was assembled with the 
coupling this condition, and inspection quite some time after the unit 
had been service, under load, showed that movement had taken place 
the coupling. 

Mr. Komora quite emphatic his condemnation the practice 
shipping the turbine shaft the generator manufacturer, there assembled 
with the generator shaft and the entire assembly checked for truth lathe. 
Such procedure advisable until the erection large units, which this 
course has not been followed, has proved satisfactory. least one case, 
such check revealed considerable inaccuracy which could corrected only 
remachining. Had the shafts reached the field, serious delay would have 
resulted. Shop assembly also makes unnecessary ream coupling bolt 
holes the field. 

Table the actual throw particular shaft given 0.0055 in. 
The author states (see heading, Erection Generator-Mechanical 
‘“The center this circle [the shaft mils south and 
mils west from the’ theoretical center line. The shaft was replumbed 
that the final position was mil west and zero mils north and 
mils extremely small for the throw shaft, and since the center the 
throw but mils away from the true center the position the shaft almost 
ideal—making any changes its location somewhat inconsistent with the 
author’s statement that much unnecessary time and effort has been spent 
alining units. this particular case the clearance the guide bearing which 
this shaft operates probably mils; therefore, the throw the shaft 
plus the error the position the center would still small enough 
within the bearing clearance. Units having shaft throws mils are operat- 
ing successfully although this is, means, desirable condition. 

Has Mr. Komora had experience plumbing the shafts and checking the 
throw units equipped with spring type thrust bearings? With this type 
thrust bearing the shaft may made absolutely plumb and centered the 
turbine guide bearing. soon the rotor turned, however, the shaft will 
longer plumb because unequal yielding the springs the thrust 
bearing. Thus, the final setting must compromise. 

The author states that period several weeks has been spent alining 
units having shrunk hubs. the surface the hub (more properly, the 
thrust block) not normal the shaft axis, how can the unit alined cor- 
rectly since amount adjustment the thrust bearing will eliminate 
this error? 
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matical treatment the problem aerodynamic stability suspension bridges 
contained this paper. appears, however, that some the data available 
the author relating the behavior the original bridge were error and 
that the author was not informed some the uniform conclusions from 
laboratory tests conducted since the failure the Tacoma Narrows Bridge. 
unfortunate that the pressure wartime has made impossible release 
the results certain very extensive investigations which, although not com- 
plete, clarify some important phases the problem. The following discussion 
refers corresponding parts the paper. 

Parts and 3.—One the most important errors the information available 
the public regarding the behavior the Tacoma Narrows Bridge lies the 
published log observed modes motion and associated wind 
most serious error arose the velocities obtained with anemometer for 
which calibration curve existed. Subsequent calibration this instrument 
revealed errors about 50% the lower velocity range. 

uniform procedure was followed obtaining these readings velocity. 
Many them were taken the top the toll house, which was 1,400 
distant from the nearest tower and where the influence the terrain was such, 
southerly winds, alter, completely, the velocity well the direction 
the wind. Many other readings were taken just beyond the first tower on. 
the main span and again offered reliable indication regarding the average 
wind the bridge whole. Furthermore, none the direction 
was obtained with any accuracy. reading was pure estimate, and 


paper Steinman was published November, 1943, Proceedings. Discussion 
this paper has appeared Proceedings, follows: February, 1944, Edward Adams Richardson; and 
March, 1944, Bernhard; and May, 1944, Robertson, and Gronquist. 

Prof., Civ. Eng., Univ. Washington, Seattle, Wash. 

Received the Secretary May 11, 1944. 

Failure the Tacoma Narrows Bridge,” report the Hon. John Carmody, Administrator 
the Federal Works Agency, Washington, C., March 28, 1941, board engineers consisting 
Ammann, Theodor von and Glenn Woodruff, Appendix VI, 19. 
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close correlation was maintained between observations velocity and direction 
and readings amplitude. short, these recordings should treated only 
very sketchy records the behavior the wind whose true characteristics must 
remain forever doubt. 

The net result the aforementioned discrepancies that all the velocity 
and direction readings the log are subject correction. general, the 
velocity error revealed the anemometer calibration tapers from approxi- 
mately 50% the lower range zero about miles impossible 
estimate the magnitude the errors due location the reading instru- 
ment, but can stated confidently that many the readings were very low 
for that reason. 

During the months following the failure the Tacoma Bridge extensive 
research was undertaken the Structural Research Laboratory the Uni- 
versity Washington Seattle under the joint sponsorship the Washington 
Toll Bridge Authority and the Public Roads Administration. program 
involved the building wind-tunnel laboratory capable testing 
scale model suspension bridge the proportions the one lost Tacoma 
Narrows. date (June, 1944) two complete models have been built, the first 
duplicate the original bridge and the second connection with the proposed 
new design. conjunction with these full model studies, three section models 
have been constructed part the investigation the proposed new design 
for bridge Tacoma Narrows. 

Paralleling this work the University Washington, program was 
sponsored the Washington Toll Bridge Authority the California Institute 
Technology under the direction Theodor von Am. E., 
using section models. 

The difference between section and full should clarified. The 
full model complete duplication the prototype, every significant detail, 
and includes towers, cables, suspenders, and suspended structures constructed 
appropriate model scale. The section model composed part the 
suspended structure mounted springs such manner that the various 
degrees freedom existing the prototype can duplicated. 

Several significant conditions have been revealed the tests embraced 
the two programs, and, although full correlation between section and full model 
has not yet been evaluated, certain basic facts there full agreement. 
With respect the original Tacoma Narrows Bridge the most pertinent 
these facts may outlined follows: 


(a) linear relationship exists between wind velocity V., corresponding 
the velocity which motion starts given mode, and the frequency the 
several vertical modes. 

(b) Each vertical mode appears twice but the same frequency. The 
V./(N for the second appearance twice that the first. From this 
follows that for given configuration the critical velocity for the second 
appearance any mode twice that associated with the first appearance. 

(c) For any given shape the suspended structure and angle attack 
the wind, the value constant for all vertical modes their first 
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appearance and constant twice that value for all vertical modes their 
second appearance. 

(d) The first appearance each vertical mode very weak and difficult 
isolate the full model. These first appearances are readily observed the 
section model. 

(e) Each configuration and angle attack the wind has its own constant 
value b). 

(f) The torsional modes also appear twice and the same frequency. The 
first appearance weak and damps out with increasing velocity, but the second 
appearance which occurs almost immediately following the damping the first 
becomes catastrophic with increasing wind velocity. general, the critical 
velocity for the first appearance torsional mode half that the second 
appearance. Some the tests involving alterations configuration suggest 
that this relationship may not always hold. 


Oscillation Coefficients.—In aeronautical literature associated with the wind 
flutter problem the dimensionless ratio V/(N sometimes called the 
duced This ratio also referred “reduced wave length,” 
the wave occurring the wake. Some investigators have found convenient 
refer the reciprocal this ratio “reduced frequency.” 

Since the ratio dimensionless and may viewed the ratio two 
velocities, two distances, two times, depending the point view the 
individual the requirements the problem, seems poor practice limit 
the meaning the label making too specific. proposed that the 
developing field suspension bridge aerodynamics the V/(N called 
the parallel practice has already been adopted 
Theodore Theodorsen when refers the critical ‘‘flutter 

from the results tests under several radically different conditions. view 
this evidence, Table should not accepted valid since involves 
variation V./(N over range 2.00 0.25. Laboratory results indicate 
that this value constant the vicinity 2.20 for the second appearance 
the vertical modes and the vicinity 1.09 for the first appearance. 

The linearity critical velocity with frequency completely confirmed 
section model tests the California Institute Technology (unpublished 
report), which course confirm well the constant value for 
given configuration. 

Example the author’s check for high velocity good, but the minimum 
velocity check error factor least five. The writer has not found 
the time for critical examination into the possibility making this mathe- 
matical approach fit the known facts more exactly. hoped that the 
author may able so. 

the tests reported Figs. and had been extended higher wind 
velocities and means had been provided prevent the development the 
torsional mode, successive higher vertical modes would have developed. (In 


Mechanism Flutter, Theoretical and Experimental Investigation the Flutter 


Theodore Theodorsen and Garrick, Report No. 685, National Advisory Committee for Aero- 
nautics, 1940. 
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Fig. 16, frequency cycles per second.) Following this procedure, mo- 
tions involving eleven nodes have been noted with indication that higher 
modes were not possible. The highest mode involved frequencies excess 
360 cycles per min (since all field frequencies must multiplied the square 
root the scale factor this would correspond approximately cycles per 
min the field). These high frequencies were especially severe the model. 
Therefore, they were never maintained for any long period time and thus 
were never open careful analysis. 

The writer feels that there evidence substantiate statement (a) (see 
heading Part “Conclusions Relative Vertical Aerodynamic Oscilla- 
every indication that vertical modes could have developed 
well above the eight-noded mode observed the field miles per hr. 
assumed 0.85 cycles per sec for eleven nodes the prototype and 
per sec miles per hr. 

view the accumulated evidence, the writer cannot agree with state- 
ment (b) either, since the evidence from both section and full model tests 
indicates that torsional tendencies are present this structure from the time 
the velocity reaches about miles per hr. 

should also emphasized that torsional motions can develop sec- 
tion having positive slope the torque graph. Fig. shows torque graphs 
for large model the original bridge which was tested both the California 
Institute Technology and the University Washington. These curves show 
that full circular fairing applied the leading and trailing girders reversed 
the slope the torque graph both tests; yet when the full model was tested 
this form the amplitude response curves were similar those Fig. 16. 
Recent section model tests the University Washington sections involv- 
ing trusses have shown decided torsional motion under certain small angles 
attack although the torque graph showed definite positive slope. 

The chief differences between the faired and the unfaired sections seem 
that the faired section the b)-value for the vertical modes con- 
siderably lower, and the maximum amplitudes developed the various modes 
somewhat lower, than the unfaired bridge. The degree torsional insta- 
bility was somewhat altered, and there reason believe that the second 
appearance torsion may have damped out again after reaching angle 
about 10° and wind velocity approximately miles per hr. These tests 
were run just before the model was dismantled. They were rather hurried and 
were never repeated for confirmation. Therefore, they should taken 
indications only. The test was not extended far enough determine whether 
torsion appeared again higher wind velocities. 

attempt has been made recently explain the vertical motions observed 
the basis the static wind tunnel tests—that is, the basis negative 
aerodynamic damping. Experiment has demonstrated that only very slight 
differences showed vertical modes between section that had definitely 
negative slope the lift line (say, d/b-ratio 0.35) and the section the 
original bridge d/b-ratio 0.2) that showed positive slope. Thus other 


7 


June, 1944 FARQUHARSON AERODYNAMIC STABILITY 1007 


Second 


per 


Cycles 
w 


Frequency, 


Double Torsional Amplitude, Degrees 
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Double Vertical Amplitude, Inches 


LEGEND 
UNLOADED BACKSTAY 
Open Deck,20% Light 
Open Deck,Weight Included 65.5% 
Closed Deck,Weight Normal 
Closed Deck,Weight Included 50% 
COMPLETE BRIDGE 
Open Deck,20% Light 
Closed Deck,Weight Normal 
SPRING REPLACING SIDE SPANS 
Open Deck,20% Light 


Closed Deck 1.09 


Open Deck 1.86 


Critical Velocity, Feet per Second 


OSCILLATIONS 


2.4 
Complete Bridge 
2.0 20% Light 
Open Deck 


1-NV Denotes One-Noded 
Vertical Motion 


1-NT Denotes One-Noded 


Velocity, Feet per Second 
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possible sources have been investigated anticipating more satisfactory 
explanation. 

Early studies had led the setting aside the discharge” theory; 
but subsequent investigations have indicated that this theory tenable 
connection with the excitation vertical modes the basis 
is, has been demonstrated that excitation given mode 
can occur under starting conditions involving vortex discharge rate that 
shows some multiple the frequency the resulting oscillation. 


> 


LEGEND 


-10 
Values 


Numerous investigations the flow were conducted using hot-wire 
conventional form and recording the resultant velocity fluc- 
tuations through the use piezoelectric oscillograph. Typical 
obtained this manner are shown Fig. 18. 


Hill Book Co., Inc., New York, Y., 1934. 


47a Correction for Transactions: the six subcaptions Fig. the expression etc., should 
read 1.52 per etc. 
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result these the hot-wire tests made simultaneously the University 
Washington and the California Institute Technology (the latter used 
different means recording), now felt that reasonable explanation 
the phenomenon available. 


(a) MODEL FIXED; SUBHARMONIC EXCITATION, 93%; 


(b) SUBHARMONIC EXCITATION STILL PRESENT; DOUBLE AMPLITUDE =0.6"; 


(c) SUBHARMONIC EXCITATION STILL PRESENT; DOUBLE AMPLITUDE=1.0"; 


MOTION MODEL CONTROLS VORTICITY; DOUBLE 


DOUBLE 


(f) DOUBLE 


the wind strikes the section that rest, the rate vortex discharge 
from the leading girder will indicated multiple some natural frequency 
the structure, and the structure will begin oscillate appropriate mode. 
the velocity maintained value close the critical velocity for that 
mode, the amplitude will remain small and excitation will continue sub- 
harmonic basis. the velocity increased, however, the section begins 
take charge the flow, and subharmonic excitation yields slowly harmonic 
excitation.* some velocity, generally well below that associated with the 
maximum steady state, the excitation becomes harmonic and the vorticity 
completely phase with the motion. words, the section has taken 
complete charge the flow. velocity increased still further, the section 


Experimental Results Wing Flutter,” William Bollay and Brown, Journal the 
Aeronautical Sciences, June, 1941, pp. 313-318. 
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begins lose its domination over the flow, the amplitude falls off, and the 
section approaches rest far that particular mode concerned. 
general, there exact rest point the amplitude response picture for any 
configuration since there certain overlapping the modes. 

The domination the flow the section seems somewhat tenuous even 
the maximum steady state since frequently noted that amplitude will 
decrease markedly without any change velocity only increase constant 
velocity maintained. the case the lower vertical modes especially, 
has often been difficult build the maximum amplitude under what appeared 
exactly similar conditions. For example, the response curves for the 
fundamental (0-N and single-noded (1-N modes recorded Fig. 
show amplitudes which were only noted two occasions. During all other 
attempts the maximum steady amplitude indicated considerably lower value. 
general, the higher vertical modes could built maximum amplitude 
with little difficulty. 

Fig. 18(d), the motion the model full control the flow, and the 
maximum double amplitude has reached 1.4 in. (70 in. the prototype). 
the same run and the same velocity few seconds later, the amplitude has 
fallen off shown Figs. 18(e) and 18(f), and the record indicates that the 
control the flow the section not complete—that is, some excitation 
the subharmonic level present again. 

Only the basis the subharmonic theory excitation can the double 
appearance each these modes explained. There some evidence, not 
yet well substantiated, that the first appearance harmonically excited, 
whereas the second appearance the excitation subharmonic. 

Part 4.—The author’s general statements under the heading, 
Aerodynamic are not full accord with experimental evidence 
derived from the full model which was tested with many variations the con- 
figuration the suspended structure. These configurations are listed 
Table 

Lift, drag, and moment tests fixed sections the wind tunnel are avail- 
able for all these configurations with the exception sections and 11, 
Table significant note that only section did negative slope 
develop the lift graph. Likewise, all these sections showed both vertical 
and torsional motions with the exception section which torsional 
motion developed velocities value greater than 110 miles per 
(prototype velocity). With great difficulty and only after very long runs, slight 
motion was observed the two-noded and three-noded vertical modes 
section 10. The maximum amplitude developed the end approximately 
min (212 min, prototype) was 0.14-in. double amplitude the two-noded mode. 

Section developed the most vigorous motion (see Fig. 16) both vertical 
and torsional modes, and significant note that the slope its lift graph 
was strongly positive all angles and that its moment graph was negative 
slope over the middle range. All tests with girders any depth showed nega- 
tive slopes the torque graph, and the only section that showed torsional 
stability (section 10) had positive slope the torque graph. respect 
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the torque graph for girder bridges, the laboratory tests are complete agree- 
ment with the author’s claims. 

However, recent tests section models and full models new design 
using truss for stiffening show positive slope the torque graph between 


TABLE GIRDER SECTIONS 


No. Section Description 


Roadway slab completely removed leaving stringers and 


to 


Section with opening covered screen proper mesh 
Solid roadway slab with sidewalks and curbs removed 


Both sidewalks and roadway slab removed, stringers still 
place 


section 
Full circular fairing both girders, open deck 
Full circular fairing both girders, closed deck 
section 


the +12° and the angles attack; yet the model exhibits strong torsional 
instability when the angle attack approximately +7° greater and 
weaker manifestation torsional instability when the angle attack —7° 
greater. This positive slope the torque graph and resultant instability 
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the section were likewise demonstrated fixed model and section model tests 
California Institute Technology (unpublished report). 

Use Center stays are only effective inhibiting asym- 
metric torsional modes, and although they have been installed number 
bridges was never with this mind. They were first conceived 
means reducing the longitudinal pendulum-like motion that observed 
connection with asymmetric vertical modes, and there evidence that they 
have ever been designed with regard the stress that might arise either the 
vertical the torsional modes. There some evidence indicate that the 
center stays the Tacoma Bridge failed (or yielded) strong eight-noded 
vertical motion which was active for some time before the transition the first 
asymmetric mode torsion. possible that this yielding during vertical 
motion permitted the development torsion with its accompanying much 
higher loading the stays, which caused further yielding and ultimate failure 
(possibly accompanied failure the lateral system midspan), leaving the 
structure quite unprotected against full development catastrophic torsion. 

Open-Floor writer not optimistic the author 
regarding the beneficial effect open-floor construction. Tests the full 
model have been confined the sections listed Table but significant 
that Messrs. Prandtl and showed strong three-noded vertical mode 
with amplitude great that shown tests Messrs. Theodorsen and 
(the lower modes were very weak) and with single-noded torsional 
motion excess 15° amplitude. much higher velocity was required 
bring each these modes but their amplitude was disquieting. 

was confidently expected that the complete removal sidewalks and 
roadway slab (section Table would vastly improve the stability condition; 
but tests involving two complete runs showed three-noded and four-noded 
vertical modes large amplitude and also developed three-noded torsional 
mode strong character. The three-noded torsion was especially violent and 
reached amplitude excess 15°. 

writer complete agreement with the author’s con- 
cluding statement that the most scientific approach this problem calls for 
the elimination the cause the motion. 
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